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EXECUTIVE SUMHIRI 



QBJECTI¥.E 

The complexity of Eulerian long range transport models 
presently being developed in Canada and the U.S. has resulted 
in questions about applying such large computer models to 
derive seasonal and annual sulphate or nitrate deposition. 
Comprehensive Eulerian models are designed tc simulate the 
physical chemical and transport processes of individual 
pollution episodes in detail. These models require very 
detailed meteorolcgical input data. The processing of that 
meteorolcgical data and the subsequent running cf the Eulerian 
model requires significant a:r.ounts cf time en large, fast 
computers . 

Several methods have been suggested for obtaining long-terrr. 
deposition results for Eulerian models. These proposals 
include: 

1« A brute force approach where the comprehensive Eulerian' 
model is run with a year long meteorological data set. 
This would require substantial computer and manpower 
r e s o^ u r c e s . 

Mm Development of simplified "engineering models" which 

would be calibrated against the comprehensive Eulerian 
model for a number of individual pollution episodes. The 
engineering models would then be used to derive annual 
and seasonal deposition. 

3. Development of a climatology of precipitation episodes 

for a given geographic region with subsequent derivation 

of physi cal /sta t i s t ical relationships between the 



precipitation climatology and observed precipitation 
chemistry. The comprehensive Eulerian model would be run 
for a number of representative episodes during a year. 
Seasonal and annual deposition totals would then be 
derived by using the cliraatological data for the year to 
extrapolate from the limited number of model runs. 

The objective of this report was to investigate the feas- 
ibility of the third method for deriving long-term deposition. 



APPROACH 

A compreiiensive climatology of daily precipitation events for 
four regions of Ontario for 1982 and 1983 was produced. This 
data set included 16 discrete and 18 continuous precipitation, 
climatology descriptors derived from surface and 850 mbar 
weather maps and trajectory plots, regionally representative 
hourly surface meteorological data, nearby upper air soundings 
and precipitation data from hourly and daily cl imat ologi cal 
records. Using 1982 as the base year, statistical comparisons 
were made to the APIOS (Acid Precipitation in Ontario Study) 
observed daily precipitation chemistry data. For each of the 
four regions the best cl imatologi cal parameters for 
representing precipitation chemistry within that region were 
selected and statistical models for estimating S0^ and NO3 
precipitation concentrations were produced. 

These statistical models were then used to classify the 
precipitation periods in each region into eight categories 
ranging low to high SO^ or NO3 concentrations. The category 
ranges were selected to give approximately the same number of 
precipitation events in each category. 



The final step in the feasibility study documented in this 
report was to compare predicted annual and seasonal SO^/NCj 



deposition at individual stations with the observed deposition 
at that station. The predicted values make use of 85 selected 
days during the year (from 3 to 11 days in each month of the 
year). Data for these days were used to produce precipitation 
weighted averaged concentrations for each cl imat ol ogi cal cate- 
gory. Using the climatological classification of each event 
in the year along with the above average S0^/N03 concentra- 
tions for the eight categories, seasonal SO^/NOj depositions 
were estimated. These estimated deposition results were then 
compared with the observed deposition to see how useful the 
climatological stratification was. 



CONCLUSIONS AND FUTURE WORK 

The principal conclusions of this report can be summarized as 

f oil ows I 

lis Although many of the meteorological parameters were cor- 
related with either SC^ or II 3 precipitation chemistry, 
the correlation coefficients were usually less than 0.^. 
Also a climatological descriptor that was significant foi 
SO^ in a particular region would not necessarily be 
important for NO 3 in that same region^' 



The statistical models incorporating the best two dis- 
crete and the best two continuous meteorological descrip 
tors gave correlation coefficients for each region which 
ranged from Q.^4 to 0.83 for SO., in warm months and 0.58 
to 0.74 for S0^ during cold months. Corresponding 
figures for NOj were 0.55 to 0.82 during warm months and 
0.56 to 0.70 during cold months. The statistical models 
had the most difficulty in southwestern Ontario where 
average concentrations were the highest. 



Analyses indicated that the best el imat ologi cal variables 
for SO^ were not necessarily the best for NOg. Also, 
better statistical correlations were found when each 
region was examined separately when determining the best 
set of variables. In some cases such differences might 
be related to the relationship of the receptor region to 
source regions. In other cases it might sifiply be a 
matter that many of the meteorological variables were 
inter-related with one slightly better for a particular 
data set. A list of the best overall predictors of pre- 
cipitation chenistry was produced. 

Comparisons of seasonal or annual deposition extrapolated 
from a limited data set versus actual observations showed 
good overall agreement. Extrapolated seasonal 
depositions were derived in two ways in the report. In 
one set of cal caula t i ons the data was stratified into 
eight climatic categories before doing the extrapolat- 
ions. The objective was to explain as much of the 
variance as pcssible with the stra t If i ea t i on and thus 
Improve the overall prediction. In the second set of 
calculations no data stratification was performed so the 
deposition predictions were based only on extrapolating 
to the total precipitation in the season. The two 
methods for predicting seasonal deposition performed 
equally well . 



The main factor affecting the accuracy of the extra- 
polated deposition was the percentage of the precipita- 
tion cases and of the precipitation amount included in 
the limited data set. Seasonal deposition predictions 

for overall regions were usually within 5 to 151 of the 
observed values. For individual stations the worst 
results were produced for stations where the limited data 
set most poorly sampled the precipitation events at that 

station . 



As the region for which you are attempting to estimate 
deposition gets larger the selection of modelling days 
becomes more difficult. This means that any method of 
extrapolating from limited model runs to seasonal 
depositions can only be expected to produce reasonable 
results over a limited area. 

The lack of improved performance when climatological 
classification was included when doing the extrapolations 
could be caused by the small number of events included in 
each category (i.e. eight categories wer^e used). As 
futur-e wor-k, this stage of the exercise will be r'epeated 
using different climatological groupings* 

As an addendum to this report analyses were carried out with 
1983 data. The predictions were based on 71 days of cata. 
Overall the err cr were 10 to 15* and again the performance 
was net improved by climatological classification. 



ARB-3L0XAM 
EXEC-SUM. 1 



SOMMAI RE 
OBJECT IF 

La complexite des modeles euleriens de transport a grande distance 
actuelleiient elabores au Canada et aux Etats-Unis a suscite des 
questions sur la fagon de mettre en application ees grands modeles 
informatiques pour calculer les depots saisonniers et annuels de 
sulfate ou de nitrate. Les modeles euleriens complets permettent 
une simulation detaillee des processus physico-cliimiques et des 
processus de transport lies a chaque eas de pollution. Ces inodeles 
necessitent I'obtention de donnees meteorologiques tres precises. 
Pour traiter ces donnees meteorologiques et executer par la suite le 
modele eulerieni il faut disposer de gros ordinateurs rapides et de 
beaucoup de temps. 

On a propose plusieurs methodes pour obtenir des resultats 
concernant les depots a long terrae a partir de ces modeles. Ces 
propositions sont les suivantes ; 

I* Une methode brute prevoyant 1 'execution du modele eulerien 
complet a partir de 1' ensemble des donnees meteorologiques 
d ' une annee. Ce procede necessiterait des ressources 
inf orraatiques et humaines importantes. 

1;» La mise au point de modeles techniques simplifies qui seraient 
etalonnes d'apres le modele eulerien complet et appliques a un 



certain nombre de cas de pollution pris separement . Ces modeles 
serviraient ensuite a calculer les depots annuels et 
saisonniers . 

3.. L 'elaboration d ' une climatologie des cas de precipitations pour 

une region geographique donnee permettant par la suite de trouver 
des liens physiques ou statistiques entre 1 'aspect climatologique 
des precipitations et 1 'aspect chimique des precipitations 
observees. Le modele eulerien complet serait execute pour un 
certain nombre de cas representatif s durant 1 ' annee . Le total 
des depots saisonniers et annuels serait alors obtenu a I'aide 
des donnees climatologiques de 1 'annee en question par 
extrapolation a partir du nombre restreint d 'executions du 
mod e 1 e . 

Le present rapport examine la possibilite d 'employer la troisieme 
methode pour calculer les depots a long terme. 

METHODE 

On a produit une climatologie detaillee des precipitations 
quotidiennes de 1982 et 1983 pour quatre regions de 1 'Ontario. Get 
enseiible de donnees comprend 16 d^escripteurs climatologiquei 
discrets et 18 descripteurs continus obtenus a partir des cartes 
meteorologiques en surface et a 850 millibars et des graphiques des 
trajectoires, des donnees ffliteorologiques horaires des observations en 



surface reprasentatives des regions * des sondages aerologiquas et das 
donnees sur les precipitations tirees de releves climatologiques 
lioraires et quotidiens. En prenant 1982 comme annee de reference, on 
a procede a des comparaisons statistiques avec las donnees chimiques 
sur les precipitations quotidiennes figurant dans 1 'Etude sur les 
precipitations acides en Ontario, Pour chacune des qoatre regions, on 
a releye les parametres climatologiques representant le mieux la 
coroposition chimique des precipitations et on a produit des modeles 
statistiques permettant d'estimer les concentrations de SO4 et de 
NO3 . 

Ces modeles statistiques ont ensuite servi a classer les cas de 
precipitations de chaque region en huit categories, salon la 
concentration de SO4 ou de MO3 . Les valeurs limit es des 
categories ont iti cliolsias afin qu'll y ait environ le mime nombre 
de cas de precipitations dans chaque categorie. 

La derniere etape de 1 'etude, consignee dans la rapport, consistait 
a comparer les depots annuals at saisonniers de SO4 et de NO3 
prevus a chacune des stations avec les depots observes a ces mimes 
stations. Le s ch i £ f r e s pr e v u s por t ai en t sur 86 j ours cho i s i s dur ant 
1' annee (de 3 a 11 jours pour chaque mois de 1' annee). Les donnees 
concernant ces journees ont permis d'etablir las concentrations 
moyennes ponderees pour chaque categorie climatologiqua. On a pui 
ainsi astimar las depots saisonniers de SO4 et de NO3 d'apres la 



classification climatologique de chaque cai de precipitations et 
d'apres les donnees sor les concentrations moyennes de SO4 et de 
NO3 pour les huit categories. Par la suite, on a compare ces 
estimations aux depots observes poar determiner I'utilite de la 
strat i f icat ion cl imatologique . 



CONCLUSIONS ET TRAVAIL A ACCOMPLIR 

Voici les principales conclusions du rapport s 

I* Bien qu'un grand nombre de parametres meteorologiques aient eti 
mis en correlation avec les concentrations de SO4 ou de 
NO31 les coefficients de correlation etaient habituellement 
infer ieurs a 0,4. De plus, un descripteur climatologique qui 
etait important pour les concentrations de SO4 dans une 
region donnee ne 1 'etait pas necessairement pour les 
concentrations de NO3 dans la meme region. 

Mm Leg modeles statistiques comprenant les deux meilleurs 

descrlpteurs meteorologiques discrets et les deux meilleurs 
descripteurs continus ont donne, pour chaque region, des 
coefficients de correlation allant de 0,44 a 0,83 durant la 
saison chaude et de 0,58 a 0,74 durant la saison froide, en ce 
qui a trait au SO4. Pour le NO3, les coefficients 



correspondants variaient de 0,56 a 0,82 durant la saison chaude 
et de 0,56 a 0,70 durant la saison froide. C'est dans le 
sud-ouest de 1 'Ontario, ou I'on retrouve les concentrations 
moyennes les plus elevees, qu'il a eti le plus difficile 
d'etablir des modeles statistiques. 

i« Selon les analyses, les variables climatologiques les plus 

appropriees pour le SO4 ne le sont pas necessaireiient pour le 
NO3 . En outre, on a pu etablir de Beilleures correlations 
statistiques en examinant chaque region separement pour 
deteriiiner les variables les plus pertinentes. Dans certains 
cas, les differences pouvaient etre reliees aux rapports entre 
la region receptrice et les regions emettrices. Dans d'autres 
cas, les differences pouvaient etre dues au fait qu'un grand 
nombre de variables meteorologiques etaient reliees entre elles 
et que I'une d' entre elles etait un peu plus appropriee pour un 
ensemble particulier de donnees. On a alors produit une liste 
des variables servant le mieux a la prevision de la conposition 
cliimique des precipitations, 

4i II y a eu, dans 1' ensemble, une bonne concordance entre, d'une 
part, les depots saisonniers ou annuels extrapoles d'apres un 
ensemble de donnees limite et, d' autre part, les observations 
reelles. L' extrapolation des depots saisonniers a ete obtenue de 

deux fagons dans le rapport. Dans une premiere serie de calculs, 
on a classe les donnees selon huit categories climatiques avant 
d'effectuer les extrapolations. Cette methode visait a expliquer 



le mieux possible les ecarts grace a la stratification et a 
ameliorer ainsi la prevision globale. Dans la deuxieme serie de 
calculs, aucune stratification n*a ete effectuee; les previsions 
d© depots etaient uniquement fondees sur 1' extrapolation des 
donnees portant sur les precipitations totales dans la saison. 
Ces methodes permettent toutes deux une bonne prevision des 
depots saisonniers. 

Le pourcentage des cas et de la hauteur de precipitations compris 
dans 1* ensemble de donnees limite const itue le principal facteur 
influant sur 1 'exactitude des depots extrapoles. Les previsions 
de depots saisonniers variaient habit uellement de 5 a 15 pour 100 
par rapport aux depots observes. Dans chacune des stations, les 
pires resultats sont apparus dans les cas ou les precipitations 
etaient le plus mal representes dans 1' ensemble de donnees 
lifflitea 

Plus la region dans laquelle on essaie de prevoir les depots 
s'elargitf plus il devient difficile de choisir des journees 
types. Cela signifie que toute methode visant a extrapoleri a 
partir de donnees limitees, des previsions de depots saisonniers 
peut vraisemblablement produire des resultats valables uniquement 
s ' i 1 s'agit de secteurs restreints. 



Le petit nombre de cas inclus dans cliaque categorie (on a utilise 
huit categories) pourrait expliquer le fait que las resultats 
n'aient pas etl meilleurs lorsqu'on a ajoute la classification 
climatologique a 1 'extrapolation. Dans un prochain travail, on 
recommencera cette etape en se fondant sur des groupes 
climatologiques dif f irents» 

Vous trouverez en supplement des analyses menees a partir de donnees 
de 1983. Les previsions sont fondees sur les donnees de 71 jours. 
Dans 1 'ensemble, le taux d'erreur varie de 10 a 15 pour 100; la 
encore, la classification climatologique n'a pas permis d'obtenir da 
mailleurs resultats. 
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1, Introductloii 



Thtre are in existence models for the long range transport, trans formation, and 
deposition of airborne pollutants which are necessarily complex and expensive to 
run. If they were to be applied for the daterinination of seasonal and/or annual 
averages of pollutant deposition amounts, the continued operation of these models 
over the time period of the averaging calculations would be prohibitively 
expensive. There is, however, and alternative to this brute force approach. If 
there was a elleiatology of precipitation episodes for one year or longer. It may 
be possible to relate a limited set of complex model runs on an episodic basis to 
the precipitation climatology for the same geographical area. This could 
Goncelvably permit the extension of the results from the limited set of model runs 
to a much longer time frame by means of physical/statistical relationships to the 
overall weather statistics. 

With a view towards accomplishing the objective described above in principle, this 

study was carried out to test the feasibility of using meteorological parameteri 

to produce SO, and NO. deposition statistics In Ontario. 
'^ 4 3 - 

The basic approach had several aspect to the work. It was first necessary to 
divide the province into approximately homogeneous regions according to the 
meteorology and the distribution of pollutant deposition stations within the Acid 
Precipitation in Ontario Study (APIOS) network. The base year was then selected 
for meteorological analysis and comparison to the APIOS measurement data for SO, 
and NO, during the same time period. Statistical analyses were carried out on 
selected meteorological parameters and deposition data for the possible 
development of a model which could predict deposition values on the basis of the 
selected storm parameters. A second year's analysis could then be carried out and 
the results from the first year's study could be tested on the second year's 
analyses as an independent data set. 



2i Meteorological Analysis 

The base year aelected for study was 1982, with 1983 to be used as a test year for 
the results for the first year of study. Following this decision, it was 
necessary to: i) develop storm classification schemes to be used for analysis, 
11) acquire the appropriate meteorological data for review, 111) analyze the 
available data according to selected parameters and develop the analysis database, 
iv) perform tiumerical and statistical experiments to determine relationships (if 
any) between the meteorological and deposition data, and v) generate deposition 
statistics from meteorological data using the relationships found in iv) , 

A survey was carried out on the various methodologies available. Based on a 
ilterature review as well as discussions with other meteorologists and physical 
scientists, a set of procedures for use in the storm analysis and classification 
phase of this study was developed. These are reproduced below In 9 specific 
steps . 

%m Fix boundaries for all regions. 

2. From the APIOS data set for 1982, determine the days In each region when 
precipitation fell and measurements of sulphate and nitrate 
concentrations were recorded and deemed reliable. 

3. Determine the meteorological stations which lie In each region. There 
were 4 groupings of surface stations per region to reflect the 4 
different types of observations considered. They were: 

- hourly and synoptic observations 

- recording raingauge data 

- daily precipitation totals 

- snowfall amounts as measured with a Nipher gauge 

4» Acquire hourly data froni apfropriate stations for the days of interest. 

1* Froii the data in Step 4, define the hours which comprise the 
precipitation event for each region for each day of interest. The 
average temperature, dewpolnt, wind, and predominant clouds and weather 
for each precipitation event was determined, 

I, From synoptic surface and upper air maps, characterize the type of 
synoptic system, the track of the controlling low pressure system, and 
the airmassCes) involved. 

7, From the trajectory analysis, determine the sector and distance of the 
air parcel trajectory 24 and 48 hours before the precipitation event. 



8, Determine representative upper air souEdlngs for each precipitation 
event and calculate precipitable water, a stability Index, and analyze 
moisture profiles to get cloud base height, cloud top, and cloud 
thickness. 

9. For each region, obtain the following precipitation data. Analyze the 
recording rain gauges to get average rainfall rates, maximum rainfall 
rates, and average duration of precipitation events for each day of 
interest. From the dally precipitation amounts, get the average dally 
precipitation amount and the percentage number of stations receiving 
precipitation. From the snowfall data (stations with Nlpher gauges), 
determine the average daily snowfall and snowfall rate (snow amount per 
6 hours at the synoptic station). 

The province of Ontario was divided into four regions (Southwestern, Central, 
Southeastern, and Northwestern) for the purpose of this study. A summary of the 
APIOS, hourly, synoptic, and climatological stations within each of the study 
regions Is given in Appendix A. Appendix B lists the 25 different meteorological 
parameters abstracted from the data for the statistical experiments. It was 
necessary to resolve which stations could be considered representative of the 
APIOS station (if rainfall is not measured close by), or alternatively, 
representative of the region. This could only be done after examination of the 
data. 



3. Cllmatologlcal Data Analysis 



This section of the report deals with the specifics of the data analysis and 
the statistical results subsequently carried out. 



3.1. Synoptic Analysis 

In carrying out the synoptic analysis of the meteorological data, the 
following procedures were followed. 

1) All significant hourly weather reports for 1982 for the AES 
stations in the 4 regions were printed out using the MOE IBM 
computer, 

ii) Using this printout and the Ontario Weather Centre surface analyses 
(every 6 hrs), the U.S. Weather Bureau dally surface analysis and 
500 mb map (once per day at 1200Z) , and the CMC 850 mb analyses (00 
and 12Z) the following was obtained for each region for each day of 

interest: 

• synoptic type (as in the planning report with the addition of 
a 7th type, namely a warm and cold front) 

- track of the controlling low (as In the planning report with 
the addition of a Hudson Bay low category) 

- predominant alrmasses (up to 2 of the 5 specified) 

- the time of the precipitation event and the most 
representative hourly station 

- the most representative upper air station And the 
corresponding time of observation 

ill) From the representative times determined from the previous step, 
surface and 850 mb trajectories were analyzed for the following: 

- the sector from which the air originated 24 and 48 hr previous 
(octants numbered counterclockwise beginning In the NNW) 

- the range of the air parcel 24 and 48 hrs previous (divided 
into 200 km segments) 



3.2. Surface Analysis 

1) The data from the most representative hourly station for the time 
of the precipitation event, as previously determined, was displayed 
using the IBM AT with software specifically developed for this 
purpose (program MDASHR) . 



11) Subjective Judgments were then made to try to characterize the 
precipitation event In terms of the nature of the precipitation 
(rain, snow» freezing precipitation and mixed precipitation) « the 
type of precipitation (continuous or convectlve) , and the 
predominant cloud type (towering cumulus, altocumulus, altostratus, 
cumulonimbus, cumulus, nlmbostratus ,stratocumulus , or stratus) with 
its reported base . The judgments focussed on the reports which 
indicated that the most significant precipitation was being 
observed. 

ill) The time of the most significant preclpttatidn was reassessed and 
for that time period, the average surface temperature, dewpoint, 
and wind were calculated. 



3.3. Upper Air 



I) Data from the most representative upper air station, as previously 
determined, was displayed using the IBM AT and program MDASUP (only 
standard levels were available) . 

II) The 850 mb temperature and wind were selected, preclpitable water 
(qAp/g) and totals index (tg5Q-t ^t^g^^-t^QQ+t^ ) were 
calculated where q Is the average mixing ratio in Ap, p is 
pressure, g is gravity, and t and t. are the temperature and 
dewpoint at the Indicated pressure levels in mb) * 

ill) Cloud top was subjectively determined from an examination of the 
temperature-dewpolnt spreads where a general rule of t-t, of less 
than 4'*C indicated the presence of cloud. 



3.4. Rainfall 



I) The 4 regions were grldded by latitude and longitude and the 
climatologlcal stations assigned their appropriate grid within the 
region. The Southwest region was divided into 10 subregions, the 
Central Into 12 subregions and the Southeast and Northwest regions 
Into 7 and 6 subregions respectively (see Figures 3.1 and 3.2). 

II) The subregions were weighted according to their areal extent. A 
weighting of 1,0 indicated that the subreglon was a degree of 
latitude by a degree of longitude in dimensions. A weighting of 
less than one Indicated that part of the degree of latitude by 
degree of longitude area fell outside the region In question (e.g. 
over water). Table 3.1 gives the weightings and the number of 
hourly and daily reporting rainfall stations in each of the 
subregions 

III) Eight precipitation fields were examined, namely maximum rainfall 
rate in mm/15 min, mm/30 mln, mm/hr, duration of precipitation 
(hrs), average hourly rate of the precipitation, dally amount, 
maximum 6 hr amount, and the proportion of stations receiving 
rainfall. For each field for each region for each day the averages 




Figure 3.1 NorthwesteTn region of the APIOS 
study area. 
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Figure 3.2 Ceitral, Soutihwestem, and 
Soufteasteni regions, of the 
APIOS study area. 
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Table 3.1: Weight Ings for the subreglons of the 4 main regions and the number of 
hourly (H) and daily (D) rainfall stations in each subreglon. 



for the subregions were first determined. Then using the welghtlnga 
for the subregions, the regional average was obtained. In doing so 
standard deviations of the various precipitation fields by 
subregion and by region were also obtained. 



3,5, Transfer to HOE IBM Computer 

All the various fields were merged into 4 data sets (1 per region) and 
transferred to the IBM Computer In Downsvlew. The record length was 173 
characters and the list of variables and the format was the following s 

Ix 

month 12 

day 13 

synoptic type i3 

track of low 12 

airmass(es) l3C,a6 

24 hr sector of 850 mb trajectory lx,i3 

24 hr range of 850 mb trajectory 13 

48 hr sector of 850 mb trajectory 13 

48 hr range of 850 mb trajectory 13 

24 hr sector of surface trajectory 13 

24 hr range of surface trajectory 13 

48 hr sector of surface trajectory 13 
48 hr range of surface trajectory * 13 

type of precipitation 12 

nature of precipitation 12 

predominant cloud type 14 

height of cloud base (m) f6.0 

surface temperature (^C) f6,l 

dewpoint (°C) f6.1 

surface wind direction (deg) fS.O 

surface wind speed (m/s) f5.1 

850 mb temperature ("C) f6.1 

850 mb wind direction (deg) f5.0 

850 mb wind speed (kts) fS.O 

cloud top height (m) fS.O 

precipl table water (mm) f7. 2 

total totals Index f7.2 

average maxlmuni rate miD/15 min f7,2 

average maximum rate mm/30 mln f7,2 

average maximum rate mm/hr f7.2 

average duration of precipitation (hr) f7,2 

average hourly rate of precipitation (mm/hr) f7.2 

average precipitation amount (mm) f7.2 

average maximum 6 hr amount (mm) f7.2 
percentage of stations reporting precipitation f7.2 

Note; missing data is indicated by 9*s. 



3.6. Problems iti Data kn&lfsis 

The following Is a list of the most apparent problems that occurred In trying 
to relate the meteorological data to the APIOS data collected. 

i) Imprecise precipitation times between neighbouring APIOS stations 

II) Non-uniform on/off times of charts between APIOS stations and 
between days at the same station. 

III) APIOS observations not always taken dally when precipitation 
observed. 

Iv) A good deal of hourly data was missing. 

v) Spurious values for the hourly data occasionally appeared (values 
not correct but not specified as missing) 

vi) No special weather reports at the hourly stations were available 

vii) The upper air data only contained the standard levels and not the 
significant levels as well. 
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4, Result i and Discussion 



The 33 meteorological variables available for this study can essentially be 
divided into 6 types. The first 3 are categorical variables.: 

- syxioptic situation (synoptic type, track of low) 

- a description of the previous history of the air (airraass type, 
trajectories) 

- a characterization of the clouds and weather (type and nature of 
precipitation, cloud type) 

The next 3 are continuous variables: 

- a description of alrmass characteristics (temperature, 
dewpointjprecipltable water, totals index, wind speeds) 

- a description of the clouds (clouds base, top, and thickness) 

- rainfall characteriitics 

The goal of the analysis is to find which subset of these variables is strongly 
associated with the sulfate and nitrate ineasurements taken In the regions. If a 
sufficiently strong relationship is found, then a statistical model can be 
constructed to predict the sulfate and nitrate measurements, 

4,1 General Results 

A summary of the statistics derived from available data for the four regions is 
given in Appendix C for each of the meteorological parameters abstracted. 

A general comment regarding the differences between the synoptics of the northwest 
region and the other 3 is that the northwest region has many more occluded systems 
and fewer fronts. It Is more affected by the Alberta lows and less influenced by 
the Colorado systems than the other 3 regions. An analysis of variance was 
performed on both the sulfate and nitrate measurements in the four regions to test 
for statistically significant differences. For SO, , the F value (the ratio of the 
between- column variance to the between-row variance in a matrix of SO. readings by 
region) was 29.40, giving a significance probability of 0.0001, I.e. this F value 
would be exceeded only 1 in 10000 times if the SO, means of the 4 regions were in 
fact equal. To further examine which means were statistically significant, 
Duncan's Multiple Range Test (run using the ANOVA routine on SAS User's Guide: 
Statistics 1982 Edition pg 122) indicated that the means which were significantly 
different at a significance level of 0.05 were the northwest region with the other 
3 regions, and the central and southwest regions. 

For the nitrate measurements, the F value was 33.73 Indicating a significance 
value of 0,0001. Duncan's Multiple Range Test In this case showed that the nitrate 
readings in the northwest and the central were each different from the other 3 
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regions while there was no difference between the southeast and southwest 
measurements (at a 0.05 slgnlfleanee level). 



4.2, Southwest legion 

Due to the higher density of rainfall statlotis, closer proximity of upper air 
stations, and the accessibility of a 24 hr reporting hourly station amongst the 
APIOS sites (London), the relationship between the meteorological variables and 
the sulfate and nitrate measurements were examined in great detail for the 
southwest region. It was felt that any meaningful relationships between the SO^ 
and NO. and the meteorological variables would most likely be found in this region 
and greater confidence could be given to these results. 

4.2,1 Seasonal Variation In Sulfate and Nitrate Readings 

The SO, and N0„ averages for 1982 are given by month In Table 4.1, 



MONTH 


SO, 
4 


^°3 


JANUARY 


2.13 


0.75 


FEBRUARY 


3.62 


1,26 


MARCH 


4.04 


1.13 


APRIL 


4.83 


0,78 


MAY 


6.50 


1,20 


JUNE 


5.04 


0.68 


JULY 


5.61 


0.80 


AUGUST 


4.01 


0.51 


SEPTEMBER 


3.45 


0.59 


OCTOBER 


3.48 


0.50 


NOVEMBER 


2.71 


0.55 


DECEMBER 


2.61 


0.48 



Table 4.1: Sulfate and Nitrate readings by month. 



The data were divided into 2 sets (April to September vs the other 6 months) to 
test for a seasonal variation. For SO, , the average for the colder months was 2.99 
mg/1 while for the warmer months it was 4.75 mg/l. To test the statistical 
significance of this an analysis of variance was done (ANOVA routine in SAS User's 
Guide: Statistics 1982 Edition). The F test statistic had a value of 8,95 which 
gave a significance probability of 0.0001. Hence there is compelling evidence that 
the SO, measurements are higher in the summer than in the winter. For N0„ on the 
other hand, the 2 averages are 0.79 mg/1 for the colder months and 0.73 mg/1 for 
the warmer months. The F statistic had a value of 0.45 in this case and a 
significance probability of 0.50 when testing the difference of means. For N0„ 
then there was no evidence of a seasonal variation. 
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4.2.2 Relationships Among the Discrete Variables 

Two way contingency tables were constructed among the discrete nieteoro log leal 
variables to exainine their relative relationships. This was useful In order to 
strengthen our coiifidence in the quality of the Deteorologlcal analysis and as an 
aid in deteriiining an independent set of predictors for the SO, and HO^ 
measuirements. A chl-square test was then performed on each table (SAS User's 
Guldei Basics 1982 Edition FREQ procedure pg 515). The significance probability 
of those tests are given in Table 4.2. The sparseness of data in parts of the 
tables in some cases may invalidate the chl-square test but it Is still felt to be 
a useful start in looking for relationships. 

As expected there was a very high degree of relationship among the meteorological 
variables. For example, It was reasaurlng to see that there Is a strong 
relationship between cloud type and precipitation type. As a broad 
generalizatioii, the trajectory ranges were the least correlated with the synoptic 
variables and clouds and weather variables. This table will be re-examined when 
considering the best variables to Include In the statistical model. 

4.2.3 Correlations Among Continuous Meteorological Variables 

Tables 4.3 a) and b) display the correlation coefficient among the meteorological 
variables. It was obtained by applying the COER procedure in the SAS User's 
Guide: Basics (1982 Edition pg 503). Also displayed are the significance 
probability of the correlation and the number of observations. Again here, there 
was a high correlation among the continuous meteorological variables. At a 
correlation coefficient above 0,5, the following groups of variables were highly 
related: 

- surface temperature, surface dewpoint, precipitable water, totals 
Index, 850 mb temperature, and rainrates of mm/15 mln and mm/30 
min, 

- surface and 850 mb wind apeed. 

- cloud top and cloud thickness. 

- rainfall rates and rainfall amounts. 

- duration of precipitation and proportion of stations receiving 
pracipitation* 

These relations are to be expected and will be reconsidered in light of the 
particular variables which are most highly correlated with the SO, and N0„ data. 

4. 2.4 Relationship Between SO, NO- and the Discrete Meteoroloeical Variables 

In an effort to identify those discrete meteorological variables which are 
important to the study of acid rain the following 2 tables were constructed, 
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SIGNIFICANCE PROBABILITY 
SYNOPTIC TRACK AIRMASS 
TYPE OF LOW 



TYPE OF NATURE CLOUD 
PCPN OF PCPN TYPE- 



TRACK 






0.0001 












AIRMASS 




0.0001 


0.0001 










850 24 


HR 


SECTOR 


0.0001 


0.0001 


0.0001 


0.1118 


0.0035 


0.0001 


850 24 


HR 


RANGE 


0.0001 


0.0472 


0.1543 


0.5670 


0.0811 


0.5881 


850 ■48 


HR 


SECTOR 


0.0001 


0.0001 


0.0001 


0.1981 


0.0343 


0.2349 


850 48 


HR 


RANGE 


0.0323 


0.0001 


0,0008 


0.6106 


0.9356 


0.9487 


SFC 24 


HR 


SECTOR 


0.0001 


0.0001 


0.0001 


0.0085 


0.0158 


0.0270 


SFC 24 


HR 


RANGE 


0.3511 


0.2423 


0,7923 


0.7391 


0.7352 


0.8387 


SFC 48 


HR 


SECTOR 


0.1997 


0.0001 


0.0559 


0.0596 


0.0565 


0,0594 


SFC 48 


HR 


RANGE 


0.2529 


0.0737 


0.5133 


0.8504 


0.0095 


0.4181 


TYPE OF PCPN 


0.0001 


0.3531 


0.0554 








NATURE 


OF 


PCPN 


0.1230 


0.0001 


0,0001 


0.0001 






CLOUD TYPE 


0.0004 


0.0001 


0.0035 


0.0001 


0.0001 




SFC WIND DIR 


0.0001 


0.0001 


0,0007 


0.0550 


0,0040 


0.1000 


850 WIND DIR 


0.0010 


0.0001 


0.0051 


0.1500 


0.0002 


0.0010 



Table 4.2: Significance probabilities among the meteorological variables using 
the chi-square test. 



'm 



CQi(RcL£T!ON COFFFICIINTS / PROS > IftI UNDER H0tRHO«O / NUHBEfl OP OBSERVATIONS 

ChT TAVO IQiVG WAV2 CAT3* CAT33 C*T16 TMR PM TT 

tnl l.C'COOO w.lOi;? C.C§!J42 -O.J'b9«J 0.1<r943 -0.02863 0.07342 -0.13950 0.14407 0.CT346 

Cl: -^ii nT ..CCcC C.1741 0.-.3<,2 O.o?79 &.076S 0.7417 W.3S66 0,12*.4 0.0183 0.3S33 

lit lib Iff 157 1<.3 135 123 122 141 143 

7-/- :.lJ--i7 I.l'wOO r.-iiSai -0.292«0 0.93302 -0.12223 -0.157S0 -0.l322i 0.35673 0.72355 

S-C Ttrr 4.17^1 :-..i:0C C.OOGl 0.0001 0,0001 0.1352 G.07:d 0.14fej O.JOOl 0.0001 

i;i 146 163 1*3 154 146 130 121 152 154 

1.-Vi '-.:5:i.; 0.9BST1 l.TJCOP -0.2<..?'!'a 0.93;'65 -0.13S1V -0.15!27 -0.11831 O.CeSZ} 0.75655 

;<■; J:;WV1 :.<.-t2 C.CCM v.uCOr O.OOOl 0.0001 0.CS3S; 0.0777 C.1962 O.OCOl f.OOOl 

156 J66 13'. 1<S 1S4 146 130 121 152 154 

UiVv -Cvi-^Ss -C.Z*26C' -0. ?*•:•£« l.oriOOC) -0.352el O.305»fc'. -O.OiOt* -O.CO^fS -0.32717 -0.3437» 

irC ''iH^ i^^ ■..6:7i e.ocji ;.oj.oi c.oooo o.oooi o.oooi o.piis 0.9160 o.ooci o.oooi 

1:1 lift Ic,? 159 155 147 131 liJ 153 155 

■:-T.-«. C.:,S«J C-. 93 = 02 C,'?i'<:;- -O.aiSfl 1.00000 -0.183li -3.10*60 -O.C£77£ C. 37130 0.73^23 

.li 7:"r 0,1)70^ O.CCGl t.wOCl O.rOCl O.OOOC 0.021* 0.2100 C,4;<»h 0.0001 O.OOOI 

143 iS-. ii^ 155 U8 l*i !40 122 iS4 Ifb 

CiT^j -'J.0?a&: -:.123?tt -C.!,i:^ 0.^':j;4 -0.13312 1. 00000 -0,0?«13 0.011^3 -0.250*5 -0,19^11 

ii-. *.'»- :,".. C;.7417 0.1332 :.w-:i C.'OOl 0.0/13 O.OOGC 0.5201 C-.r-O*! 0.0014 0.0l2p 

13J :**ft t«& l*-? us IS-i 130 113 156 15a 

CJTl£ :.J7-i2 -:.J5?^: -r.li.527 -0.02044 -0.10660 -0.02013 1,00000 0.^7213 0.01214 0.05147 

cl: T:* 0.3366 0.0728 U.C777 D.'^lfS C.2100 0.S201 C.OOCO 0.0001 0.3503 0.545f 

123 13U 130 131 1^0 130 140 122 Hi 140 

TnH -f. 13553 -C. 13221 -O.ll-il -0.0CJ«5 -0.06775 0.0U33 0.97218 1.00000 0.02775 0.06395 

Cw- u:?Th U.1254 C.!4fil3 C.l<ifc2 C^lfO 0.4*34 0,9052 0.0001 0,0000 0.7535 0.4641 

122 121 121 122 122 |iS 122 122 120 122 

f» C. 14407 :.a>573 D.-ri-SJO -0.2">7r7 0.67190 -0.25065 0.D129*. 0.02775 1.00000 0.32810 

rt.fClr t>^llk :.08d3 COOra o.oooi O.OCOl o.oooi 0.0C16 0.8303 0.7635 0.0000 0.0001 

141 152 I SI 1 93 166 156 133 120 - 165 166 

Tl t.C7's= l,lZiii C.liiiT -0,34379 0.73823 -0,19311 0,05147 0.06395 0.B2S10 1,00000' 

T;,T«L T^TiLs :,3633 J,Cli;C': COOOl O.O&Ol 0,0001 0,0126 0.5459 0,4641 O.OCOl 0,0000 

\*,l 1<H 134 155 168 15B 140 123 165 1*8 
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Table 4.3a, Correlation coefficients, significance probabilities, and number 
of observations among the continuous meteorological variables. 



Table Dl (in Appendix D) gives the mean and standard deviation (SD) of the SO, and 
NO* measurenients as a function of category for the discrete variables. Table El 
(found in Appendix E) displays the distribution of observations among the classes 
is given in frequency and percentage of observations for the top 30% of the SO. 
and N0» readings (above 4.55 mg/1 and 0,90 mg/l respeGtlvely) . 

From Tables Dl and El, soniie inferences regarding the meteorological factors which 
Influence the SO, and N0„ measurements can be drawn. In order for a factor to be 
considered important, the mean of the SO, or NO. readings for that category had to 
be the highest or second highest among the means for all categories of that 
variable. Also the frequency (In percentage) of that factor had to Increase when 
considering the top 307. of the SO. or NO. readings as compared to the frequency in 
the whole data set. 

Mlth the above criteria, the following is a characterization of the high episodes 
of SO, . (Since the factors were considered separately no claim is made that all 
the factors happen simultaneously, although many of the fields mentioned are 
consistent with each other.) 



SO, 



NO, 



SYNOPTIC TYPE 

TRACK OF LOW 

AIRMASS 

850 24 HR SECTOR 

RANGE 
48 hr SECTOR 

RANGE 
SFC 24 HR SECTOR 

RANGE 
48 HR SECTOR 

RANGE 
TYPE OF PCPN 
NATURE OF FGPC 
CLOUD TYPE 
SFC WIND DIR 
850 wind dir 



stationary front 

Colorado low 

mP mT 

WN¥ WSW 

200-400 km 

WSW 

<800 km 

ESE 

none 

SSW 

<200 km 

convective 

rain 

stratus 

SSW SSE 

WSW 



none 

none 

cA cm A mA 

WNW WSW 

200-400 km 

WSW WNW 

<800 km 

SSE ENE 

none 

WSW 

<200 km 

none 

snow 

stratus 

none 

none 



Table 4.3c: Characterization of high 
meteorological variables. 



SO, 



and NO. readings by the discrete 



In examining Table 4.3c, the trajectories associated with high SO^ and NO- 
measurements were similar but for the other meteorological fields the important 
factors were dissimilar. Synoptic conditions had some effect on the high SO^ 
readings but not on the high N0„ readings. The airmasses associated with them 
were also different, mP and oT air with SO, and cA, coA, and mA with N0». High 
SO, readings were associated with convective situations and rainfall while high 
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io, 







1982 


WARM 


MONTHS 


COLD 


MONTHS 




f 




P>F 


i 


P>F 


F 


F>F 


Synoptic type 


1 


35 


0.24 


1.16 


0.33 


1.88 


0.09 


Track of low 


i. 


57 


0.18 


1.19 


0.32 


1.18 


0.32 


Airmaa3 


2 


.90 


0.003 


1.81 


0.09 


1.42 


0.20 


850 mb 24 hr sector 


1. 


78 


0.11 


1.92 


0.10 


2.09 


0.09 


range 


2, 


17 


0.03 


0.44 


0.84 


0.76 


0.64 


48 hr sector 


1 


.07 


0.38 


1.17 


0.33 


0.58 


0.68 


range 


1. 


68 


0.09 


1.09 


0.38 


0.59 


0.80 


ifc 24 hr sector 


1 


.75 


0.10 


2.29 


0.04 


1.47 


0.19 


range 


1 


65 


0.13 


0.51 


0.77 


1.34 


0.25 


4# hr sector 


1 


.51 


0.17 


1,43 


0.21 


1,01 


0.44 


range 


2 


,47 


0.01 


1,90 


0.07 


0.76 


0.66 


type of pcpn 


5 


18 


0.006 


5.92 


0.02 


0.27 


0.76 


nature of pcpn 


8 


.29 


0.0001 


1.64 


0.21 


5.75 


0.001 


cloud type 


1 


.82 


0.09 


0.86 


0.53 


5.36 


0.0003 


sfc wind dir 


1 


.51 


0.16 


1.47 


0.19 


0.78 


0.61 


850 mb wind dir 


1 


.06 


0.39 


1.39 


0.23 


0.21 


0,97 



NO, 



1982 

F P>F 



WARM MONTHS 
F P>F 



COLD MONTHS 
F P>F 



Synoptic type 



0.56 0.76 



2.02 0.07 



Track of low 


0.43 


0.79 


0.90 


0.47 


Airmass 


1.50 


0-15 


0.81 


0.59 


850 mb 24 hr sector 


2.28 


0.04 


1.06 


0.38 


range 


1.10 


0.36 


0.82 


0.55 


48 hr sector 


1.99 


0.07 


1.30 


0.27 


range 


0.69 


0.73 


0.79 


0.61 


sfc 24 hr sector 


0.89 


0.52 


1.38 


0-22 


range 


1.52 


0.17 


0.27 


0,92 


48 hr sector 


0.82 


0.57 


0.73 


0.65 


range 


1.78 
0.02 
1.93 


0,07 
0.97 
0,12 


2.29 
4.49 
0.35 


0.03 


type of pcpn 


0.04 


nature of pcpn 


0.55 


cloud type 


3.69 


0.001 


1.12 


0.36 


sfc wind dir 


1,25 


0.27 


1.27 


0.28 


850 mb wind dir 






2.41 


0.03 



0.75 
0.83 
2 . 23 
5.46 
1.82 
1.94 
0.73 
1.24 
1.65 
1.06 
0.77 
0.2 5 
1.76 
2.88 
0.84 
0.41 






61 


0. 


51 





04 





001 





09 





11 





.68 





.29 





,15 





.40 





.66 





.78 





.16 





.01 





.56 





.87 



Table 4.4: An analysis of variance between SO, and NO^ and the discrete 
variables. Given are the F statistic and the significance probability 
(P>F). Underlined values are those deemed significant enough to 
warrant further statistical testing. 



MO. readlngi were associated 
precipitation. 



with snow but had no relationship to the type of 



A one way analysis of variance was done between the SO, and NO^ measurements and 
the discrete meteorological variables Indivlduariy to assess the slgnlflcanc© of 
the characterizations developed. Table 4.4 gives the F value and the significance 
probability for tests on the whole data set and tests dividing the data into cold 
and warm periods (procedure ANOVA in SAS User's Guide: Statistics 1982 Edition). 

From Table 4.4 the following list of variables (Table 4,5) having a significance 
level below 0.10| by season, and will be tested in further statistical analysis. 



SO. 



WARM MONTHS 



COLD MONTHS 



WARM MONTHS 



COLD MONTHS 



alrmass 

850 24 hr sector 
sfc 24 hr sector 
sfc 48 hr range 
type of pcpm 



synoptic type 
850 24 hr sector 
nature of pcpn 
cloud type 



synoptic type 
sfc 43 hr range 
type of pcpn 
850 wind dir 



alrmass 

850 24 hr sector 
850 24 hr range 
cloud type 



Table 4,5, Variables having a significance probability less than 0,10. 



4.2.5, Relationship of SO, 
Meteorological Factors 



and 



NO. 



Measurements to the Continuous 



Table 4.6 gives the correlation coefficients and significance probabilities 
between SO, and NO^ and the continuous meteorological variables (procedure CORR in 
SAS User's GuldeiBaslcs 1982 Edition). From this table, the variables with which 
the SO, had a correlation coefficient above 0.20 (either plus or minus) were 
surface temperature, precipi table water, duration of precipitation, and 
precipitation amount. For NO^ the significant correlations were only with the 
rainfall variables, namely maximum rainrate in mm/hr, duration of precipitation, 
precipitation amount, maximum rainfall in 6 hr, and the proportion of stations 
receiving precipitation. 



In a manner similar to the method of the previous section on discrete variables, 
the continuous meteorological variables behaviour for the high SO, and N0« 
readings were examined. Table 4.7 gives the averages and standard deviation for 
the continuous meteorological variables for the top 307, readings for SO, and N0«, 



Table 4,6i Correlatioii coefficients, slgnlflcaiice probafeillties, and number of observations between SO,, HO-, and t!it 
continuous meteorological variables, 

€0RRILATI01I COEFFICIENTS/PROB > R UTOER HOiRHO-0/NlJMBlE OF OBSllVATIONS 



id, 



m. 



so. 



Wo, 



tt 



CLOUD BASE HEIGHT 



SURFACE TEMPERATURE 



SURFACE DEWPOINT 



SURFACE WIND SPEED 



850 MB TEMPERATURE 



850 MB WIND SPEED 



CLOUD TOP HEIGHT 



CLOUD DEPTH 



PRECIPITABLE WATER 



-0.033 


-0.115 


0.69 


0.16 


150 


150 


0.320 


-0.128 


0.0001 


0.11 


161 


161 


0.297 


-0.151 


0.0001 


0.06 


161 


161 


-0.090 


-0.083 


0.25 


0.30 


162 


162 


0.337 


-0.083 


0.0001 


0.30 


162 


162 


-0.081 


-0.077 


0.31 


0.34 


153 


153 


-0.131 


-0.047 


0.13 


0,59 


134 


134 


-0.061 


-0.019 


0.52 


0.84 


116 


116 


0.289 


-0.121 


0.0002 


0.12 


160 


160 



TOTAL TOTALS INDEX 



RAINIATI (MM/15 MIN) 



RAINRATE (MM/30 MIN) 



RAINRATE (MM/HR) 



DURATION OF PRECIP 



AVERAGE PRECIP RATE 



AMOUNT OF PRECIP 



MAX 6 HR AMOUNT 



AREA OF PRECIP 



0,174 


-0.110 


0.03 


0.16 


162 


162 


-0.021 


-0.149 


0.83 


0,13 


106 


lOi 


-0.046 


-0.170 


0.64 


. 08 


106 


106 


-0.101 


-0.215 


0.30 


0.03 


106 


106 


-0.295 


-0.324 


0.002 


0.0008 


104 


104 


-0,068 


-0.100 


0.50 


0.31 


104 


104 


-0.213 


-0.380 


0.005 


0.0001 


174 


174 


-0.098 


-0.333 


0.22 


0.0001 


159 


159 


-0.1 38 


-0.328 


0.07 


0.0001 


176 


176 





SO^ 




NO -J 




MEAN 


SD 


MIAN 


SD 


CLOUD BASE HEIGHT 


603 


421 


609 


482 


SURFACE TEMPERATURE 


10.4 


8.7 


2.6 


10.3 


SURFACE DEWPOINT 


7.7 


8,3 


0,1 


9,7 


SURFACE WIND SPEED 


8.7 


4.8 


8.4 


5.0 


850 MB TEMPERATURE 


6.7 


7.5 


-0.1 


10.1 


850 MB WIND SPEED 


25 


11 


25 


11 


CLOUD TOP 


2973 


1893 


3417 


1870 


CLOUD THICKNESS 


2493 


1892 


2932 


1794 


PRECIPITABLE WATER 


22.9 


11.8 


16.3 


11.4 


TOTAL TOTALS 


15.0 


18.5 


9.5 


18.0 


RAINRATE min/15 min 


2.1 


1.7 


1.5 


0,9 


mm/ 30 mln 


2.7 


2.3 


1.9 


1.1 


mm/hr 


3.3 


2.8 


2.2 


1.3 


DURATION OF PCPN 


5.7 


3.3 


5,5 


3,7 


AVERAGE RAINRATE ram/hr 


1.0 


0.9 


0.9 


0,8 


AMOUNT OF PCPN 


4 . 9 


3.7 


3.5 


1.9 


MAX 6 HR AMOUNT 


3.5 


3.2 


2.3 


1,9 


7. STNS WITH PCPN 


0.61 


0.27 


0.58 


0.23 


Tab 1 e 4.7: Average s 


and s 


tandard 


devla 


it ions 



(SD) of the eontiauous 

meteorological variables for the top 307, readings of the SO, and N0«. 

4 3 



Wheii the data are split by season (warm vs cold months) the significant parameteri 
displayed some interesting variation. For example, airnass was significant for 
the SO, in the warm months while for the N0«, it was significant in the cold 
months. This is in keeping with the previous statement that high SO, was 
associated with mP and mT air and high NO- with cA, cmA and mA air. The nature of 
the precipitation was significant for SD. In the winter when there were many 
examples of rain and snow, as opposed to the summer when only rain was seen so no 
useful discrimination would be possible. For the type of precipitation the 
reverse was true. The distinction between continuous and convective precipitation 
would be much less clear and more difficult to observe in winter. For NO* the 
synoptic type and the type of precipitation, which previously were found not to be 
related to the N0„ readings, were significant In the warm months. During this 
time high NO- measurements were associated with warm fronts and stationary fronts, 
and with convective precipitation. During the cold months, high N0» readings are 
associated with other factors such that when viewed over the year, the synoptic 
type and precipitation type were not highly related to the NO- readings. The 
synoptic type was significantly related to SO, in winter. In this case there was a 
higher frequency of synoptic types for which the SO, readings were low (cold 
fronts, cold lows, occluded lows) 



a. 



k comparison of Table 4.7 with Table CI In Appendix C (the means and standard 
deviation for the whole data set) gives the variables with somewhat higher or 
lower averages In the top 307. of the 
Table 4.8. 



SO. and NO. readings and Is summarized In 



SO. 



HIGHER MEANS 



LOWER MEANS 



HIGHER MEANS 



LOWER MEANS 



surface temperature 
surface dewpolnt 
850 mb temperature 
preclpl table water 
total totals Index 



cloud base height 

surface wind speed 

cloud top height 

cloud thickness 

duration of 

precipitation 

precipitation 

amount 

max 6 hr precipitation 

amount 

proportion of stations 

with precipitation 



none 



cloud base height 
surface temperature 
surface dewpolnt 
surface wind speed 
350 mb temperature 

preclpitable water 

total totals index 

all precipitation 
variables except 
average pcpn rate 



Table 4.B. Variables which had means somewhat different from means in the whole 



data set when 
measurements. 



stratified into the top 30X of the SO, 



NO, 



The general trend from Table 4.8 is that high SO, readings were associated with 
high temperatures and moisture and low precipitation events, low cloud bases, and 
low windspeeds. For N0„ high readings were associated with low temperatures and 
moisture, low wind speeds, low cloud bases, and low precipitation events. 

To further examine this, least squares estimates were used to construct linear 
regression models of SO, and N0« values with the continuous meteorological 
variables individually tprocedure REG In SAS User's Guide: statistics 1982 
Idition). The nature of the relationships and their significance can be better 
assessed. Table 4.9 gives the slope parameter bl (y=b +b *x) and Its significance 
probability for the whole year and the warm and cold months individually. 

In Table 4.9, those variables which have the highest correlations with the SO, and 
N0„ measurements are underlined and itemized In Table 4.10. 



tt 



BQ, 





1982 




WARM MONTHS 


COLD MONTHS 






m 


P>F 


bl 


P>F 


m. 


P>F 


CLOUD BASE HEIGHT 


0.001 


0.26 


0.002 


0.24 


-0.001 


0.15 


SURFACE TEMPERATURE 


0.13 


0.02 


0.08 


0.39 


-0.01 


0.89 


SURFACE DEWPOINT 


0.14 


0.02 


0.08 


0.44 


0.01 


0.91 


SURFACE WIND SPEED 


-0.11 


0.14 


-0.12 


0.35 


0.01 


0.89 


850 MB TEMPERATURE 


0.14 


0.02 


0.13 


0.15 


-0.0002 


0.99 


850 MB WIND SPEED 


-0.04 


0.15 


-0.01 


0.86 


0.01 


0.65 


CLOUD TOP 


-0.0001 


0.82 


-0.0001 


0.79 


-0,00003 


0.89 


CLOUD THICKNESS 


-0.0001 


0.67 


-0.0001 


0.65 


0.00002 


0.90 


PRECIf ITABLE WATER 


0.12 


0.001 


0.12 


0.025 


0.0004 


0.99 


TOTAL TOTALS 


0.08 


0.01 


0,11 


0.05 


-0.001 


0.98 


RAINRATE mm/ 15 mln 


0.28 


0.25 


0.14 


0.65 


-0.26 


0.51 


mm/ 30 mln 


0.17 


0.36 


0.09 


0.68 


-0,26 


0.33 


mm/hr 


0.07 


0.65 


0.05 


0-80 


-0.21 


0.25 


DURATION OF PCPN 


-0.22 


0.004 


-0.22 


0.08 


-0.05 


0.51 


AVERAGE RAINRATE mm/hr 


0.27 


. 54 


0.04 


0.94 


-1.24 


0.13 


AMOUNT OF PCPN 


-0.14 


0.03 


-0.12 


0.25 


-0.11 


0.05 


MAX 6 HR AMOUNT 


-0.14 


0,17 


-0.17 


0.21 


-0.12 


0.23 


X STNS WITH PCPN 


-3.04 


0,12 


-2.95 


0.34 


-0.95 


0,54 



NO. 





1182 




warm i 


lonths 


esld months 




W 


P>P 


bl 


P>F 


m 


P>F 


CLOUD BASE HEIGHT 


. 00001 


0.95 


0,0001 


0.58 


-0.003 


0.09 


SURFACE TEMPERATURE 


0.004 


0.57 


-0.01 


0.58 


-0.01 


0.32 


SURFACE DEWPOINT 


. 004 


0.58 


-0.01 


0.46 


-0.01 


0,45 


SURFACE WIND SPEED 


-0 . 02 


0.02 


-0.02 


0.13 


-0,01 


0.40 


850 MB TEMPERATURE 


0.001 


0.87 


-0.004 


0.67 


-0.01 


0,23 


850 MB WIND SPEED 


-0.01 


0,06 


-0.003 


0.68 


-0.002 


0,59 


CLOUD TOP 


-0.00001 


0,71 


-0,00001 


0,80 


-0.00001 


0.70 


CLOUD THICKNESS 


-0.00001 


0.70 


-0.00001 


0,74 


-0.000003 


0,93 


PREC IP ITABLE WATll 


0.004 


0.43 


0.001 


0.90 


-0.01 


0.28 


TOTAL TOTALS 


. 004 


0.33 


0,01 


0.32 


-0,01 


0,28 


RAINRATE mm /1 5 mln 


-0.001 


0,95 


-0.01 


0.73 


-0.01 


0,11 


mm/ 30 rain 


-0.01 


0.77 


-0.01 


0.73 


-0.09 


0.06 


mm/hr 


-0.01 


0.49 


-0.01 


0,66 


-0.07 


0.05 


DURATION OF PCPN 


-0.03 


0.002 


-0.03 


0,02 


-0,01 


0.60 


AVERAGE RAINRATE mn/h 


r-0.02 


0.76 


-0.03 


0.61 


-0,03 


0,04 


AMOUNT OF PCPN 


-0.02 
-0.03 


0,01 
0.03 


-0.02 
-0.03 


0,10 
0.09 


-0,02 
-0.04 


0.06 


MAX 6 HR AMOUNT 


0.05 


I STNS WITH PCPN 


-0,39 


0,11 


-0,30 


0.39 


-0.24 


0.43 



Table 4.9: Slope estimate and significance probability between SO,, N0«| and the 
continuous meteorological variables. 



S 



so, 



1982 



WARM MOHTHS 



COLD MOITTHS 



surface temperatore 

surface dewpolnt 

850 mb temperature 

preclpitable water 

total totals 

duration of precipitation 

precipitation amount 



preclpitable water 

total totals 

diiratlon of precipitation 



cloud base height 
precipitation amouat 



HO. 



1982 



WARM MONTHS 



COLD MONTHS 



surface wind speed surface wind speed 

duration of precipitation duration of precipitation 

precipitation amount precipitation amount 

max 6 hr amount max 6 hr amount 



precipitation amount 
max 6 hr amount 
cloud base height 



Table 4.10. Variables slgntficantly correlated with SO. and 10-, 



For SO, , splitting the data into warm and cold months reduces the number of 
significant variables. In essence the splitting removes the Information which was 
contained In many of the correlations, for example surface and 850 mb temperature, 
surface dewpolnt. For NO-, except for surface wind speed, which Itself was very 
marginal, only precipitation variables show significant correlations. 



,24 



5. Statistical Model 

5.1 SQUthwest Region 

Using the ideas and results of the previous section, a statistical modal can be 
developed to predict SO. and N0„ values. A separate model can be developed for the 
warm and cold seasons. The teclnlqiie to be used is the analysis of covariance 
(ANCO¥A) where 2 discrete and 2 continuous meteorological variables will be used 
(procedure GLM in SAS User's Guides Statistics 1982, Edition pg 140). The choice of 
2 discrete variables was arbitrary but done in order that a sufficient number of 
cases would be present for each combination of the discrete variables. For the 
continuous variables, the correlations among the meteorological variables and 
between the SO. and NO- and the meteorological variables indicated that 2 continuous 
variables would be sufficient. In theory every conceivable combination of 
meteorological variables could be tested. But the findings of the previous section 
have provided the basis for eliminating many of the meteorological variables from 
consideration and resulted in a more manageable set of variables for statistical 
modelling tests. 

Based on Tables 4.5 and 4.10 the following ANCOVA analyses were carried out. Table 
5.1 gives the 4 variables plus the interaction term between the 2 discrete 
variables, their individual F value and significance probability (P>F), the F value 
for the set of variables, Its significance probability, the proportion of the 
variance of the dependent variable (SO, or HO,) explained by the model (R), and the 
mean squared .error of the difference between the model's predictions and 
observations (E ), The sum of the errors would by design be 0. The optimal situation 
would be low significance probabilities, high R, and low E . 



Table 5,1: Results of the ANCOVA. 

WARM MONTHS 

SO, 
4 

VARIABLl i P>F F P>F R E^ 

type of pcpn 3.78 0,06 1,02 0.46 0.27 7.62 

airmass 

Interaction 

preclpltable water 

duration of pcpn 

sfc 48 hr range 1.18 0.35 0.57 0.92 0.35 6.79 

airmass 

interaction 

preclpltable water 

duration of pcpn 



3.78 


0,06 


1.13 


0.37 


0.38 


0.82 


0.09 


0.77 


0.98 


0.33 


1.18 


0.35 


0.66 


0.70 


0.18 


0.99 


0.14 


0.71 


. 09 


0.77 



c ont ' d 



ii 



Table 5.1 cont'd 

850 24 hr sector 2.21 0.07 1.44 0.19 0.29 7,62 

type of pcpn 

interaction 

precipitable nater 

duration of pcpn 

sfc 48 hr range 1.05 0.42 0.85 0.62 0.27 8.12 

type of pcpn 

interaction 

precipitable wate^ 

duration of pcpn 

type of pcpn 3.78 0.06 1.02 0.46 0.28 7.62 

airmass 

interaction 

precipitable water 

duration of pcpn 

type of pcpn 4.97 0.03 2.69 O.Oi 0.47 5.96 

sfc 24 hr sector 

interaction 

precipitable water 

duration of pcpn 



2.21 


0,07 


2.10 


0.16 


0.69 


0,57 


0.05 


0.82 


0.60 


0,44 


1.05 


0.42 


2,31 


0.14 


0,48 


0.75 


0.24 


0.62 


0.02 


0.88 


3.78 


0.06 


1.13 


0.37 


0.38 


0.82 


0.09 


0,77 


0,98 


0.33 


4.97 


0.03 


3.13 


0.01 


2.32 


0.09 


0.68 


0.41 


0.82 


0.37 


NO3 




2.51 


0.12 


0.87 


0.54 


0.35 


0.84 


1.43 


0.24 


3,48 


0.07 


2.11 


0.08 


1.16 


0.36 


1.12 


0.38 


3.10 


0.09 


1.65 


0.21 


2.53 


0.12 


0.88 


0,53 


0,36 


0,84 


1.44 


0.24 


3.80 


0.06 



type of pcpn 2.51 0.12 1,07 0,42 0.29 0.21 

sfc 48 hr range 

interaction 

duration of pcpn 

sfc wind a peed 

synoptic type 2.11 0.08 1.43 0.17 0.56 0,14 

sfc 48 hr range 

interaction 

duration of pcpn 

amount of pcpn 

type of pcpn . 2.53 0.12 1.10 0.39 0,30 0.21 

sfc 48 hr range 

interaction 0,36 0,84 0.30 

duration of pcpn 

amount of pcpn 

cont'd 



Table 5.1 cont'd 



COLD 


MONTHS 
















SO, 


















8.02 


0.0002 


4.81 


0,0001 


0.59 


2.19 


2.59 


0.03 
















3.65 


0.001 
















5.11 


0.03 
















18.63 


0.0001 
















7.73 


0,0001 3 


.91 





.0001 


0. 


,62 


1. 


,99 


1.63 


0,15 
















2.54 


0.01 
















4.28 


0.04 
















10,72 


0.002 
















9.15 


0,0001 7 


.22 


_0 


.0001 


_0_ 


.64 


1, 


.93 


6,87 


0.0005 
















5.50 


0.0001 
















5,18 


0.03 
















9.09 


0.004 
















m,^ 


















3,65 


0.004 2. 


17 


0. 


01 0. 


56 


0. 


20 




1.06 


0.41 
















1.29 


0.24 
















0.48 


0.49 
















17.90 


0,0001 

















naturs of pcpm 
synoptic type 
interaction 
cloud base height 
amount of pcpn 

cloud type 
synoptic type 
interaction 
cloud base height 
amount of pcpn 

cloud type 

850 24 hr sectoi: 

interaction 

cloud base height 

asiount of pcpn 



cloud type 

850 24 hr range 

interaction 

cloud base height 

amount of pcpn 

cloud type 3.75 0.003 2.63 0.001 0.52 0.22 

airnass 2.06 0.06 

interaction 1.36 0.22 

cloud base height 0.07 0,79 

amount of pcpn 15.13 0.0003 

850 24 hr range 1-37 0.23 1.56 0.07 0.51 0.22 

airmass 1.89 0.09 

interaction 1.08 0,40 

cloud base height 1.41 0.24 

amount of pcpn 9.41 0,003 

cloud type 5.37 0.0001 6.96 0.0001 0.64 0.17 

850 24 hr sector 12,32 0.0001 

interaction 6.05 0.0001 

cloud base height 1,09 0.30 

amount of pcpn 11.76 0,001 



From Tabla 5.1, the variables which conprised the best set of variables, by tinie of 
year, were the following (their values are under lined lii table 5.1)i 



WARM MONTHS 



for SO, 



type of precipitation 
surface 24 hr sector 
preclpltable water 
duration of precipitation 



for H0> 



synoptic type 
surface 48 hr range 
duration of precipitation 
amount of precipitation 



COLD MONTHS 



for both SOx AND NO^ 



cloud type 
850 24 hr sector 
cloud base height 
amount of preclpitatlom 

The combinations of discrete variables which perforoed the best Invariably had very 
strong interaction terms. Appendix F, which contains contingency tables for the 
pairs of discrete variables which were chosen (Tables Fl to F3), shows the 
relationship between the different categories of the 3 pairs of discrete variables. 
Referring back to Table 4.2, the type of precipitation was highly related to surface 
24 hr sector as was cloud type to 850 mb 24 hr sector. Synoptic type and surface 48 
hr range on the other hand were not so highly related. Also given in Appendix A are 
the scatter plots of SO, or NO, versus the continuous meteorological variables 
chosen (Figures Fl to F8). They show the nature of the relationship, as summarized 
In Table 4.9, and establish that there were not some extreme points which controlled 
the regression. 

Summarizing, the situations which will characterize high acidic measurements are the 
following. For SO, In the warm months, convectlve situations with the 24 hr surface 
trajectory from the ESE and preclpltable water high and duration of the 
precipitation event low were the aspects of the 4 chosen variables associated with 
high SO. readings. For the cold months, for both SO, and NO^, they are 850 mb 24 hr 
trajectories from the WNW with overcast stratus and low cloud bases and low 
precipitation amounts. For NO- in the warm season, warm or stationary fronts with 
surface 48 hr trajectories originating less than 200 km away giving low 
precipitation duration and amount were the aspects of the 4 chosen meteorological 
variables associated with high N , readings. 

The preceding judgments are somewhat open to conjecture. However it was reassuring 
to note that generally there was little departure from noriiallty for the errors in 
prediction. But the test of the model on an independent data set, e.g. 1983, will 
ultimately determine the model's worth in trying to be able to estimate annual 
deposition based on a limited set of Eulerlan model runs and a specified time period 
of observed deposition and meteorological data. 



5,2- lorthwest Regioti 

In order to see how conditions vary, an analysis of variance was done between the 
SO, and N0<. measurements and the meteorological variables for the northwest region, 
the one most dissimilar from the southwest region. Table 5,2 gives the results of 
the analysis of variance (la SAS User's Guide: Statistics 1982 Edition AN0V4 
procedure), again split into warm and cold months (April to September vs the rest). 



Table 5,2i Analysis of variance for the northwest region. 

SO, 

4 

WARM MONTHS COLD MONTHS 
f P>F F F>F 



Synoptic type 


0.75 
0.17 


0.61 
0.68 


2.44 
0.73 


0.04 


Track of low 


0,57 


Airmass 


2.17 


0,04 


4-03 


0.002 


850 mb 24 hr sector 


2,63 


0-02 


0.76 


0,63 


range 


0.70 


0,65 


1.97 


0.08 


48 hr sector 


3.86 


0,002 


1.32 


0.27 


range 


1,26 


0.28 


0.84 


0.58 


sfc 24 hr sector 


0.77 


0.62 


0.95 


0,48 


range 


1.97 


0.09 


2.81 


0.02 


48 hr sector 


3.03 


0.008 


0,88 


0.53 


range 


3.41 


0.004 


0.23 


0.98 


type of pcpn 


2.49 


0.12 


0.59 


0.45 


nature of pcpn 


2,09 


0,13 


. 96 


0.39 


cloud type 


2.07 


0-08 


0.26 


0,93 


sfc wind dir 


1.00 


0,44 


0-81 


0.60 


850 mb wind dir 


0.39 


0.88 


0.92 


0.50 



'CO'Dt d 
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Table 5.2 cont'd 



NO. 



WARM MONTHS GOLD MONTHS 



P>F 



P>F 



Synoptic type 


0.81 


0.56 


3.64 


0. 


004 


Track of low 


0.06 


0.81 


0.49 





.74 


Alrmass 


3.03 
1.82 


0.005 
O.li 


2.43 
2.23 


0. 


04 


850 mb 24 hr sector 


0, 


05 


range 


0.36 


0.90 


0.89 





,52 


48 hr sector 


1.92 
0.98 


0.09 
0.46 


4.46 
0.57 





001 


range 


0. 


80 


sfc 24 hr sector 


1.64 


0.14 


2.00 





.07 


range 


2.69 


0.03 


1.40 


0. 


23 


48 hr sector 


3.19 


0.006 


1.04 





.42 


range 


1.56 


0.16 


0.48 





.84 


type of pcpn 


0.69 


0.41 


0.06 


0. 


80 


nature of pcpn 


2.03 


0.14 


0.35 





.70 


cloud type 


1.12 


0.36 


0.46 


0. 


80 


sfc wind dir 


1.05 


0.40 


0.88 





.53 


850 mb wind dlr 


0,66 


0.68 


1.14 


0, 


36 



Also, linear correlations between SO, and N0„ readings and the continuous 

meteorological variables were done by time of year for the northwest region 

(procedure REG in SAS User's Guide; Statistics 1982 Edition). The slope parameter of 

the regression along with the significance probability of the relationship is given 
in Table 5.3. 






Table 5.3: Slope estimate and slgnlfieance probability betifeen SO. NO. 
continuous meteorological variables. 



and t'he 



SO. 



WARM MONTHS 



COLD MONTHS 



iii 



P>F 



bl 



F>F 



CLOUD BASE HEIGHT 
SURFACE TEMPER ATURl 
SURFACE DEWPOINT 
SURFACE WIND SPEED 
850 MB TEMPERATURE 
850 MB WIND SPEED 
CLOUD TOP 
CLOUD THICKNESS 
PRECIPITABLE WATER 
TOTAL TOTALS 
RAINRATE mm/ 15 min 

mm/ 30 mill 

mm/hr 
DUIATION OF PCPN 
AVERAGE RAINRATE mm/hr 
AMOUNT OF PCPN 
MAX 6 HR AMOUNT 
I STNS WITH PCPN 



0.0009 


0.02 


-0.0003 


0.46 


0.64 


0.01 


0.019 


0.18 


0.067 


0.01 


0.018 


0,17 


0.024 


0.54 


-0.021 


0,53 


0.034 


0.13 


0.02 


0.26 


0,02 


0.22 


0.013 


0.39 


0.0001 


0.40 


0.0001 


0.52 


0.00003 


0.77 


0.0002 


0.26 


0.032 


0.08 


0,036 


0.19 


0,011 


0.34 


0,008 


0-50 


0.017 


0.72 


-0.039 


0.83 


0.016 


0.67 


-0.03 


0.81 


0.012 


. 69 


-0.039 


0.65 


0.014 


0.52 


0-012 


0.59 


0.026 


0.81 


-0.024 


0,92 


•0.032 


0.14 


-0.026 


0.38 


■0,023 


0.31 


-0.041 


0,40 


■1.004 


0.03 


-0.452 


0,38 



NO, 



WARM MONTHS 



COLD MONTHS 



||:|: 



P>F 



bl 



P>F 



CLOUD BASE HEIGHT 
SURFACE TEMPERATURE 
SURFACE DEWPOINT 
SURFACE WIND SPEED 
850 MB TEMPERATURE 
850 MB WIND SPEED 
CLOUD TOP 
CLOUD THICKNESS 
PRECIPITABLE WATER 
TOTAL TOTALS 
RAINRATE mm/15 min 

ram/ 30 min 

mm/hr 
DURATION OF PCPN 
AVERAGE RAINRATE mm/hr 
AMOUNT OF PCPN 
MAX 6 HR AMOUNT 
7. STNS WITH PCPN 



0.0001 


0.13 


-0.000002 


0.98 


0.009 


0.01 


-0,003 


0.32 


0.009 


0.01 


-0,003 


0.32 


0.0024 


0.67 


-0.009 


0.24 


0.003 


0,40 


-0.004 


0,37 


0.001 


. 54 


0.002 


0.62 


0.000002 


0.90 


-0.000003 


0.92 


0.00002 


0.90 


0.00001 


0,74 


0.004 


0.16 


-0,007 


0.24 


0.001 


0.47 


-0.002 


0,43 


0.001 


0,91 


-0.012 


0.71 


0.00013 


0.98 


-0,013 


0.57 


0.007 


0.86 


-0,012 


0.43 


0-002 


0.52 


-0.006 


0.15 


0.008 


0.54 


-0.043 


0.40 


0.003 


0.28 
0.40 


-0.009 
-0.014 


0.10 


0.002 


0.13 


0.158 


0.02 


-0.133 


0.21 



Using Tables 5.2 and 5.3, ttie variables whicli were deemed to show enough promise in 
their relation to SO. and NO, (and whose significance probablilties are underlined) 
to warrant further testing In ANC07A are given in Table 5.4* The significance 
probabilities for the continuous variables in the cold season were quite high and 
made selection of any variables then somewhat dubious. 



SO. 



NO, 



WARM MONTHS 



COLD MONTHS 



WARM MONTHS 



COLD MONTHS 



airmass 

850 48 hr sector 
sfc 48 hr range 
cloud type 

cloud base height 
surface temperature 



airmass 
synoptic type 
sfc 24 hr range 



airmass 

sfc 48 hr sector 



airmass 

850 48 hr sector 



surface temp. surface temp. anount of pcpn 
precipi table water stations with pcpn 



Table 5.4: Meteorological variables significantly related to SO. and NO^, 

For the variables listed in Table 5.4 the ANCOVA was run. The results are given In 
Table 5,5. 

Table 5.5: Results of ANCOVA for northwest region. 

WARM MONTHS 

to. 



VARIABLE 



W 



f>F 



P>F 



sfc 48 hr range 
airmass 
interaction 
cloud base height 
sfc temperature 

sfc 48 hr range 
cloud type 
interaction 
cloud base* height 
sfc temperature 

airmass 

850 48 hr sector 

interaction 

cloud base height 

sfc temperature 



10.18 


0.0001 


6.89 


0.0001 


1.93 


0.072 


0.09 


0,772 


9.75 


0.004 


3.98 


0.004 


0,65 


0.662 


0.50 


0.892 


. 34 


0.563 


1.36 


0.25 


2.73 


0.023 


3.78 


0.005 


0.54 


0.864 


0.07 


0.792 


0,61 


0,441 



5.27 0.0001 0.83 0.39 



1.43 0.17 0.51 1.09 



1.86 0.04 0.58 1.00 



Gont' d 



Taile 5.5 cont'd 



cloud type 


3.04 


0.020 


2.46 


0.007 0.56 


850 48 hr sector 


3.19 


0.012 






interaction 


1.00 


0.453 






Cloud base height 


0.11 


0.742 






sfc temperatiire 


9.26 


0.004 






alrmas s 


3.43 


0.005 


2.78 


0,002 0.62 


cloud type 


3.82 


0.0061 






interaction 


2.21 


0.032 






cloud base height 


0.16 


0,693 






sfc temperature 


1.93 


0,172 







1.03 



0.84 



NO. 



atraiass 7.60 0.001 

sfc 48 hr sector 3.45 0.007 

interaction 2.45 0.012 

stns with pcpn 0.00 0.992 

sfc temperature 1.83 0.186 



3.67 0.0002 0.81 0.01 



COLD MONTHS 



VARIABLE 



SO, 

4 



iW' 



alrmas s 


2.55 


sfc 24 hr range 


1.74 


interaction 


0.85 


sfc temperature 


1.04 


precipltable water 


0.79 


sfc 24 hr range 


2.67 


synoptic type 


1.30 


interaction 


0.91 


sfc temperature 


0.03 


precipltable water 


3.69 



0.048 
0.158 
0.585 
0.318 
0.382 

0.041 
0.297 
0.539 
0.866 
0.067 



1,50 0.17 0.61 0.66 



1.53 0.16 0.61 0.65 



NO. 



airmass 


3.49 


0.010 


850 48 hr sector 


2.05 


0.091 


interaction 


2 . 36 


0.044 


pcpn amount 


1.13 


0.296 



2.53 0.01 0.65 0.01 
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Frotn Table 5.5 the set of meteorological variables whleb provide the best 

statistical model (and are imderllcied) are the following: 



SO, 



NO, 



WARM MONTHS 



COLD MONTHS 



WARM MONTHS 



COLD MONTHS 



airmass 

sfc 48 hr range 

cloud base height 
sfc temperature 



sfc 24 hr range 
synoptic type 



airmass 

sfc 48 hr sector 



airmass 

850 48 hr sector 



sfc temperature sfc temperature amount of pcpn 
precipltable water stations with pcpn 



A comparison of this list of variables with the list composed in the southwest 
region analysis reveals that for the discrete variables, airmass was a more 
significant variable in the northwest region while precipitation type and cloud type 
were more significant in the southwest region. Both areas relied on a trajectory 
variable as a significant variable. For the continuous variables, precipitation 
variables were more predominant in the southwest (which is not surprising considered 
the differences in the rain gauge networks) while temperatures were more prominant 
in the northwest. 

It is interesting to compare these results with those that would have been obtained 
if the significant variables from the southwest region were used in the ANCOVA 
analysis for the northwest region. In this case the results of the ANCOVA are shown 
in Table 5.6. 



Table 5.6j ANCOVA results using the 
southwest region. 



same variables as for the analysis in the 



VARIABLE 



P>F 



P>F 



WARM MONTHS 



SO, 



type of pcpn 
sfc 24 hr sector 
interaction 
precipltable water 
duration of pcpn 



4.97 0.03 

0.89 0.52 

3.17 0.04 

0.00 0.94 

0.11 0.74 



1.60 0,13 0.37 0.92 



NO, 



synoptic type 
sfc 48 hr range 
interaction 
duration of pcpn 
amount of pcpn 



2.16 


0.75 


5.08 


0.0007 


1.44 


0.198 


0.01 


0.91 


0.36 


0.55 



2.42 O.Ol 0.69 0.01 



cont'd 



Table 5.6 cont'd 
COLD MONTHS 



cloud type 

850 24 hr sector 

interaction 

cloud base height 

amount of pcpn 



SO. 



0,25 


0.94 


0.60 


0.75 


1.07 


0.40 


0.11 


0.74 


0.03 


0.87 



0.60 0.88 0.29 1.12 



NO, 



cloud type 

850 24 hr sector 

interaction 

cloud base height 

amount of pcpn 



0.48 0.79 

1.53 0.20 

1.29 0.29 

0.20 0.66 

0.28 0.60 



1.07 0.43 0.42 0.03 



A comparison of Tables 5.1, 5.4, and 5.6 shows that the variables picked out in the 

analysis of the southwest region did not perform nearly so well in the northwest 

when compared to the ANCOVA results of these variables in the southwest or when 

compared to the ANCOVA results for the northwest when the variables used were 

determined by the analysis of the northwest data (Table 5.7). In general the ANCOVA 

analyses for the 2 regions were comparable for the cold months (when the variables 

used in the northwest analysis were chosen from the analysis of the northwest data) 

but for the warm months the ANCOVA analysis was better in the northwest than the 

southwest regions in that it had smaller errors and explained more of the SO. and 

NO^ variance. This effect would largely be attributable to the smaller range and 

means of the SO, and N0„ in the northwest in the warm months, 
4 3 



R (from ANCOVA) 
SW NWCSW) NW 



WARM MONTHS 



COLD MONTHS 



Table 5,7: 



SO, 



SO, 



0.47 0.37 0.83 



NO- 0.56 0.69 0.81 



0.64 0,29 0.61 



10^ 0.64 0.42 0,65 

A comparison among the ANCOVA analyses for the southwest region (SW), 
the northwest region when the variables used in the analysis were the 
same as those used in the southwest analysis (NW(SW)), and the 
northwest region when the variables used were determined by an analysis 
of the northwest region (NW). Displayed is the amount of variance in 



the dependent 
(R). 



variable (SO, 



or N0„) 
3 



explained by the ANCOVA analysis 
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Hoiiever, the net resylt Is that the variables found to be significant In one 
region are not necessarily the best in another region. This should be neither 
surprising nor disturbing. The quality of the meteorological data in the 2 
regions Is not the same. The southwest region has a much more dense rain gauge 
network, the APIOS stations are more closely grouped , and the main hourly site 
(London) is well located amongst the 4 stations. This is not the case in the 
northwest sector. Secondly, the meteorological conditions are significantly 
different. As was previously stated, fronts play a larger role in the southwest 
region and occluded systems in the northwest. Consequently the type of alrmass 
was a more significant factor In the northwest while cloud type and precipitation 
type were bigger factors In the southwest. The result then is that each region 
should be considered on its own, with its meteorological peculiarities and data 
base taken into account. Thus each region would have its own tailored statistical 
model. However this is not to say that the processes involved in acid rain are 
different in each region. As was shown for the southwest region, the different 
meteorological variables were highly related. The fact that one may be more 
appropriate In a particular region is due to the uniqueness of that region to the 
meteorological conditions and the non-uniform degree of the data between regions, 
A second year of data would help further clarify this point. 



5,3 Central Region 

An analysis of variance was performed between the discrete meteorological 
variables and the SO, and N0„ measurements, in a manner similar to the northwest 
region. The results are given in Table 5.8, where the underlined variables are 
the ones which were retained for the covariance analysis. In this regard, the 
variables which were kept were those which had low significance probabilities and 
most consistent with the previous analyses. Thus, the selection process was more 
stringent here than previously. 

Table 5.9 contains the results of the linear regression between the continuous 
meteorological variables and the SO, and NO^ readings. Here also, the desire to 
have consistency with the variables chosen In the first two regions studied 
strongly affected the choice of variables to be used in the ANCOVA analysis. In 
addition, cloud top and cloud thickness were thought to have too many missing 
values to be useful. 



WAEM MONTHS COLD MONTHS 
P>F F P>F 



Synoptic type 


1.01 


0.43 


1.44 


0.21 


Track of low 


1.48 


0.22 


2.75 


0.02 


Alrmass 


1.84 


0.09 


2.56 


0.02 


850 mb 24 hr sector 


1.93 


0.09 


0.52 


0.67 


range 


2.35 


0,04 


0.62 


0.78 


48 hr sector 


3.87 


0,004 


1.86 


0.13 


range 


2,22 


0.03 


1.23 


0.28 


sfc 24 hr sector 


3,21 


0.005 


2.06 


0.06 


range 


1.86 


0,11 


0.62 


0.71 


48 hr sector 


3.51 


0.003 


7.56 


0.0001 


range 


0.96 


0.47 


0.50 


0.85 


type of pcpn 


3.71 


0.06 


0.15 


0.86 


nature of pcpn 


1.01 


0.37 


1.39 


0.25 


cloud type 


0.53 


0.75 


0.45 


0.87 


sfc wind dir 


2.69 


0.01 


1.44 


0.19 


S50 mb wind dlr 


1-53 


0.18 


0.66 


0.68 



N03 

WARM MONTHS COLD MONTHS 
f P>F F P>F 



Synoptic type 


1.33 


0.26 


1.74 


0.11 


Track of low 


0.35 


0.84 


0.83 


0.53 


Alrmass 


1.67 


0.13 


0.70 


0.67 


850 mb 24 hr sector 


1.91 


0.09 


1.26 


0.29 


range 


2.22 


0.05 


1.19 


0,31 


48 hr sector 


3.05 


0.01 


2.22 


0.07 


range 


2.31 


0.02 


1.00 


0.45 


sfc 24 hr sector 


3.94 


0.001 


3.05 


0,007 


range 


1.61 


0.17 


0.94 


0.47 


48 hr sector 


2.92 


0.01 


4.26 


0.001 


range 


0.62 


0.74 


0.75 


0.65 


type of pcpn 


2.33 


0.13 


0,43 


0.65 


nature of pcpn 


1.45 


0.24 


2.32 


0.06 


cloud type 


0.50 


0.78 


0.92 


0.50 


sfc wind dir 


1.53 


0.16 


1.11 


0,36 


850 mb wind dir 


2.40 


0.04 


1.00 


0.43 



Table 5.8: Analysis of Variance of the discrete variables for Central legion 



if 



so. 



CLOUD BASl HEIGHT 
SUEFACE TEMPERATURE 
SimFACE DEWPOINT 
SURFACE WIND SPEED 
850 MB TEMPERATURE 
850 Ml WIND SPEED 
CLOUD TOP 
CLOUD THICKNESS 
PRECiPITABLE WATII 
TOTAL TOTALS 
RAINIATI mm/15 mln 

mm./ 30 rain, 

mm/hr 
■DUIATION OF PCPN 
AVEEAGE RAI'NRATE nm 
AMOUNT OF PCPN 
MAX 6 HR AMOUNT 
Z STOS WITH PCPN 



/hr 



CLOUD BASE HEIGHT 

SURFACE TEMPERATURE 
SURFACE DEWPOINT 
SURFACE WIND SPEED 
850 MB TEMPERATUIE 
850 MB WIND SPEED 
CLOUD TOP 
CLOUD THICKNESS 
PEICIPITABLE WATER 
TOTAL TOTALS 
RAINIATE mm/is min 

mm/ 30 min, 

mm/hr 
DUIATION OF PCPN 
AVEMGl lAIMATE mm/ht 
AMOUNT OF PCPN 
MAX 6 HI AMOUNT 
X STNS WITH PCPN 



WARM MONTHS 


COLD MONTHS 


bl 


P>F 


bl 


P>F 


0.0001 


0.85 


0,0006 


0.24 


0.05 


0.29 


0.09 


0.002 


0.002 


0.97 


0.08 


0,003 


-0.011 


0,87 


0.015 


0.76 


0.004 


0,93 


0.130 


0.001 


0.024 


0.42 


-0.016 


0.59 


-0.0002 


0.40 


-0.0004 


0.01 


-0.0001 


0.52 


-0.0003 


0.05 


0.02 


0.53 


0.100 


0.02 


-0,005 


0.82 


0.019 


0.42 


0.201 


0,38 


1,175 


0.26 


0.098 


0.54 


0,795 


0.28 


0.034 


0.79 


0.346 


. 54 


-0.104 


0.09 


-0.262 


0.06 


0.338 


0.44 


3.959 


0.02 


-0.024 


0,70 


-0.009 


0.91 


0.020 


0.74 


0.121 


0.29 


0.504 


0.66 

NO 
MONTHS ■" 


-0.971 


0.52 


WARM 


COLD MONTHS 


bl 


P>F 


M 


P>F 


0.00004 


0.72 


0.0003 


0.01 


-0.002 


0.82 


-0.012 


0.09 


-0.003 


0.69 


-0.012 


0.06 


-0.011 


0.33 


-0.004 


0.73 


-0.004 


0.55 


-0.003 


0.69 


0.001 


0.84 


-0.005 


0.34 


-0.00001 


0.71 


-0.0001 


0.05 


-0.00001 


. 76 


-0.0001 


0.05 


-0.003 


0.59 


-0.006 


0.42 


-0.001 


0,86 


-0.003 


0.53 


0.002 


0.97 


0,044 


0,76 


-0.005 


0.85 


0.029 


0.78 


-0.014 


0.49 


-0.004 


0.96 


-0.016 


0.09 


-0.037 


0.05 


0.004 


0.95 


0.389 


0.11 


-0.008 


0.37 


-0.045 


0.002 


-0.007 


0.50 


-0.046 


0.02 


0.16 


0.37 


-0,52 


0.05 



Table 5,9: Linear Regressions between the contin;uO'Us meteorolo'glcal variables 
and the SO, and NO. measurements for the Central Region 
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For the variables underlined in Tables 5.8 and 5.9 the ANCOVA routine was run 
(Table 5.10). 



Table 5,10j Results of the ANCOVA for the Central Region, 



VARIABLE 



P>F 



F 



P>F 



WARM MONTHS 



SO. 



type of precipitation 

sfc 48 hr sector 

interaction 

duration of precipitation 

sfc temperature 

airmass 

sfc 48 hr sector 

interaction 

duration of pcpn 

sfc temperature 



5.25 


0.026 


4.12 


0.001 


0,07 


0.998 


3.57 


0.065 


0.24 


0.627 


1.58 


0.166 


2.86 


0.015 


0.81 


0.643 


2.30 


0.136 


0.00 


0.990 



2.39 0.100 0,44 3.94 



1.54 0.100 0.50 3.77 



10, 



sfc 48 hr sector 
airmass 
interaction 
duration of pcpn 
sfc wind speed 

sfc 48 hr sector 
850 wind direction 
interaction 
duration of pcpn 
sfc wind speed 



2.80 


0.017 


0.50 


0.827 


0.54 


0.864 


3.59 


0.065 


0.02 


0.902 


6.49 


0.0001 


7.93 


0.0001 


2.74 


0.0089 


0.64 


0.4285 


0.08 


0.7747 



1.21 0.284 0.44 0.09 



4.44 0.0001 0.80 0.03 



cont'd 
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Table 5.10 cont'd 



COLD MONTHS 



SO. 



airmass 

sfc 48 hr sector 

Interaction 

sfc temperature 

track of low 
sfc 48 hr sector 
interaction 
sfc temperature 



5.11 


0.001 


4.59 


0.001 


5.00 


0.0001 


2.72 


0.107 


3.96 


0.005 


4.34 


0.001 


2,39 


0,010 


7,12 


0.011 



4.86 0.0001 0.79 1.50 



3.30 0.0001 0.67 2.35 



NO. 



sfc 48 hr sector 
nature of pcpn 
interaction 
cloud base height 
pcpn amount 

sfc 24 hr sector 
nature of pcpn 
interaction 
cloud base height 
pcpn amount 



6.65 


0.0001 


3,27 


0.019 


10.3 


0.0001 


1,58 


0.214 


5.12 


0,028 


2,76 


0.020 


1.38 


0.253 


0.59 


0.673 


0.75 


0.389 


8.13 


0.006 



6.32 0,0001 0.70 O.ll 



2.08 0.022 0.37 0.21 



from the AHCOVA analyses, the variables in Table 5.11 whose ANCOVA results are 
underlined in Table 5.10 provided the best set In explaining the variation in SO. 
and HO-. Here again, the set of variables chosen for the SO^ for the warm raonths 
was done with a view to being consistent with the results from the southwest 
region. 



SO. 



NO. 



WARM MONTHS 



COLD MONTHS 



WARM MONTHS 



COtD MONTHS 



precipitation type airmass sfc 48 hr sector sfc 48 hr sector 

sfc 48 hr sector sfc 48 hr sector 850 wind direction nature of pcpn 

pcpn duration sfc temperature pcpn duration cloud base height 

sfc temperature sfc wind speed pcpn amount 



Table 5.11: Variables which comprised the best set for the central region, as 
deduced from the ANCOVA, 



Table 5,12 contains the results of the AHCOVA whefi the variables which Here 
slgnif leant in the southwest analysis were used. 



?ARIABLI 



yARM MONTHS 



P>F 



SO. 



type of pcpn 
sfc 24 hr sector 
interaction 
precipi table water 
duration of pcpn 



4 


92 


0.032 


2 


.59 


0.026 





.42 


0.860 





.06 


0,807 


1 


.08 


0.304 



1.67 0.092 0.39 4.23 



NO, 



synoptic type 
sfc 48 hr range 
interaction 
duration of pcpn 
pcpn amount 



cloud type 

850 24 hr sector 

interaction 

cloud base height 

pcpn amount 



1.99 


0.090 


0.48 


0.841 


1.07 


0.413 


0.74 


0.395 


0.11 


0.743 




COLD MONTHS 




SO, 


0.56 


0.787 


1.11 


0.351 


2.87 


0.015 


3.75 


0.057 


1.31 


0.257 



1.07 0.412 0.48 0.09 



1.64 0.076 0.31 4.37 



NO, 



cloud type 

850 24 hr sector 

interaction 

cloud base height 

pcpn amount 



1.11 


0.369 


2.27 


0.087 


1.65 


0.149 


2.23 


0.141 


6.42 


0.014 



1,84 0,039 0.34 0.21 



Table 5.12: ANCOVA using the variables from the southwest analysis. 



■tt 



Gorapiring Table 5.12 with the results of the best set of variables in Table 5.11 
reireals that the amount of variance explained by using the variables that were 
best In the southwest region was much less than when variables determined by the 
analysis of the Central region were used. For SO, , 0.39 vs 0.44 and 0.79 vs 0.31 
were the differences in variance explained In tie warm aonths and cold months 
respectively. For N0„, the corresponding nunibers are 0.48 vs 0.80 and 0.34 vs 
0.70 (Table 5.13), 



WAIM MONTHS 



COLD MONTHS 



SO^ 
HO, 



SO^ 
NO, 



E (from 


ANCOVA) 


CE 


CE(SW) 


0.44 


0.39 


0.80 


0.48 


0.79 


0.31 


0.70 


0.34 



Table 5.13: The amount of variance of the dependent variable (SO. or N0«) 
explained by the ANCOVA analysis (R) for the central region when tie 
variables were selected from by an analysis of the central region 
(CE) and when the variables used in the analysis were the same as 
those used in the southwest analysis (CE(SW))* 
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5,4 Southeast Region 

fable 5.14 contains the results of the analysis of variance and Table 5.15 the 
results of the linear regression of the appropriate variables for the southeast 
region. The values umderlined indicate those variables which will be used in the 
ANCOVA. 

SO. 
WARM MONTHS COLD MONTHS 
F P>F F P>F 



Synoptic type 


1.20 


0.32 


0.26 


0.95 


Track of low 


2.64 


0.04 


1.46 


0.22 


Air mass 


1,49 


0.19 


1.47 


0.19 


850 nb 24 hr sector 


3.07 


0.02 


1.00 


0.41 


range 


1.47 


0.21 


0.92 


0.51 


48 hr sector 


4,28 
0.96 


0.001 


0.56 
0.49 


0.76 


range 


0.48 


0,89 


sfc 24 hr sector 


2.03 


0.07 


1.57 


0.16 


range 


0.62 


0.68 


0.86 


0.53 


48 hr sector 


3.15 
0.63 


0.008 


1.56 
0.94 


0.16 


range 


0.73 


0,49 


type of pcpn 


0.17 


0.68 


2.17 


0,14 


nature of pepm 


1.64 


0.21 


3.90 


0.01 


cloud type 


1.48 


0,20 


0.26 


0.97 


sfc wind dlr 


1.87 


0.09 


1,13 


0.36 


850 Bib wind dlr 


2.51 


0.05 


0.75 
NO 


0.56 
MONTHS 




WARM 


MOUTHS 


GOLD 




F 


P>F 


F 


P>F 


Synoptic type 


1,11 


0.37 


0.64 


0.70 


Track of low 


1.43 


0.24 


0.94 


0.44 


Alroass 


2.13 
1.89 


0.06 


0.50 
2.25 


0.83 


850 mb 24 hr sector 


0.11 


0.07 


range 


1.80 


0.12 


1.19 


0.32 


48 hr sector 


2.88 
1.73 


0.02 


1.26 
0.99 


0.29 


range 


0.11 


0.46 


sfc 24 hr sector 


1.59 


0.16 


1.29 


0.27 


range 


1.01 


0.42 


1.71 


0.13 


48 hr sector 


1.95 
0.51 


0.08 


0.50 
0.78 


0.83 


range 


0.82 


0.63 


type of pcpn 


. 01 


0,92 


1.11 


0.30 


nature of pcpn 


6.45 
1.35 


0.01 


4.67 
0.71 


0.01 


cloud type 


0.25 


0.67 


ifc wind dir 


2.72 


0,02 


0.94 


0.48 


850 mb wind dlr 


2.27 


0.07 


0.39 


0.82 


Table 5.14: Analysl 


i of Vi 


ir lance 


of the 


dlscr< 



43 



so, 





WARM 


MONTHS 


COLD MONTHS 


! 




bl 


P>F 


bl 


F>F 


CLOUD BASE HEIGHT 


0-001 


0.33 


0.001 


0. 


44 


SURFACE TEMPERATURE 


-0.012 
-0-002 


0.86 
0.98 


0.086 
0.83 


0. 


003 


SURFACE DEWPOINT 


0. 


003 


SURFACE WIND SPEED 


-0.152 


0.05 


-0.049 


0, 


24 


850 MB TEMPERATURE 


0.024 


0.71 


0.053 


0, 


09 


850 MB WIND SPEED 


0,023 


0.37 


-0.026 


0. 


.14 


CLOUD TOP 


-0.0003 


0-14 


-0.0002 


0. 


09 


CLOUD THICKNESS 


-0.0004 
-0.012 


0,11 
0.76 


-0,0003 
0,017 





.04 


PRECIPITABLE WATER 


0. 


.59 


TOTAL TOTALS 


-0.010 


0.72 


0.010 





.55 


RAINRATl mm/15 mln 


-O.lOf 


0.70 


0.037 





.93 


mm/30 min 


-0.147 


0.48 


-0.057 





.84 


mm/hr 


-0.195 


0,22 


-0.136 





.52 


DURATION OF PCPN 


-0,120 


0-14 


-0.144 





.12 


AVERAGE RAINRATE mm/li] 


r -0.527 


0-22 


-0.46 





.55 


AMOUNT OF PCPN 


-0.105 


0.16 


-0.070 





.17 


MAX 6 HI AMOUNT 


-0.08 


0.32 


-0,017 





.85 


7. STNS WITH PCPN 


0.707 


0.68 


-0.760 





.36 




NO 
WARM MONTHS 


COLD MONTHS 




bl 


P>F 


bl 




P>F 


CLOUD BASE HEIGHT 


0.0001 


0.38 


0.0002 




0.38 


SURFACE TEMPERATURE 


-0.017 


0.11 


-0.006 




0.51 


SURFACE DEWPOINT 


-0.016 


0.13 


-0.007 




0.45 


SURFACE WIND SPEED 


-0.02 


0.10 


-0.026 




0.04 


850 MB TEMPERATURE 


-0.015 


0.19 


-0.011 




0.26 


850 MB WIND SPEED 


0.005 


0.27 


-0.011 




0.03 


CLOUD TOP 


-0.0001 


0.04 


-0.00001 




0.80 


CLOUD THICKNESS 


-0.0001 


0.02 


-0.00001 




0,78 


PRECIPITABLE WATER 


-0.013 
-0.006 


0.04 
0.21 


-0.018 
-0,004 




0.05 


TOTAL TOTALS 


0.38 


RAINRATE mm/ 15 min 


-0.06 


0,20 


-0.015 




0.84 


mm/30 mla 


-0.05 


0,13 


-0.022 




0.67 


mm/hr 


-0.051 


0,05 


-0.032 




0.41 


DURATION OF PCPN 


-0.010 


0.44 


-0.021 




0.21 


AVERAGE RAINRATE mm/h 


r-0.147 


0.03 


-0.167 




0.22 


AMOUNT OF PCPN 


-0.023 
-0.009 


0.05 
0.47 


-0,058 
-0-055 




0.0004 


MAX 6 HR AMOUNT 


0.04 


7. STNS WITH PCPN 


0,105 


0,71 


-0.81 




0.001 



Table 5,15: Results of the linear regressions between the 
neteorologieal variables and the SOi and NO. readings. 



contlnyotis 



Table 5.16 contains the results of the ANCOVA for the southeast region and Table 
3-17 the results when the variables which were best for the southwest region were 
used. 
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Fp» 



Table 5.16s flit AHCOVA lor SQutheast region 

WARM MONTHS 





SO, 










f 


P>F 


F 


F>F 


R 


e2 


3.15 


0.026 


2,36 


0.013 


0.56 


3.88 


3,08 


0.012 










1.68 


0.154 










0,22 


0.640 










0.40 


0.534 










3.43 


0,020 
0,005 


3.06 


0.002 


0.69 


2,47 


3,72 










3.05 


0.010 










0.11 


0.741 










0.01 


0.941 










4.08 


0.003 


2,38 


0.014 


0.46 


4.56 


1.66 


0.191 










0.43 


0.787 










3.68 


0.062 










0,80 


0,377 










4.54 


0.001 


2.56 


0.010 


0.47 


4,10 


0.51 


0.728 










2.12 


0.133 










2.27 


0.140 










0.00 


0.973 
NO3 











track of low 
sfc 48 hr sector 
interaction 
duration of pcpn 
ifc wind speed 

850 wind direction 
sfc 48 hr sector 
Interaction 
duration of pcpn 
sfc wind speed 

850 48 hr sector 
track of low 
interaction 
duration of pcpn 
sfc wind speed 

850 48 hr sector 
850 wind direction 
interaction 
duration of pcptj 
sfc wind speed 



nature of pcpn 9,66 0.003 3,82 0.001 0.51 0.12 

sfc 48 hr sector 2.20 0.054 

Interaction 12.47 0.001 

amount of pcpn 3.23 0.080 

preclpitable water 1.20 0.280 

airmass 2.40 0.047 1.64 0.100 0.60 0.09 

sfc 48 hr sector 1.09 0.399 

interaction 1.61 0,156 

amoun t © f pcpn 1.99 0,16 9 

precipitable water 0.30 0.586 

850 48 hr sector 3.95 0.004 2.24 0.017 0,61 0.09 

airmass 1.73 0,136 

interaction 1.45 0.212 

amount of pcpn 3.07 0.089 



precipitable water 0.89 0.353 



cont'd 



4.1 



Tabli 5.11 



cost d 



850 48 hr sector 

nature of pcpn 
Interaction 

amount of pcpn 
preclpltable water 

nature of pcpn 
sfc 48 hr sector 
interaction 
amount of pcpn 
cloud thlcknesi 

airaiass 

sfc 48 hr sector 

Interaction 

amount of pcpn 

cloud thickness 

850 48 hr sector 
airmass 
interaction 
amount of pcpn 
cloud thickness 

850 48 hr sector 
nature of pcpn 
interaction 
amount of pcpn 
cloud thickness 



nature of pcpn 
sfc 48 hr sector 
interaction 
sfc temperature 

nature of pcpn 
sfc 24 hr sector 
Interaction 
sfc teaperature 

sfc 24 hr sector 
nature of pcpn 
interaction 
cloud thickness 
sfc temperature 



4.47 


0.001 


5.02 


0.0001 


0.5 


0.11 


6.47 


0.015 










5.36 


0.025 










7.96 


0.007 










3.62 


0.064 










3.54 


0.070 


2.64 


0.017 


0.49 


0.11 


1.59 


0.177 










9.92 


0.004 










3.61 


0.067 










0.83 


0.370 










1.05 


0.435 


0.98 


0.525 


0,58 


0.09 


0.88 


0.545 










1.03 


0.451 










1.55 


0.230 










0.33 


0.574 










2.42 


0.60 


1.82 


0,085 


0.63 


0.08 


1.37 


0.266 










1.60 


0,194 










3.33 


0.082 










1.21 


0.283 










2.66 


0.032 


3.34 


0.004 


0.50 


0.11 


2.71 


0.109 










4.76 


0.036 










7.46 


0.010 










2.46 


0.126 










« 


COLD MONTHS 










SO, 










3.88 


0.016 


2.09 


0,023 


0.54 


2.55 


1.10 


0.380 










1.71 


0.114 










7.54 


0.009 










5.46 


0.003 
0.017 


3.32 


0.0003 


0.58 


2.10 


2.77 










2.67 


0.014 










6.60 


0.011 










1.46 


0.221 


2.87 


0.005 


0.68 


1.81 


6,26 


0.002 










2.95 


0.013 










4.14 


0.051 










0.48 


0.496 











cont'd 



Table 5.16 



cont'd 



sfc 48 hr sector 
nature of pcpn 
interaction 
cloud thickness 
sfc tenperature 



1.67 


0.165 


3.85 


0.023 


2.35 


0.052 


3,12 


0.091 


1.85 


0.187 



2.35 0.025 0,67 2.02 



NO, 



850 24 lir sector 
nature of pcpn 
interaction 
pEicipi table water 
amount of pcpn 



4.43 


0.004 


5.13 


0.004 


2.16 


0.074 


3.41 


0.071 


18.1 


0.0001 



4.67 0.0001 0.56 0,19 



Table 5.17s The AUCOVA for southeast region with 
best in the southwest region. 



the variables which were 



WARM MONTHS 



SO. 



type of pcpn 
sfc 24 hr sector 
Interaction 
precipitable water 
duration of pcpn 



0.87 


0.358 


2.13 


0.074 


0.27 


0.893 


1.66 


. 206 


0.37 


0.547 



1.29 0.265 0.32 5.51 



synoptic type 
sfc 48 hr range 
interaction 
duration of pcpn 
amount of pcpn 



3.04 
1,90 
5.50 
0.14 
0.95 



0,022 
0.111 
0.0001 
0.712 
, 340 



COLD MONTHS 



3.89 0,0004 0,823 0.04 



SO, 



cloud type 

850 24 hr sector 

interaction 

cloud base height 

ainount of pcpn 



0.28 


0.959 


0.86 


0.495 


0.39 


0.855 


4.99 


0.030 


4.83 


0.032 



0.95 0.524 0,24 3.45 



MO, 



cloud type 

850 24 hr sector 

interaction 

cloud base height 

amount of pcpn 



0.82 


0.577 


1.39 


0.248 


0.67 


0.653 


0.07 


0.786 


12.50 


0.001 



1.51 0.123 0.33 0.25 



As in the other regions, the ANCOVA explained more of the vmriance in the SO. and 
NO. when the variables ysed were based on the ANOVA and regressions of the 
southeast region. The one exception was for the HO. for the warm aonths where the 
variables from the southwest region analysis (synoptic type, sfe 48 hr range, 
duration and amount of precipitation) performed better than those variables with 
the lower individual significance probabilities. This points up the non-objective 
nature of the process of choosing the best set of variables. Another possibility 
for the surprisingly good result with the variables from the southwest analysis In 
this case Is that the ANCOVA only considers cases for which all the variables are 
available. In using only this subset, some of the cases which contributed to high 
significance probabilities for some of the variables may have been eliminated in 
the ANCOVA for this particular set of variables. Table S.18 summarizes the 
results and Table 5.19 Indicates the variables which were found to have the best 
perfornance In the ANCOVA for this region. 



WARM MONTHS 



COLD MONTHS 



NO, 



SO^ 
NO, 



E (from 


ANCOVA) 


SE 


SE(S¥) 


0.69 


0.32 


. 61 


0.82 


0.58 


0.24 


0.56 


0.33 



Table 5,18! Tlie amount of variance of the dependent variable 
explained by the ANCOVA analysis (R). 



(SO, 



or 



NO^) 



NO. 



WARM MONTHS 



COLD MONTHS 



WARM MONTHS 



COLD MONTHS 



sfc 48 hr sector sfc 24 hr sector 
850 mb wind dlr nature of pcpn 



sfc 48 hr range 
synoptic type 



850 24 hr sector 
nature of pcpn 



pcpn duration sfc temperature 
sfc wind speed 



preclpl table water 
pcpn amount 



precipi table water 
pcpn amount 



Table 5.19: The variables which performed best tn the ANCOVA analysis for the 
southeast region. 
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5.4 Summary 

Table 5.20 summarizes the amount of variaace explaliied by the best analysis of 
covarianee for the 4 regions. 



WARM MONTHS 



COLD MONTHS 



no; 



NO, 





1 (from 


ANCOVA) 




il 


il 


SE 


m 


0.47 


0,44 


0.69 


0.83 


0.56 


0,80 


0.82 


0.81 


0.64 


0,79 


0,58 


0.61 


0.64 


0.70 


0.56 


0.65 



AVERAGE 



0.58 



• 6 8 



0.66 



0,73 



Table 5,20; Summary of the amount of variance explained by the best analysis of 
eovarlance for the 4 regions and the regional averages. 



the ANCOVA is In the reverse order to the 

the SO, and NO^ for the regions, I.e. the southwest region had the 

4 , . 3 



The amount ©f variance explained by 

averages of 

highest SO. and NO^ readings and the ANCOVA was able to explain the least amount 

of variance for this region. The pattern follows for the southeast, central and 

northwest regions respectively. Overall, the amount of variance explained 

averaged 0.66* 



An examination of the ANOVA and regressions of the 4 regions (Tables 4.4, 4.7, 
5,2, 5,3, 5.8, 5.9, 5.14, and 5,15) can shed some light on those meteorological 
variables which were the most consistently highly associated with the SO^ and NO^ 
measurements. Referring to the beginning of section 4 the best variables are 
these (the first 3 are categorical variables and the last 2 continuous variables); 



VARIABLE TYPE 



BEST VARIABLE 



synoptic situation 
airnass history 
characterization of pcpn 

ai rmas s char ac t er 1 s 1 1 c s 
rainfall 



airmass 

sfc 48 hr sector 

type of pcpn In the warm months and the 

nature of the pcpn in the cold months. 

surface temperature or preclpitable water 

pcpn amount or duration* 
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I. Testa to Predict Armual DepositioQ from a Selected Sample of Days 

The predictions of SO^ and NO- values from the ANCOVA using the best set of 
variables for each region for both the warm and cold months served as the basis 
for categorizing the days. Eight classes were used such that the bounds of the 
classes divided the predictions of SO and N0„ into 8 groups with approximately 
the same number of observations. With this classification in place for both SO 
and NO- for each day, region, and season, the following three experiments were run 
to see how well annual deposition statistics could be predicted from a selected 
sample of days. 

a) Non -Weighted Average 

For each categdry, the average deposition per day from the sample times the number 
of days of that particular category in the whole data set was calculated. These 
values for all the categories were then summed up to get the overall deposition. 
As a base line comparison, the average deposition per day from the selected number 
of days (regardless of the categorizations) times the number of days in the whole 
data set was also calculated. If the categorizations determined from the ANCOVA 
are useful in this application, the first calculation, referred to as prediction 
one, should perform better than the second (prediction two) at estimating annual 
deposition. 

b) Weighted Average 

This is similar to the above case a) except the sum of the deposition for the days 
in a given category in the sample is divided by the amount of rainfall received on 
those days and then multiplied by the total amount of precipitation received by 
that particular category within the whole data set. This reduces the variability 
introduced in the calculation of deposition from concentrations due to varying 
amounts of rainfall on the different days. 

e) Weighted Average for the Individual Stations 

Whereas the previous 2 methods dealt only with regional averages, here the 
depositions are calculated and compared with each of the 16 stations individually 
using the weighted average approach. 

The criteria for selecting the days for the test were that they had to include 
examples of all the categorizations, of all months of the year, and be contained 
generally in groups of at least 3 or 4 consecutive days. Also, the days had to be 
the same for all regions. These restrictions were necessary to simulate the 
process which would be carried out if results of the Eulerian model were used to 
predict annual deposition. The days selected were the following; 
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NO of DAYS 



January 11-17 


f 


February 8-13 


i 


March 8-18 


ti 


April 16-21 


« 


May 18-23 


6 


June 18-21,28,29 


6 


July 16-18,26-30 


8 


August 18-24 


7 


September 20-27 


8 


October 14-16 


3 


November 19-29 


11 


December 14-20 


f 



In all, 86 days were chosen 

iiven in Table 6.1 for each region are the number of days of each of the 8 
categories In the whole data set being considered and the number of days in the 
selected sample for each category. 



Table 6.1: The number of days, Iby category, for the data set being considered 
(first number) and for the selected sample (second number) for each 



region. 



NUMBIR OF DAYS IN EACH 

OF THE 8 CATEGORIES 
(WHOLE DATA SIT/ SAMPLE) 



SOUTHWEST REGION 







WARM 


MONTHS 






COLD 


MONTHS 


CAT 




^ 




HO 


h 




SO, 


NO 
10 3 


1 


5 


7 




9 


2 


7 


5 


8 




6 


11 


4 


11 4 


3 


6 


2 


7 




3 


11 


3 


11 5 


4 


7 


4 


8 




3 


11 


4 


12 4 


5 


6 


4 


8 




3 


10 


6 


11 4 


6 


6 


3 


8 




4 


12 


3 


10 5 


7 


6 


2 


7 




3 • 


10 


5 


10 3 



10 3 10 



coiit*'<l 



it 



Table 6.1 



cont'd 



CENTRAL REGION 







WARM 


MONTHS 






COLD 


MONTHS 


CAT 




^ 


NO 


% 




^ 


NO, 


1 


8 


7 


9 


2 


9 


4 


9 


5 


10 


2 


f S 


3 


8 


2 


7 


3 


9 


5 


7 


4 


9 


5 


i 


3 


11 


5 


9 3 


5 


8 


3 


3 


5 


9 





7 2 


6 


8 


2 


8 


3 


10 


2 


11 2 


7 


8 


4 


6 


1 


9 


1 


5 3 


8 


8 





8 


1 


10 


4. 


9 2 



SOUTHEAST REGION 







WARM 


MONTHS 




COLD 


MONTHS 


CAT 




SO, 


_N0^ 

7 




^ 


NO 

5 \ 


1 


6 


9 


2 


7 


3 


9 5 


8 


1 


8 2 


3 


7 


3 


5 1 


6 


4 


10 5 


4 


7 


4 


9 5 


11 


4 


9 2 


5 


5 


2 


6 5 


9 


2 


6 4 


6 


8 


4 


8 2 


10 


1 


9 2 


7 


6 


1 


6 1 


7 


1 


8 4 


8 


7 


4 


8 4 


9 


2 


8 1 



NORTHWEST REGION 







WARM 


MONTHS 


CAT 




SO, 


NO 
7 ^4 


1 


6 


2 


7 


2 


10 2 


3 


3 


5 


8 2 


4 


6 


1 


9 1 


5 


8 


1 


9 2 


6 


8 





10 2 


7 


6 


2 


8 1 


8 


8 


2 


9 2 



COLD 


MONTHS 


5 1 


NO^ 
6 -^0 


6 1 


6 i 


6 2 


7 1 


6 2 


6 I 


6 


6 2 


6 2 


7 


6 2 


6 2 


1 


7 2 



Table 6,2 contains the results of the first test of estimating the annual 
deposition using non-weighted averages for the 4 regions. The unit for deposition 
in all the following tables is mg m" . 



m 



OBSERVED PREDICTION PREDICTION NO OF DAYS 
DEPOSITION ONE TWO SAMPLE SET 



MONTHS 



REGION 



1699.535 
271.007 

1074.615 
241.371 

1391.583 
154.670 
702.781 
183.760 

1495.058 
212.401 

1010.660 
204.027 

646.121 

105.368 

245,217 

53.575 



2014.667 
299.023 

1050.385 
282.649 

1271-018 
183.356 
793.232 
188.870 

1770.407 
235.813 
693.914 
154.757 

437.727 

108.877 

153.094 

53.081 



2034.546 
316.928 

1059.523 
274.855 

1271.381 
160.352 
692.399 
180.715 

1540.353 
247.264 
814.526 
144,087 

450.752 
90.567 

165.797 
51.348 



27 
29 
32 
32 

23 
21 
25 

16 

23 
23 

18 
21 

15 

16 
10 
12 



50 


NO. 

SO 

NO^ 


WARM 


southwest 


60 


WARM 




84 


COLD 




85 


COLD 




66 


SO, 
NO 

NO3 


WARM 


central 


61 


WARM 




77 


COLD 




63 


COLD 




53 


NO^ 
NO3 


WARM 


southeast 


58 


WARM 




69 


COLD 




63 


COLD 




57 


SO, 

NO 

SO, 
NO^ 


WARM 


northwest 


70 


WARM 




47 


COLD 




51 


COLD 





Table 6,2: Deposition results using non-weighted averaging. 

The following Table 6,3 gives the same results as Table 6,2 using the weighted 
calculations. 



OBSERVED 


PREDICTION 


PREDICTION 


NO. OF 


DAYS 


MONTHS 


REGION 


DEPOSITION 


ONE 


TWO 


SAMPLE 


SET 








1699.530 


1808,050 


1756,010 


27 


50 


SO, 

NO. 

SO, 
NO^ 


WARM 


southwest 


271.016 


292,813 


285.101 


29 


60 


WARM 




1074.610 


1109.290 


993.714 


32 


84 


COLD 




241.373 


255,227 


260.318 


32 


85 


COLD 




1391.590 


977.407 


1104,570 


23 


66 


SO, 
NO^ 
SO 
NO3 


WARM 


central 


154.667 


154,793 


131,031 


21 


61 


WARM 




702.783 


815.510 


688.642 


25 


77 


COLD 




183.750 


210,161 


198.971 


16 


63 


COLD 




1495.060 


1826.030 


1681.020 


23 


53 


SO 
NO. 
SO, 
NO3 


WARM 


southeast 


212.409 


270,116 


269.026 


23 


58 


WARM 




1010.660 


957,081 


1032,380 


18 


69 


COLD 




204.032 


317,278 


215.674 


21 


63 


COLD 




646,136 


504,768 


527.749 


15 


57 


SO, 
HO 

so; 


WARM 


northwest 


105.357 


108,000 


102.626 


16 


70 


WARM 




245,207 


237,790 


244.141 


10 


47 


COLD 




53.576 


67.999 


72.145 


12 


51 


COLD 




Table 6.3: 


Deposition resu 


Its using w 


eighted j 


averae: 


Ine. 
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The followlag Table 6.4 eontalns the results of the weighted averaging for the 16 
Individual s tations . 



Tshli i.4j 


Deposition results using weigl 


OBSERfED 


PREDICTION 


PREDICTION 


m 


DEPOSITION 


ONE 


TWO S 


AMP] 


1120.540 


1275.235 


1120.000 


19 


203.858 


205.918 


206.802 


22 


849.895 


890.425 


802.802 


23 


115.988 


179.613 


172.599 


23 


1561.920 


1569.070 


1615.230 


21 


225.512 


224.273 


220,497 


22 


854.680 


747.379 


778.719 


23 


185.855 


185.654 


210.850 


23 


972.620 


1072.080 


1030.620 


19 


151.306 


178.553 


162.656 


20 


829.165 


859.081 


703.964 


16 


160.039 


168.656 


152.701 


16 


1642.630 


1765,010 


1721.430 


21 


213.167 


223.313 


213.948 


21 


899.555 


910.459 


907.509 


24 


197*701 


228.103 


241.801 


24 


1237.210 


1038.207 


1049.660 


19 


114.572 


132.278 


113.210 


17 


573.895 


542.270 


599.593 


21 


155.003 


203.160 


190.944 


14 


1048.310 


900.324 


809.637 


14 


120.818 


133.317 


114.105 


14 


645.620 


672.770 


561.736 


13 


156.832 


155.918 


126' 620 


8 


1131.220 


891.301 


861.826 


19 


110,286 


109.929 


96.522 


17 


437.295 


523.663 


526.940 


14 


104.343 


150.997 


136.616 


10 


1177.680 


1037.542 


976,166 


20 


122.757 


129.974 


111.289 


17 


540.810 


543.769 


574.829 


14 


141.332 


203.423 


173.889 


10 



SET 



32 


SO 
N0„ 


WARM 


38 


WARM 


54 


COLD 


53 


COLD 


42 


NO3 

SO4 

NO3 


WARM 


47 


WARM 


64 


COLD 


62 


COLD 


31 


^^3 
SO, 


WARM 


38 


WARM 


36 


COLD 


36 


COLD 


38 


SO, 

no;! 

sol 

NO3 


WARM 


41 


WARM 


56 


GOLD 


58 


COLD 


53 


so, 

SO^ 
NO3 


WARM 


47 


WARM 


66 
55 


COLD 
COLD 


37 


NO3 

SO4 
NO3 


WARM 


36 


WARM 


50 


COLD 


41 


COLD 


44 


SO, 

N03 

S04 

N03 


WARM 


39 


WARM 


42 


COLD 


34 


COLD 


49 


so 

NO^ 
SO, 


WARM 


41 


WARM 


46 


COLD 


39 


COLD 



STATION 



1011 



1021 



1031 



2011 



3011 



3021 



3031 



3041 



eont'd. 



## 



Table 6.4 


cont'd 


1412.450 


1668.980 


171.038 


226.004 


836.845 


809.845 


160.186 


249.280 


1230.930 


1247.744 


182.686 


228.152 


889.785 


923.838 


196.775 


293.937 


1032.180 


1153.450 


161.252 


164.850 


937.880 


865.137 


168.000 


158.890 


624.835 


833.982 


85.701 


93.476 


344.935 


190.102 


49.658 


5.684 


310.445 


174.372 


44.878 


48.744 


178.035 


230.351 


36.140 


62.445 


214.005 


40.229 


30.850 


8.091 


18.590 


16.363 


6.811 


5.561 


403.975 


349.248 


71.662 


69,266 


122.090 


84.096 


35.195 


35.462 


472.250 


584.340 


84.293 


95.144 


184.485 


226.876 


40.083 


57.591 



1522.970 


20 


207.359 


19 


909.082 


17 


173.856 


20 


1339.810 


17 


226.397 


17 


833.310 


12 


199.638 


13 


900.626 


17 


186.804 


18 


997.450 


11 


190.556 


13 


682.128 


15 


100.685 


16 


182.951 


5 


5.928 


5 


163.492 


5 


40.200 


5 


253,167 


6 


68.406 


7 


39.525 


2 


6,052 


3 


16.363 


3 


5.561 


5 


380.424 


11 


71.314 


11 


85.696 


7 


27.485 


10 


638.911 


12 


87.071 


12 


228.618 


3 


55,426 


6 



42 


NO3 

SO 

NO-j 


WARM 


43 


WARM 


57 


COLD 


52 


COLD 


35 


SO, 
NO* 


WARM 


38 


WARM 


39 


COLD 


35 


COLD 


33 


SO4 
NO 

NO3 


WARM 


37 


WARM 


50 


COLD 


44 


COLD 


32 


SO, 
NO^ 

sol 

NO3 


WARM 


33 


WARM 


33 


COLD 


29 


COLD 


25 


SO4 
NO 

sol 

NO3 


WARM 


32 


WARM 


24 


COLD 


29 


COLD 


12 


NO, 
SO 

N03 


WARM 


16 


WARM 


6 


COLD 


8 


COLD 


39 


so^ 

NO 
SO 


WARM 


42 


WARM 


28 


COLD 


35 


COLD 


40 


SO, 
NO 

^\ 
NO^ 


WARM 


46 


WARM 


25 


COLD 


30 


COLD 



4011 



4011 



4031 



4041 



6051 



6061 



6071 



6081 



From Tables 6.2 to 6.4 the absolute fractional error (defined as the absolute 
value of the difference between the observed and predicted annual deposition 
divided by the observed annual deposition) was calculated for the 2 types of 
predictions (with and without the categorizations being used) from the 3 tests 
performed. The results are displayed in Tables 6.5 to 6.13. Shown in the final two 
columns are the number of times one prediction was better than the other, i,e. a 
smaller absolute fractional error. Table 6,5 gives the results for all the data, 
Tables 6.6 to 6.9 divide the data by regions, Tables 6.10 and 6.11 stratify the 
data by warm and cold months and Tables 6.12 and 6.13 divide the data depending on 
the number of observations in the selected sample. 



IS 



ABSOLUTE FRACTIONAL EREOl 
SO NO- 

PREDICTION PREDICTION PREDICTION PEEDICTION 
ONE TWO ONE IW 



BEST PREDICTION 
ONE TWO 



NON WEIGHTED 


0.20 


0.15 


0.11 


0.13 


i 


10 


WEIGHTED 


0.13 


0.08 


0.18 


0.18 


J 


11 


STATIONS 


0.17 


0.17 


0.22 


0.20 


m 


29 



Table 6.5i Deposition errors for all the data 



ABSOLUTE FRACTIONAL ERROR 
SO, NO. 

PREDICTION PREDICTION PREDICTION PREDICTION 
ONE TWO ONE TWO 



BEST PREDICTION 
ONE TWO 



NON WEIGHTED 


0.10 


0.11 


0.14 


0.15 


2 


WEIGHTED 


0.05 


0.05 


0.07 


0.07 


2 


STATIONS 


0.07 


0.06 


0.07 


0.07 


f 



Table 6.6: Deposition errors for the southwest region. 



ABSOLUTE FRACTIONAL ERROR 
SO, NO^ 

PREDICTION PREDICTION PREDICTION PREDICTION 
ONI mB ONE Hi 



BEST PREDICTION 

ONE "'tis 



NON WEIGHTED 0.11 
WEIGHTED 0,23 
STATIONS 0.12 



0.05 

0.11 
0.15 



0.11 
0,07 
0.19 



0.03 
0.12 
0.16 



Table 6. 7 1 Deposltldn eftofs for the centfal tegioti. 



ABSOLUTE FRACTIONAL ERROR 
SO, NO^ 

PREDICTION PREDICTION PREDICTION PREDICTION 

ONE TWO ONE li© 



BEST PREDICTION 

ONE XiO 



NON WEIGHTED 


0.25 


0,11 


0,18 


0.23 


2 


WEIGHTED 


0.14 


0.07 


0.41 


0.16 





STATIONS 


0.11 


0.13 


0.33 


0.24 


8 



§_ 



Table 6.8: Deposition errors for the southeast region. 



5fi 



ABSOLUTE FRACTIONAL EtlOl 
SO NO 

PREDICTION PREDICTION PREDICTION PlEDIGTION 

ONE TWO ONE TIS 



BEST PREDICTION 
ONE TWO 



NON WEIGHTED 


0.35 


0.31 


0.02 


0.09 


i 


WEIGHTED 


0.12 


0.0f 


0-15 


0.10 


1 


STATIONS 


0.32 


0.35 


0-29 


0.33 


it 



Table 6.9: Daposition errors for the northwest regioti. 



ABSOLUTE FRACTIONAL ERROR 
SO, NO 

PREDICTION PREDICTION PREDICTION PREDICTION 
ONE TWO ONE TWO 



BEST PREDICTION 
ONE TWO 



NON WEIGHTED 


0.19 


0.15 


0.11 


0.13 


i 


4 


WEIGHTED 


0.2Q 


0.14 


0.09 


0.21 


% 


6 


STATIONS 


0.20 


0.19 


0.14 


0.13 


i:f 


15 



T«hl© I, 10 I Beposltlon errors for the warm months. 



ABSOLUTE FRACTIONAL ERROR 
SO NO 

PREDICTION PREDICTION PREDICTION PREDICTION 
ONE TWO ONE TWO 



NON WEIGHTED 0.21 
WEIGHTED 0.07 
STATIONS 0.13 



0,14 
0.03 
0,16 



0.11 
0.26 
0.30 



0.12 
0.14 
0.26 



BEST PR] 

ill 


SDICTION 
TWO 


1 

i 
m 


i 

s 
m 



Table 6.11: Deposition errors for the cold momthi. 



ABSOLUTE FRACTIONAL ERROR 
SO NO^ 

PREDICTION PREDICTION PREDICTION PREDICTION 
ONE TWO ONE TWO 



BEST PREDICTION 

.'Oil TWO 



NON WEIGHTED 


0.25 


0.18 


0.10 


0,11 


i 


1 


WEIGHTED 


0.15 


0.09 


0.20 


0.14 


1 


f 


STATIONS 


0.22 


0.25 


0.36 


0.30 


m 


1* 



Table 6,12s Deposition errors when the ratio of sampled to total obs eryatloiis 
were less than 0,35, 



if 



NON WEIGHTED 


0.13 


0.08 


0.13 


0.16 


3 


WEIGHTED 


0.11 


0.08 


0.14 


0.13 


2 


STATIONS 


0-11 


0.10 


0.10 


0.11 


19 



ABSOLUTE FRACTIONAL ERROR 
SO. NO 

PREDICTION PIEDICTION PREDICTION PREDICTION BEST PREDICTION 
ONE TWO ONE TWO ONE TWO 
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Table 6.13: Daposition errors when the ratio of sampled to total observations 
were less than 0.35. 



from the previous 9 tables the following generalizations can be made: 

i) The weighted averaging performed better for the SO, depositions but not 
for the NO. depositions (Table 6.5). 

11) Predictions using the categorizations performed marginally better for 
the individual stations (35 to 29) while predictions not using the 
categorizations were better for the regional averages » both weighted and 
non-weighted (Table 6.5). 

lii) For the individual stations, the southwest region was the best while the 
northwest region was the poorest in estimating annual deposition (Tables 
6.6 to 6.9). 

Iv) By season, SO, estimates were better In the cold months while the HO. 
estimates were better in the warm months (Tables 6.10 and 6.11), 

v) The estimates improved with the number of cases in the sample. When the 
selected sample had less than one-third of the total observations, the 
errors were around 207. while they were around 10% when more than one- 
third of the days were in the sample (Tables 6.12 to 6.13). 

Overall however, there is no compelling evidence that the categorizations based 
upon the ANCOVA for each region are useful In estimating annual deposition. The 
variabilities within the categories are too large. 



7, Summary and Conclusions 

Meteorological data were analysed on those days In 1982 when SO^ or MO^ 



measurements were recorded by the any of the 16 APIOS st 
divided into 4 regions (southwest, central, southeast a 



at ions. The province was 
nd northwest) to reflect 



the grouping of the APIOS stations. Data used in the meteorological analysis 
included surface weather maps from the Ontario Weather Centre (every 6 hours)), 
U.S. Weather Bureau surface and 500 rab maps (daily), CMC 850 mb maps (every 12 
hours), representative surface hourly reports and upper air soundings for the 
particular precipitation event, and cliraatologlcal precipitation data, both dally 
and hourly when available. Fields extracted were the synoptic situation, the track 
of the controlling low, the airmasses involved, 24 and 48 hr surface and 850 mb 
trajectory information (sector and range), various surface and upper air 
temperatures and winds, and intensity amount and duration variables of the area- 
wide precipitation. 

The relationships between the SO, and N0„ readings for the southwest region were 
examined. For the discrete variables (Table 4.3c), the trajectories associated 
with high SO and NO^ measurements were similar but for the other meteorological 
fields, the important factors were dissimilar. Synoptic conditions had some effect 
on the high SO, readings but not on the high NO^ readings. The airmasses 
associated with them were also different, mP and mT air with SO , and cA cmA and 
mA with N0„. High SO, readings were associated with convective situations and 
rainfall, while high N0„ readings were associated with snow but had no 
relationship to the type of precipitation (convective or continuous). 

For the continuous variables, the general trend (Table 4.8) is that high SO^ 
readings were associated with high temperatures and moisture and low precipitation 
events, low cloud bases, and low windspeeds. For NO^, high readings were 
associated with low temperatures and moisture, low wind speeds, low cloud bases, 
and low precipitation events. 

The data were subdivided into 6 warm months (April to September) and 6 cold montbs 
to reduce the variability of the SO. and NO^ readings while still maintaining 
enough data in each set to perform statistical tests. An analysis of variance was 
performed between the SO, and NO^ measurements and the discrete meteorological 
variables while linear regressions were done between the SO, and NO^ readings and 
the continuous meteorological variables to search for significant variables. Those 
variables which showed promise individually were then combined using the analysis 
of covariance technique to try to explain as much of the variance In the SO^ and 
10^ readings as possible. 

tests indicated that better results were obtained by examining each region 
separately in determining the best set of variables. Different meteorological 
conditions among regions and the different observing networks in each region (both 
meteorological and APIOS stations) with regard to the type and density, and 
spatial distribution of stations were the underlying causes. From the ANCOVA the 
following were the variables which were deemed the best set in each region. 



SOIJTHWEST REGION 



SO. 



NO. 



lAiH MONTHS 



GOLD MOIfTHS 



ilAlM MONTHS 



COLD MONTHS 



precipitatioTi type cloud type synoptic type 

sfc 24 hr sector 850 24 hr sector sfc 48 hr range 

precipitable water cloud base height duration of pcpn 

duration of pcpn amount of pcpn amount of pcpn 



cloud type 

850 24 hr sector 

cloud base height 
anount of pcpn 



CENTRAL REGION 



SO. 



WARM MONTHS 



COLD MONTHS 



NO, 



WARM MONTHS 



COLD MONTHS 



precipitation type airmass sfc 48 hr sector sfc 48 hr sector 

sfc 48 hr sector sfc 48 hr sector 850 wind direction nature of pcpn 

pcpn duration sfc temperature pcpn duration cloud base height 

sfc temperature sfc wind speed pcpn amount 



SOUTHEAST REGION 
SO. 
WARM MONTHS 



COLD MONTHS 



WARM MONTHS 



COLD MONTHS 



sfc 48 hr sector sfc 24 hr sector 
850 mb wind dlr nature of pcpn 



sfc 48 hr range 
synoptic type 



850 24 hr sector 
nature of pcpn 



pcpn duration 
sfc wind speed 



sfc temperature 



precipitable water precipitable water 
pcpn amount pcpn amount 



NORTHWEST REGION 
SO. 
MAIM MONTHS 



airmass 

ifc 48 hr range 



COLD MONTHS 

sfc 24 hr range 
synoptic type 



tO> 
WAIM MONTHS 

airmass 

sfc 43 hr sector 



COLD MONTHS 

airmass 

850 48 hr sector 



cloud base height 
sfc tempe r a tur e 



sfc temperature sfc temperature amount of pcpn 
precipitable water stations with pcpn 



In general, the ANCOVA with these variables was able to explain about two-thirds 

of the variance in the SO, and NO., measurements. 

4 3 

From the ANOVA and the linear regressions performed on the 4 regions, the 
following meteorological variables were the most consistently strongly associated 
with the SO and NO variations for the respective type of variable. 



synoptic sltuatloTi alrniass 
alrmass history sfc 48 hr sector 

characterization of pcpa type of pcpn in the warm months and the 

nature of the pcpn In the cold months, 
airmass characteristics surface temperature or preclpltable water 
rainfall pcpn amount or duration. 

For more conclusive results In future experiments or studies in statistical 
modelling of deposition amounts, these are the meteorological factors which should 
be considered and designed Into the data collection efforts. 

The predictions of SO. and NO. based on the ANCOVA were then used to classify the 
days in an effort to estimate annual deposition from a selected sample of days. 
These results showed that annual deposition could be estimated to within 107. to 
207, when approximately one-third of the days were sampled' . But the categorizations 
used were not significantly better than treating all the days the same in this 
regard I.e. stratifying the days sampled by categories showed no significant 
improvement. The variation in daily deposition within the categories is too large. 

To Improve upon this situation a number of problems associated with the data must 
be rectified. The most important are: 

1. The precipitation events need to be better and more uniformly defined in 
their timing. Drastically different times for the beginning and end of 
precipitation events at neighbouring APIOS stations were noted. 

2. The aeteorologleal observations should be taken at the same locale as 
the APIOS sites. Differences in distances from the APIOS sites to the 
closest AES hourly station among the 4 regions (and the different hours 
of operations of such stations) added noise to the data set, 

3. Better time resolution on the meteorological and precipitation chemistry 
observations are necessary. During significant rainfalls, the synoptic 
situation is often rapidly changing so that any depiction over 24 hrs 
can be misleading. 

4. More detailed meteorological cloud observations (e.g. types and bases) 
would be useful. The analysis in the southwest region showed these to be 
important but the other regions, where these observations may not have 
been as good or as representative, did not support this. If several 
years of precipitation chemistry measurements are available, it may be 
worthwhile to compare the precipitation chemistry (perhaps concentration 
contours by month or season) to the similar cloud and precipitation 
climatology stratified by frequency, nature, type, and/or amount. 
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Appendix A. REGIONS FOR MOE APIOS CLIMATOLOGY 



Surface Stations 



Southwest em Region 

1. APIOS StatloQd 

- Melbourne 

- Longwoods 

- North Eaathope 

- Welles ley 

2. Hourly and Synoptic Stations 

- Windsor (YQG) 

- Sarnla (YZR) 

- London (YltJ) 

- Slmcoe (WMK) 

« Haallton (YHM) 

- Toronto International Airport (YYZ) 

- Toronto Island (YTZ) 

- Mount Forest (WMN) 

- Kitchener (TKF) 

- St. Catherines (YSN) 

- Goderlch (WGD) 



3. CliBi«tolQ|leal Stations s districts 613 ,614 , 612 amth ©f 44,5'N, and 613 
vest of 80 'W or south of 44*N. 



Cgmtgal Regloti 

I 

1. AflOS Stations 

- Raven Lake 
» lalsani Lake 

- Nithgrove 

- Dorset 

2. Hourly and Synoptic Stations 

- WiartOE (YW) 

- Jfuskoka (YQA) 

- §0Tm Bay (YZl) 

- Sodbury (YSB) 

- Horth Bay (YYB) 

- Peterborough (YPQ) 

3. Glimatological Stations: districts 609, 611 , 616, 612 iiorth of 44.5°N, 
606 south of 46.75'N, 608 south of 46. 75"'^ and 615 between 78'¥ and 
79'¥ and north of 44®Ef between 79''¥ and 80°W, 



Southeastern Region 

1, API OS Stations 

- yhltman Creek 

- Eailton 

- Charleston Lake 
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- arahaa Laka 

2. iourly and Symoptle Statloaa 

- TrtntoQ (m) 

- Kingston (YGK) 

- Ottawa lateraatlonal Airport (YOW) 

- Petawawa (Y¥A) 

3. Cllaatologlcal StatloQs: districts 610, aod 615 east of 78'W. 



Northwestern Region 

1. AFIOS Stations 

- Forbes Township 

- Quetico Centre 
»■ Lac La Croix 

- Fernberg 

2. Hourly and Synoptic Stations 

- Thunder Bay (YQT) 

- Armstrong (YYW) 

• Sioux Lookout (YXL) 

- Drydan (YHD) 

- Atlkokan (WIS) 

- International Falls, U.S. (INL) 

- Hibbing, U.S. (HIB) 

3. Cllmatologlcal Stations: district 604 between 88*W and 93®y and south of 
50.5*N. 



tipper Air Stat ions 



Horth Reg too 

1. iQCernatioBal Falls ( 

2. Trout Lake (WTL) 



South Region 

1. FUmt (FNt) 

2. Buffalo CBUF) 

3. Manlwaki (YHW) 

4. SauLt St© Marie (SSM) 
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Appendix B: SUMMARY OF THE DATA TO BE COLLECTED 

1. Synoptic type 

- waEm froQt 

- cold f^oat 

• stationary front 
^ occluded low 

• cold low 
« alraass 

2. Track of eontrolllng low 

« Alberta low * 

- Colorado low 

- Texas low coming north of the Appalachians 

- Texas low going south of the Appalachians 

- Low forming in the Great Lakes region 

3. AlrniassCts) Involved 

- maritime Tropical air 

- maritime polar air 

- maritime arctic air 

f cold maritime arctic air 
" continental arctic air 

4. Sector (octants) of air 24 and 48 hours previous at surface and 850 mb. 

5. Distance of the origin of the 24 and 48 hour trajectories from the main 
APIOS observing site. 
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6. 850 ab. tinptratura 

7. 330 nb. wind 

3. Pteelpltablt tfater 

9. Stability Index 

10. Cloyd base height 

11. Cioiid top 

12. CloMd thickness 

13. Surface teiapeiratpse 

14. Surface dewpoint 

15. Surface wind 

16. Predominant type of cloud 

17. Nature of precipitation 

- rain 

- snow 

- ice 

- freezing precipitation 

18. type of precipitation 

- intermlttant 

- continuous 

- showery 

19. Average rainfall rate (mm/hr) 

20. Maxinun rainfall rate (mm/hr, aim/30 mln, mm/15 mln) 

21. Average duration of precipitation event 
22* Average precipitation amount 

23. Percentage of stations receiving precipitation 

24. Average snowfall amount 

25. Average snowfall rate (co/6 hr) 



67 



Appendix C 



The following la a summary of the data available for the 4 regions (Soothwest 
(SW), Central (CE), Southeast (SE), and Northiiest (!W)). Sulfate, nitrate, and 
AFIOS rainfall data are averages from the APIOS stations which reported 
precipitation on a particular day. For the discrete variables the first number Is 
the number of observations of a given category followed by Its percentage of the 
available observations In brackets. For the continuous variables, the mean and 
standard deviation (SD) is given. Missing data (indicated by MSG) refers to either 
the fact that data were not available, a reliable determination of the variable 
could not be made with the data available (e.g. predominant cloud type), or there 
were no data (e.g. no rainfall reports). 



Table CI: A sunmary of the meteorological data collected. 



TOTAL NO OF OBS 



183 



189 



SE 
151 



NW 
161 



SYNOFTIC TYPE 



WARM FRONT 


40 (22) 


38 (20) 


30 (20) 


28 (17) 


COLD FRONT 


47 (26) 


42 (22) 


32 (21) 


18 (11) 


STATIONARY FRONT 


16 (9) 


10 (5) 


8 (5) 


15 (9) 


OCCLUDED LOW 


40 (22) 


47 (25) 


40 (26) 


61 (38) 


COLD LOW 


8 (4) 


13 (7) 


10 (7) 


12 (7) 


AIRMASS 


21 (11) 


20 (11) 


15 (10) 


13 (8) 


WARM AND COLD FRONT 11 (6) 


19 (10) 


16 (11) 


14 (9) 



TRACK OF LOW 



ALBERTA LOW 


55 


(30) 


61 (32) 


45 (30) 


106 (66) 


COLORADO LOW 


85 


(46) 


88 (47) 


70 (46) 


51 (32) 


TEXAS LOW NORTH 


13 


(7) 


7 (4) 


10 (7) 


2 (1) 


TEXAS LOW SOUTH 


11 


(6) 


8 (4) 


11 (7) 





GREAT LAKES 


19 


(10) 


22 (12) 


14 (9) 


1 (1) 


HUDSON BAY LOW 







3 (2) 


1 (1) 


1 (1) 



US 



SI 



m 



AIRMASS 










CA 


33 (18) 


25 (13) 


19 (13) 


29 (18) 


CMA 


3 (2) 


3 (2) 


2 (1) 


5 (3) 


MA 


13 (7) 


27 (14) 


23 (15) 


26 (16) 


MP 


4 (2) 


5 (3) 


4 (3) 


2 (1) 


MT 


2 (1) 


1 (1) 


1 (1) 





CA CMA 


6 (3) 


12 (6) 


7 (5) 


11 (7) 


CA MA 


35 (19) 


3 2 (17) 


30 (20) 


26 (16) 


CMA MA 


17 (9) 


15 (8) 


14 (9) 


10 (6) 


MA MP 


50 (27) 


56 (30) 


40 (26) 


40 (25) 


MP MT 


20 (11) 


12 (6) 


11 (7) 


12 (7) 


MSG 


.0 


1 


Q 






850 MB TRAJECtORT 



24 HR SECTOR 



Nm? 


16 


(9) 


12 (7) 


8 (5) 


25 (16) 


mvj 


47 


(26) 


39 (21) 


25 (17) 


35 (23) 


¥S¥ 


72 


(40) 


85 (46) 


67 (45) 


39 (25) 


ssw 


37 


(21) 


38 (21) 


43 (29) 


33 (22) 


SSE 







5 (3) 


3 (2) 


9 (6) 


ESE 


4 


(2) 


4 (2) 


4 (3) 


1 CD 


ENE 


2 


(1) 


1 (1) 





2 (1) 


ME 


1 


(1) 








9 (6) 


MSG 


i 




5 


"1- 


1 


24 HI RANGE (KM) 












<200 


5 


(3) 


9 (5) 


6 (4) 


20 (13) 


200-400 


42 


(23) 


28 (15) 


17 (11) 


37 (24) 


400-600 


30 


(17) 


36 (20) 


31 (21) 


45 (29) 


600-800 


35 


(19) 


47 (26) 


35 (23) 


36 (24) 


800-1000 


28 


(16) 


25 (14) 


28 (19) 


8 (5) 


1000-1200 


16 


(9) 


17 (9) 


13 (9) 


3 (2) 


1200-1400 


17 


(9) 


17 (9) 


14 (9) 


3 (2) 


1400-1600 


5 


(3) 


1 (1) 


3 (2) 





1600-1800 


2 


(1) 


3 (2) + 


1 (1) 


1 (1) 


1800-2000 







1 (1) 


2 (1) 





MSG 


3 




5 


1 


S 



69 



ft SI Ht 



A8 HR SECTOR 



NKW 


16 


(10) 


IS 


(10) 


6 


(4) 


21 (15) 


imw 


51 


(31) 


48 


(27) 


30 


(21) 


41 (30) 


tfSW 


63 


(38) 


69 


(39) 


67 


(48) 


32 (23) 


ssw 


28 


(17) 


34 


(19) 


28 


(20) 


31 (22) 


SSE 


4 


(2) 


4 


(2) 


2 


(1) 


7 (5) 


ESE 







2 


(1) 


4 


(3) 


2 (1) 


ENE 


1 


(1) 







2 


(1) 





ME 


1 


(1) 







1 


(1) 


4 (3) 


MSG 


19 




14 




11 




23 


48 HI RANGE (KM) 
















<200 


3 


(2) 


2 


(1) 


1 


(1) 


2 (1) 


200-400 


7 


(4) 


6 


(3) 


4 


(3) 


11 (8) 


400-600 


11 


(6) 


13 


(7) 


7 


(5) 


22 (16) 


600-800 


21 


(12) 


16 


(9) 


12 


(8) 


32 (23) 


800-1000 


21 


(12) 


21 


(12) 


21 


(14) 


21 (15) 


1000-1200 


23 


(13) 


24 


(13) 


18 


(12) 


22 (16) 


1200-1400 


30 


(17) 


23 


(13) 


26 


(18) 


13 (9) 


1400-1600 


28 


(16) 


21 


(12) 


14 


(10) 


6 (4) 


1600-1800 


19 


(11) 


26 


(15) 


19 


(13) 


7 (5) 


1800-2000 


12 


(7) 


19 


(11) 


10 


(7) 


2 (1) 


2000-2200 


4 


(2) 


7 


(4) 


10 


(7) 





2200-2400 


Q 




1 


(1) 


3 


(2) 





MSG 


4 




10 




6 




23 


SIIIFACE TIAJECTORY 
















24 HR SECTOR 

















ssv 

SSE 

ESE 
ENE 
KNE 
MSG 



12 


(7) 


14 (8) 


10 (7) 


18 (12) 


21 


(13) 


15 (9) 


10 (7) 


9 (6) 


30 


(19) 


37 (21) 


26 (18) 


a (5) 


44 


(27) 


37 (21) 


46 (33) 


33 (22) 


14 


(9) 


32 (18) 


17 (12) 


36 (24) 


15 


(9) 


20 (12) 


23 (16) 


18 (12) 


11 


(7) 


9 (5) 


5 (4) 


13 (9) 


14 


(9) 


10 (6) 


4 (3) 


12 (8) 


22 




15 


10 


m 



70 



Stf Si SE NV 

24 HR RANGE (KM) 



<200 


28 


(17) 


27 


(15) 


19 (13) 


22 (15) 


200-400 


43 


(27) 


36 


(21) 


31 (22) 


48 (32) 


400-600 


40 


(25) 


48 


(27) 


46 (32) 


38 (26) 


600-800 


24 


(15) 


34 


(19) 


23 (16) 


26 (17) 


800-1000 


16 


(10) 


19 


(11) 


8 (6) 


6 (4) 


1000-1200 


3 


(5) 


7 


(4) 


14 (10) 


8 (5) 


1200-1400 


3 


(2) 


4 


(2) 


2 (1) 


1 (1) 


MSG 


21 




14 




8 


12 


SURFACE TRAJECTORY 














48 HI SECTOR 














NNW 


20 


(15) 


16 


(10) 


12 (9) 


16 (12) 


WNW 


16 


(12) 


28 


(17) 


28 (21) 


8 (6) 


wsv 


28 


(20) 


28 


(17) 


22 (17) 


7 (5) 


ssw 


26 


(19) 


34 


(21) 


28 (21) 


27 (21) 


SSE 


12 


(9) 


17 


(10) 


15 (11) 


28 (22) 


£S£ 


10 


(7) 


17 


(10) 


13 (10) 


15 (12) 


KNE 


14 


(10) 


9 


(6) 


5 (4) 


16 (12) 


NNE 


12 


(9) 


14 


(9) 


8 (6) 


12 (9) 


MSG 


45 




26 




20 


32 


48 HR RANGE (KM) 














<200 


9 


(6) 


9 


(5) 


6 (4) 


14 (11) 


200-400 


25 


(18) 


23 


(14) 


12 (9) 


32 (24) 


400-600 


27 


(19) 


35 


(21) 


36 (27) 


27 (21) 


600-800 


31 


(22) 


28 


(17) 


18 (13) 


25 (19) 


800-1000 


24 


(17) 


20 


(12) 


28 (21) 


11 (8) 


1000-1200 


7 


(5) 


30 


(18) 


17 (13) 


16 (12) 


1200-1400 


5 


(3) 


13 


(8) 


13 (10) 


3 (2) 


1400-1600 


5 


(3) 


6 


(4) 


3 (2) 


3 (2) 


1600-1300 


3 


(2) 


2 


(1) 


2 (1) 





1800-2000 


2 


(1) 







1 


i 


2000-2200 


1 


(1) 













MSG 


44 




23 




16 


30 


TYPE OF PRECIPITATION 


88 


(54) 


68 


(41) 


73 (52) 




CONTINUOUS 


71 (51) 


CONVECTIVE 


74 


(46) 


99 


(59) 


68 (48) 


68 (49) 


MSG 


21 




22 




10 


22 



71 



sv 



CI 



SE 



NATURE OF PRECIPITATION 



RAIN 


99 


(61) 


103 (62) 


95 


(67) 


95 (68) 


SHOW 


54 


(33) 


58 (35) 


35 


(25) 


36 (26) 


FREEZING 


2 


(1) 


2 (1) 


3 


(2) 





MIXED 


7 


(4) 


4 (2) 


8 


(6) 


8 (6) 


MSG 


21 




22 


10 




W 


PREDOMINANT CLOUDS 














TOWERING CUMULUS 


3 


(2) 


11 (7) 


5 


(4) 


11 (9) 


ALTOCUMULUS 


2 


(1) 


2 (1) 


2 


(1) 


2 (2) 


ALTOSTRATUS 


4 


(3) 


3 (2) 


1 


(1) 


2 (2) 


CUMULONIMBUS 


21 


(13) 


11 (7) 


18 


(13) 


21 (17) 


CUMULUS 


14 


(9) 


22 (14) 


8 


(6) 


12 (9) 


NIMBOSTRATUS 


29 


(18) 


12 (8) 


12 


(9) 





STIATO CUMULUS 


70 


(45) 


82 (51) 


77 


(56) 


60 (47) 


STRATUS 


14 


(9) 


17 (11) 


15 


(11) 


19 (15) 


MSG 


26 




29 


13 




M 


SUIFACE WIND SECTOR 














NNM 


16 


(10) 


8 (4) 


12 


(8) 


16 (11) 


WNW 


31 


(18) 


32 (18) 


11 


(7) 


18 (12) 


WSW 


28 


(17) 


28 (16) 


23 


(15) 


16 (11) 


SSW 


25 


(15) 


25 (14) 


28 


(19) 


13 (9) 


SSE 


26 


(15) 


34 (19) 


20 


(13) 


9 (6) 


ESE 


16 


(9) 


22 (12) 


15 


(10) 


21 (14) 


ENE 


14 


(8) 


10 (6) 


25 


(17) 


47 (31) 


NNE 


11 


(7) 


13 (7) 


16 


(11) 


10 (7) 


CALM 


2 


(1) 


6 (3) 







X 


MSG 


14 




11 


1 




M 


850 MB WIND SECTOR 














NNW 


15 


(9) 


22 (14) 


3 ( 


[2) 


18 (14) 


WNW 


40 


(25) 


40 (25) 


26 


(19) 


28 (22) 


WSW 


60 


(38) 


63 (40) 


75 


(55) 


24 (18) 


SSW 


24 


(15) 


21 (13) 


27 


(20) 


28 (22) 


SSE 


8 ( 


:5) 


5 (3) 


3 ( 


[2) 


16 (12) 


ESE 


4 1 


:3) 


3 (2) 







6 (5) 


ENE 


5 ( 


:3) 


4 (2) 


2 


[1) 


3 (2) 


NNE 


2 ( 


'D 










8 (5) 


MSG 


25 




31 


15 




ja 



72 



'■^'% 



^ CE SE NW 
CLOUD BASE HEIGHT (M) 



MEAN 


666 


658 


718 


749 


SD 


500 


498 


431 


516 


MSG 


26 


If 


1:4 


39 


SURFACE TEMPERATURE CO 










MEAN 


5.5 


6.5 


6.9 


5.8 


JD 


10.7 


10.0 


9.4 


11.0 


mm 


15 


10 


1 


10 


SURFACE DEWPOINT CO 










MEAN 


3.4 


3.9 


5.1 


3.3 


SD 


10.7 


10.3 


9.3 


11.2 


MSG 


15 


10 


1 


10 


SURFACE WIND SPEED (m/s) 










• MEAN 


9.2 


8.4 


8.7 


l.i 


SD 


5.7 


4.7 


5.3 


4.2 


MSG 


14 


11 


1 


10 


850 MB TEMPERATURE CC) 










MEAN 


2.7 


2.4 


3.9 


%,* 


SD 


10.4 


10.2 


8.6 


10.5 


MSG 


15 


23 


13 


28 


850 MB WIND SPEED (kts) 










MEAN 


26 


24 


28 


21 


SD 


13 


12 


15 


11 


MSG 


24 


31 


15 


31 


CLOUD TOP (m) 










MEAN 


3553 


3520 


3666 


3603 


SD 


1982 


2015 


2052 


1847 


MSG 


43 


14 


ts 


55 


CLOUD THICKNESS (in) 










MEAN 


2912 


2883 


3170 


2842 


SD 


1993 


2055 


2049 


1791 


MSG 


61 


71 


46 


75 



P'3'- 



ii Gl £E ^ 

PilCIPITABLE WATER (tmm) 



MEAN 19.5 


18,8 


20.6 


18.7 


SD 12.0 


11-0 


10.8 


9,8 


MSG 17 


U 


15 


28 


TOTALS TOTALS 








MEAN 12.9 


12.7 


16.4 


16.3 


SD 17.6 


16.6 


15.7 


15.2 


MSG 15 


24 


13 


28 


MAXIMIJM RAINRATE (nm/15 rain) 








MEAN 2.0 


1.6 


l.j*P 


1.9 


SD 1 . 7 


1.2 


1.3 


2.4 


MSG 74 


77 


57 


67 


MAXIMUM RAINRATE (mm/ 30 tain) 








MEAN 2 . 7 


2.2 


t»i 


2.5 


SD 2,2 


1.7 


1.8 


3.0 


MSG 7% 


77 


57 


i7 


MAXIMUM RAINRATE (mni/hr) 








MEAN 3 . 5 


2.9 


3.5 


3.3 


SD 2.7 


2.1 


2.3 


3.9 


MSG f* 


76 


57 


f? 


DURATION OF PRECIPTTATION (hr) 








MEAN 7.5 




7.8 7.7 


9.6 






SD 5.0 


5.4 


5.2 


6.1 


MSG 77 


77 


67 


75 


PIECIPITATION RATE (mm/hr) 








MEAN 1 . 1 


0.9 


i,i 


0.9 


SD 0.9 


0.7 


0.8 


1.2 


MSG 77 


77 


67 


76 



€S 



SE 



NV 



DAILY AMOUNT (nun) 



MEAN 








6.3 


5.8 


6.1 


6.3 


m 








5.2 


4.3 


4.6 


5.8 


m§ 






. 


2 


1 





8 


MAXIMUM 6 HI 


AMOUNT 


(mm) 










MEAN 








4 . 3 


4.5 


4.6 


i*J 


m 








3.8 


4.0 


4.0 


5.0 


mm. 








17 


17 


13 


29 


% OP STNS WITH 


PRECIP 










MEAN 








0.66 


0.72 


0.72 


0.59 


il 








0.27 


0.23 


0.26 


0.31 


S04 (mg/1) 
















MEAN 








3.75 


2.99 


3.55 


1.32 


SD 








2.82 


3,08 


2.59 


1.32 


wm 








f 


1 


# 


* 


N03 (ng/1) 
















MEAN 








0.76 


0-61 


0.75 


0.26 


SD 








0.63 


0.52 


0.58 


0.22 


:^i 








7 


3 


6 


3 


APIOS PRECIPITATION 


(mm) 










MEAN 








7.0 


6.0 


7.5 


7.1 


SD 








6.9 


5.7 


6.9 


6.2 


MSG 








1 


t 


9 


9 



Appendix D 



Table Dl: Sulfate and Nitrate data as a functioo of the categorlea of the discrete 
variables. 



S04 



NO 3 







■"' 










MEAN 


SD 


MEAN 


SD 


SYNOPTIC TYPE 










WARM FRONT 


4.20 


2.73 


0.87 


0,66 


COLD FRONT 


3.66 


2.49 


0.67 





62 


STATIONARY FRONT 


4.96 


3,49 


0.87 





82 


OCCLUDED LOW 


3.24 


3.47 


0.76 





.53 


COLD LOW 


2.51 


1.49 


0.72 





57 


AIRMASS 


3.27 


2.84 


0.66 





.60 


WARM AND COLD FRONT 


4.46 


2,81 


0.89 





.69 


TRACK OF LOW 













ALllRTA LOW 
COLORADO LOW 
TEXAS LOW NORTH 
TEXAS LOW SOUTH 
GREAT LAKES 



3.56 


2.67 


0.82 


0.63 


4.23 


2.66 


0.74 


0.65 


2.55 


1.74 


0.59 


0.43 


3.69 


4.97 


0.71 


0.46 


2.96 


2.81 


0.84 


0.73 



AIRMASS 



CMA 
MA 
MP 
MT 

CA CMA 
CA MA 
CMA MA 
MA MP 
MP MT 



2.30 


2.21 


0.76 


0.56 


2.07 


1.16 


0.52 


0.21 


2-50 


1.48 


0.45 


0.29 


4.71 


2.21 


0.61 


0.35 


3.29 


0.05 


0.45 


0.12 


2.57 


1.76 


0,71 


0.63 


3.69 


3.03 


0.95 


0.79 


4.55 


2.31 


1.10 


0.71 


4.27 


2.72 


0.67 


0.60 


5.56 


3.72 


0.72 


0.53 



76 



p 



S04 



1103 



MEAN 



SD 



HIAN 



SD 



ISO 


MB TEAJECTOIY 










24 


HE SECTOR 












mu 


2.02 


1.94 


0.44 


0.48 




vm 


4.14 


2.95 


0.95 


0.71 




wsw 


4.13 


2.80 


0.84 


0,57 




ssw 


3.70 


3.00 


0.61 


0.66 




ESI 


2.04 


1,34 


0,33 


0.22 




ENE 


3.26 


0.11 


0.72 


0.36 




NHI 


1.22 


-- 


0.60 


~ - 


24 


HE lANGE (m) 


4.02 


3.69 


1.14 






<200 


0,87 




200-400 


5.11 


3.47 


0.90 


0.70 




400-600 


3.17 


2.00 


0.72 


0.56 




600-800 


4.09 


3.15 


0.84 


0.73 




800-1000 


3.04 


2.36 


0.76 


0.60 




1000-1200 


2,63 


1.97 


0.63 


0.55 




1200-1400 


3.01 


1.36 


0.54 


0.38 




1400-1600 


3.88 


3.42 


0.44 


0,20 




1600-1800 


3.37 


2.12 


0.41 


0.30 


48 


HE SECTOR 












NHW 


3.47 


4.22 


0,49 


0.43 




HHW 


3.62 


2,84 


0.88 


0.67 




WSW 


4.51 


2.72 


0,94 


0.71 




SSW 


3.68 


2.56 


0.68 


0.44 




SSE 


2.04 


1.34 


0.33 


0.22 




IHl 


3.19 


-• 


0.47 


im^m 




NNE 


0.85 


-- 


0.17 


mm 


48 


HR RANGE (KM) 












<200 


8.01 


8.22 


0.91 


0.42 




200-400 


4.74 


3.17 


0.93 


0.79 




400-600 


4.63 


4.28 


0.90 


1.02 




600-800 


4.75 


2,53 


0.97 


0.63 




800-1000 


3,82 


2.48 


0.70 


0.34 




1000-1200 


3.95 


2.61 


0.17 


0.53 




1200-1400 


2.73 


1.98 


0.55 


0.47 




1400-1600 


3.44 


2.77 


0.75 


0.75 




1600-1800 


3.54 


2.41 


0.83 


0,72 




1800-2000 


3.41 


2.44 


0.81 


0.67 




2000-2200 


2.59 


1.55 


0.78 


0.43 



If 



S04 



N03 



MEAi 



SD 



MEAN 



SD 



SUEFACE TRAJECTORY 










24 HE SECTOR 










NliW 


2.98 


3,21 


0.57 


0.42 


WNy 


3.75 


2.45 


0.76 


0.57 


WSW 


3.50 


2.29 


0.78 


0.67 


SSW 


4.47 


2.58 


0.78 


0.57 


SSE 


4.02 


2.81 


1.01 


0.74 


ISE 


5.41 


5.00 


0.99 


0.86 


ENE 


2.90 


1.93 


0.89 


0.78 


MME 


2.53 


2.58 


0.56 


0.65 


24 HE RANGE (KM) 










<200 


4.87 


3.S7 


0.82 


0.50 


200-400 


3.44 


2.29 


0.77 


0.54 


400-600 


4.47 


3.74 


0.98 


0.91 


600-800 


3.44 


2.01 


0,62 


0.34 


800-1000 


3.44 


2.40 


0.78 


0.77 


1000-1200 


2.12 


1.31 


0.40 


0.40 


1200-1400 


2.89 


1.70 


0.40 


0.21 


SURFACE TRAJECTORY 










48 HR SECTOR 










NNW 


2.75 


2.04 


0.55 


0.33 


WNW 


2.91 


2.67 


0.78 


0.59 


WSW 


4.13 


2.94 


0.94 


0.76 


SSW 


5.32 


2.57 


0.89 


0,49 


SSE 


4.02 


1.89 


0.90 


0.76 


ESE 


3.47 


1.94 


0.65 


0,42 


ENE 


4.03 


4,65 


0.90 


0.70 


NNE 


3.76 


4.36 


0.92 


1.17 


48 HR RANGE (KM) 










<200 


6.83 


3.43 


1.38 


0,60 


200-400 


3.87 


2.25 


0.74 


0.46 


400-600 


4.62 


3.25 


0.97 


0.92 


600-800 


4.50 


3.56 


0.84 


0.57 


800-1000 


3.03 


2.43 


0.81 


0.69 


1000-1200 


2.42 


1.74 


0.66 


0.51 


1200-1400 


2.17 


1.55 


0.49 


0.18 


1400-1600 


2.60 


2.03 


0.47 


0.38 


1600-1800 


1.70 


0.84 


0.22 


0.05 


1800-2000 


0.30 


0.11 


0.01 


0.00 


2000-2200 


1,05 


-• 


1.17 


-« 



S04 103 

HEAH SO MEAN SB 



TTPE OF PRECIPITATION 



CONTINUOUS 


2.94 


2.19 


0.74 


0.58 


CONVECTIVE 


4,31 


3.26 


0.75 


0.67 
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Table El: The frequency distributloii of the categorical variables for the top 307, 

of the SO, and HO. readings. 
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Table Fl. Contingeiicy table between preGlpitation type (TYl) and syfface 
24 hour sector (TR5) for warm months. 
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Table F2. Contingency table between synoptic type (ISN) and surface 48 hour range (TR8) for warm months. 
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Table F3. Contingency table between cloud type (TY3)and 850 mb 24 hour sector (TRl) for cold months. 



Figure Fl. SO^ versus preclpl table water for warm iionths. 
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Figure F2. SO^ versus duration of precipitation for warm months. 
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Figure F3. NO^ versus duration of precipitation for warm months. 
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Figure F4. NO3 versus amount of precipitation for warm months. 
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Figure F5. SO. versus cloud base height for cold months 
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Figure F6. SO^ versus amount of precipitation for cold monehs. 
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Flgtire F7: NO versus amoimt of for cold months. 
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8. Introduction to Addendum 



The 1983 meteorological data for the four regions were analysed in the same 

way as the 1982 data. However the following variables were not analysed as 

they showed no promise in the previous analyses of being significantly 

correlated with the SO, or NO- measurements: 

4 3 

track of the low 

cloud top height 

cloud thickness 

rainfall rate (mm/15 mln) 

rainfall rate (mm/30 min) 

rainfall rate (mm/hr) 

maximum 6 hr amount 

The abstracted data were then transferred to the MOE computer. The format was 
the same as that of the previous year except that the order of the 850 mb 
trajectories and the surface trajectories was interchanged. 

The purpose of this analysis was to check on the lack of positive results 
found from the 1982 data as contained in the main report entitled Final 
Report on the Determination of APIOS Precipitation Chemistry from the 
Climatology of Regional Precipitation Episodes, July 1986 (referred to as 
Kell). If the problem was due to an insufficient amount of data, then the 
results should improve with an additional year's analysis. Otherwise, the 
basic problem would be in the quality of the data for the problem at hand. 

In this report, the 1983 data will firstly be examined in relation to the 
1982 data to ensure that it can serve as a useful supplement to the data 
base. Then the relationships between the SO, and the N0„ readings and the 
meteorological variables for 1983, and 1982-83 will be carried out. Finally, 
the statistical model developed in Kell will be tested for the 1982-83 data 
set and subsequent results from experiments in forecasting annual deposition 
from a selected sample of days for both years will be presented. 



9. Acidic Precipitation Data for 1983 

Table 9,1 gives the SO. and NO^ averages, by region, for 1983 

Table 9.1: SO, and N0« averages by region* 

sw CI m WW 

NO OF OBS 137 153 132 131 

MEAN SO (mg/1) 3.30 2.88 3.05 1.24 
MEAN NO^ (mg/1) 0.61 0.59 0.57 0.22 
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A comparison of Table 9.1 with the SO, and NO valuea In Table CI of Kell 
reveals that there were 107. to 157, fewer observations In 1983. This 



was 



particularly noticeable in the southwest where the nomber of observations 
dropped from 183 to 137. The means of the SO measureiDents were somewhat 
lower in 1983 but relative positions of the values In the 4 regions remained 
the same. Duncan's Multiple Range Test (procedure ANOVA in SAS User's Guide: 
Statistics 1982 Edition) indicated that the northwest readings were 
significantly different from the other 3 but the remaining 3 regions were not 
significantly different from each other. The means of the NO. measurements 
were almost exactly the same In 1983 as in 1982 with only the southwest 
region being slightly lower. Duncan's Multiple Range Test showed that the 
northwest readings were significantly different from the other 3 regions, 
which in turn were not statistically significant from each other. But 
statistical tests of the averages of the 2 years showed that there is no 
statistically significant difference in the readings in the 4 regions between 
the 1982 and 1983. 



Table 9.2 gives the monthly variation in the 
southwest region for 1983, 



SO, and 
4 



NO 



-, readings in the 



Table 9.2: Monthly averages in the sonthfcrest region for 1983 



MONTI 



SO, 



m. 



Cmg/l) 



JANUARY 


3.96 


1.12 


FEBRUARY 


3.66 


0.84 


MARCH 


3.89 


0.89 


APRIL 


3.48 


0.66 


MAY 


3.38 


0.46 


JUNE 


3.72 


. 50 


JULY 


3.34 


0.55 


AUGUST 


4.41 


0.70 


SEPTEMBER 


3.30 


0.50 


OCTOBER 


2.80 


0.37 


NOVEMBER 


2.48 


0.45 


BEGEMBER 


2.14 


0.4S 



The average of the SO^ values for the warmer months (April to September) was 
3.58 mg/1 while for the colder months it was 3.01 mg/1. When these seasonal 
differences were tested for statistical significance (procedure ANOVA in SAS 
User's Guide: Statistics 1982 Edition), the F value was 3.22 and the 
significance probability was 0.08. In other words, the differences by season 
were not different at the 0.05 level, A further examination reveals that the 
very high values in May, June, and July in 1982 were not seen in 1983. 

The mean for the NO^ readings in the warmer months was 0.66 mg/1 and in the 
colder months 0.57 mg/1. An F value of 1.37 and significance probability of 
0.25 reveal that there was no significant difference in the seasonal readings 
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©f the NO., as was found also to be the case in 1982. With the exception that 
some of the very high readings In 1982 were not seen in 1983, the SO^ and NO^ 
measurements In 1983 appear to be comparable with the 1982 readings. 



10. Meteorological Data for 1983 



Appendix G of this addendum summarizes the 1983 data for the 4 regions. For 
the southwest region the analysis of variance was performed on the 
meteorological data to see if there was any difference in the data between 
1982 and 1983, Table 10.1 summarizes the variables which showed a 
statistically significant difference in their values between the 2 years at 
the 107, significince level. 

Table 10.1: Variables which had statistically significant differences (at the 
107, level) between 1982 and 1983 for the southwest region. 



VARIABLE P>P 

surface temperature 0.01 

surface dewpoint 0.01 

precipitable water 0.05 

total totals Index 0.04 

X of stns with rain 0.0001 

a 1 mass 0.09 

nature of pcpn 0.002 

surface wind dir 0.02 

850 mb wind dir 0.04 



REASON 

1983 was 3.0- C warmer 
1983 was 2.9''C warmer 
1983 was 3.0 mm higher 
1983 was 4.6 higher 
1983 was 107, higher 

more mT m' air and less cA and cmA air in 

1983 

less snow and more rain In 1983 

more easterly winds in 1983 

more easterly winds in 1983 



From Table 10.1, the differences between the 2 years meteorologically when 
precipitation was recorded at the APIOS stations could be characterized in 
the following terms. In 1983, conditions were warmer and more moist with the 
rain being more widespread than in 1982. The warmer alrmasses (mT, mP) were 
more in evidence in 1983. Also, there were more easterly winds during the 
precipitation in 1983 than 1982. 



11. A Comparison between SO,, 



NO^ and the Meteorological Variables for 1983 
3 



For the discrete meteorological variables. Table Hi in Appendix H summarizes 
the SO, and N0„ data as a function of the discrete categories while Table II 
in Appendix I gives the frequency distribution of the categorical variables 
for the top 307, of the SO, and NO^ readings. When compared with the results 
of 1982 (Table 4.3 c in Kell) some of the findings were not substantiated In 
1983, In particular, the characterization of high SO, episodes with 
stationary fronts and stratus clouds was not seen in 1983. For the NO^ 
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readiags, the high values associated with snowfall, the cA, cmA, and 
airmasses, and stratys clouds In 1982 were not seen in the 1983 data. 



mA 



Table H2 gives the correlation coefficients and significance probabilities 
between the SO. and N0» values and the continuous meteorological variables 
while Table 12 gives the averages and standard deviations of the continuous 
meteorological variables for the top 307. readings of the SO, and the ND^. 
Inspection of these tables reveals that cloud base height and the duration of 



the precipitation were not related to variations In SO, and iO* In the saine 

way as in 1982 when lower clouds base heights and lower durations were 

associated with high values of SO, and NO,. 

4 3 

Combining the 2 years data, high episodes of SO, and N0„ can be characterized 



by the meteorological conditions as given in Table 



h. 



1 . In order for a 



factor to be considered to be important, the mean of the SO, or N0„ readings 
for that category had to be the highest or second highest among the means for 
all the categories for that variable. Also, the frequency (In percentage) of 
that factor had to increase when considering the top 307. of the SO, or NO^ 
readings as compared to the frequency in the whole data set. 

Table 11. li Meteorological conditions associated with high episodes of SO, 
and N0„ as derived from an inspection of the tables in Appendix 
G, H, I of this addendum and their counterparts, nanely Appendix 
C, D, E, and Tables 4.6 and 4.7 of Kell, 



,SQ. 



warm fronts 



850 mb trajectories originating in 
the f SW 



850 mb trajectories originating in 
the WSW 



24 hr 850 mb trajeetory ranges less 24 hr 850 mb trajectory ranges less 
than 400 km than 400 km 

24 hr surface trajectories originating 24 hr surface trajectories 
in the SSW to ESE originating in the SSW to ESE 

convective situations 

rain 

surface and 850 mb winds from the 
WSW to SSE 

bT and mp air 



high surface temperatures, dewpoints, low surface temperatures, dewpoints 



precipitable water, totals index 



wind speeds , rainfall amounts, and 
area of pcpn 



low winds speeds, rainfall amounts, 
and area of pcpn 
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12. Statistical Tests Between SO' , NO^ and the Meteorological Variables 

4 3 

An analysis of variance was done between SO, and N0„ and the discrete 
meteorological variables for the 1983 data set (Table 12.1). The data were 
divided into warm and cold months, as before. The underlined values indicate 
those variables whose significance probability was 0,10 or less, A comparison 
between Table 12.1 with Table 4.4 in Kell reveals that the variables which 
were significant from one year to the next were not the same ones. For the 
SO, readings, the variables which were significant for both years were, in 
the warm months, the 24 hr surface and 850 mb trajectory sectors and in the 
cold months, the nature of the precipitation. For the N0„ measurements, only 
the 24 hr 850 mb trajectory range in the cold months was common to both 
analyses . 

Table 12.2 summarizes the results for 1983 of the regressions between the SO 
and NO^ measurements and the continuous variables. In this case the only 
variables significant at the 107. level for both analyses were the cloud base 
height and the amount of precipitation with SO. in the cold months, and the 
amount of precipitation with N0« . For the cloud base height, the slope 
estimate is also of a different sign for the 2 years. 

These comparisons point out very sharply that an Increase in the quantity of 
the data to 2 years is not going to help Improve the results that were 
obtained with 1 year's data. 



13- The Statistical Model 



The statistical model based on the analysis of covarlance for the 1982 data 
set was run using the combined 1982-83 data set for all 4 regions. In order 
to assess the possible improvement In the model's performance, the variables 
which were deemed best in the 1982 analysis were used. Tables 13.1 to 13.4 
contain the results of the ANCDVA for the four regions. The significance 
probabilities (P>F) were very similar in the 2 analyses. However the amount 
of variance explained in the analyses (R) using the 1982-83 data set was 
substantially less than in the analyses with only the 1982 data. In the 
former one, R averaged 0.49 while in the latter it averaged 0.66. In all 16 
cases, SO, and N0« for the warm and cold months for the 4 regions, R was less 
using the 1982-83 data set than using the 1982 data set. In 13 out of the 16 
cases the mean square error was larger with the 1982-83 data set. This 
clearly suggests that having more data does not overcome the difficulties 
encountered in the first year analysis and that the relationship between the 



SO, and N0„ measurements and the meteorological variables presented is not 
strong enough to be used effectively in a statistical model. 
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so, and NO- and the discrete 
are the F statistic and the 



Table 12.1. An analysis of variance between 
variables for 1983. Given 
significance probability (P>F). Underlined values are those 
deemed significant enough to warrant further statistical 
testing. 



SO, 



WARM MONTES 



COLD MONTHS 



P>F 



P>F 



synoptic type 



2.00 



0.09 



alrmass 


1 


.57 


0.17 


850 mb 24 hr sector 


4 


.54 


0.01 


range 





.33 


0.89 


48 hr sector 


2 


91 


0.03 


range 





.62 


0.75 


sfc 24 hr sector 


3 


. 11 


0.01 


range 


1 


.29 


0.28 


48 hr sector 


1 


64 


0.15 


range 


1 


.34 


0.26 


type of pcpn 


1 


49 


0.23 


nature of pcpn 


3 


33 


0.07 


cloud type 


1 


91 


0.09 


sfc wind dir 





.73 


0.66 


850 mb wind dir 


2 


.36 


0.05 



,04 
.71 
.53 
.03 
,07 
.03 
.42 
.85 
.33 
.33 
.06 
.82 
.83 
.53 
.29 






40 





,12 





.18 





.41 





41 





43 





.88 





.56 





93 





.94 





08 





004 





.55 





.17 



0.27 



WARM MONTHS 



COLD MONTHS 



W>f 



P>F 



synoptic type 


1 . 70 


0.15 





.71 


0.62 


alrmass 


2.06 


. 07 





.66 


0.71 


850 mb 24 hr sector 


2.10 


0.11 





.41 


0.89 


range 


0.59 


0,71 


2 


.23 


0.07 


48 hr sector 


1.10 


0.37 


1 


.24 


0.32 


range 


0,58 


0,79 





.39 


0.88 


sfc 24 hr sector 


1.94 


0,08 





.32 


0.94 


range 


0.36 


0.87 





.53 


0.81 


48 hr sector 


0,75 


0,63 





.46 


0.85 


range 


1.14 


. 36 





.39 


0.90 


type of pcpn 


D . 78 


0.38 


2 


.58 


0,11 


nature of pcpn 


2 . 54 


0.12 





.22 


0.88 


cloud type 


. 96 


, 47 





.73 


. 63 


sfc wind dir 


0.26 


0.98 





.29 


0.96 


850 mb wind dir 


1.46 


0.22 





.71 


0.67 
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Table 12.2 i 



Slope estimate aiid slgnif Icance probability between SO, , N0«, and the 
continuous meteorological variables for 1983. 



SO, 



warm months 



cold months 



M. 



P>F 



bl 



CLOUD BASE HEIGHT 
SURFACE TEMPERATURE 
SURFACE DEWPOINT 
SURFACE WIND SPEED 
850 MB TEMPERATURE 
850 MB WIND SPEED 
PREGIPITABLE WATER 
TOTAL TOTALS 
DURATION OF PCPN 
AVERAGE RAINRATE mm/hr 
AMOUNT OF PCPN 
I STNS WITH PCPN 



0.0005 


0.18 


0.0015 


0.07 


0.047 


0.21 


0.054 


0.32 


0.026 


0.50 


. 04 


0.43 


-0.006 


0.91 


-0.14 


0.03 


0.05 


0.14 


0.083 


0.15 


0.009 


0.66 


-0.004 


0.86 


0,013 


0.49 


0.035 


0,42 


0.012 


0.39 


0.019 


0.36 


-0.07 


0.25 


0.035 


0.01 


-0.43 


0.11 


0.034 


0.01 


-0.05 


0.17 


-0.14 


0,03 


-2.48 


0.07 


-0.71 


0.61 



WBs 



warm months 



cold months 



■bl 



P>F 



bl 



P>F 



CLOUD BASE HEIGHT 
SURFACE TEMPERATURE 
SURFACE DEWPOINT 
SURFACE WIND SPEED 
850 MB TEMPERATURE 
850 MB WIND SPEED 
PRECIPITABLE WATER 
TOTAL TOTALS 
DURATION OF PCPN 
AVERAGE RAINRATE mm/hr 
AMOUNT OF PCPN 
% STNS WITH PCPN 



0.0001 


0.19 


0,0001 


0.60 


0.001 


0.85 


-0.002 


0.84 


0.001 


0.85 


-0.002 


0.88 


0,003 


0.79 


-0.03 


0.01 


0.003 


0.63 


0.008 


0.4f 


0.004 


0.36 


-0.003 


0.55 


0.0009 


0,80 


0.001 


0.89 


0.003 


0.22 


0.003 


0,38 


0.01 


0.36 


0.009 


0.0001 


0.95 


0.07 


0.008 


0,0001 


0.02 


0.02 


-0.02 


0.06 


0.67 


0.01 


-0.12 


0.65 
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Table 13.1: lesults ©f the ANCQVA for the Southwest Region for 1982-83 



fARIABLE 



WARM MONTHS 

i P>F F 



p>F' 



type of pcpn 
sfc 24 hr sector 
Interaction 
preclpl table water 
duration of pcpn 



synoptic type 
sfc 48 hr range 
Interaction 
duration of pcpn 
anount of pcpn 



4 

6.17 0.02 

4,97 0.0001 

2.15 0.05 

0.97 0.32 

0.56 0.46 

N0„ 
3 

2.47 0.03 
1.31 0.26 

1.48 0.11 
2.10 0,15 
1.42 0.24 



3.38 0,0001 0.42 4.44 



1.62 0.04 0.47 0.14 



VARIABLE 



COLD MOHTHS 

M P>F 



P>F 



cloud type 

850 24 hr sector 

Interaction 
cloud base height 
amount of pcpn 



cloud type 

850 24 hr sector 

interaction 

cloud base height 

amount of pcpn 



SO, 
4 

3.13 0.005 

2.64 0.014 

5.95 0.0001 

0,90 0.34 

24.67 0.0001 

10' 



4,87 0.0001 0.53 2,28 



3.03 


0,006 


5.65 


4.58 


0.0002 




6.54 


0.0001 




0.25 


0.62 




27.1 


0.0001 





0.0001 0.57 0.17 



Table 13*2: Results of the ANCOVA for the Central Region for 1982-83. 



VARIABLE 



WARM MONTHS 
F P>F F 



P>F 



SO. 



type of precipitation 

sfc 48 hr sector 

Interaction 

duration of precipitation 

sfc temperature 



sfc 48 hr sector 
850 wind direction 
interaction 
duration of pcpn 
sfc wind speed 



7.58 


0.007 


4.98 


0.0001 


0.30 


0.937 


3.35 


0.070 


0.41 


0.526 



NO, 



VARIABLE 



4.71 


0.001 


3.68 


0.002 


1.52 


0.107 


0.37 


0.544 


0.12 


0.729 


COLD 


MONTHS 


F 


F>F 



3.00 0.0004 0.32 5.04 



2.41 0.001 0.51 0.08 



P>F R 



E 



airmass 

sfc 48 hr sector 

Interaction 

sfc temperature 



SO, 



2.41 


0.026 


1.82 


0.104 


1.92 


0.008 


0.37 


0.545 



1.95 0.0033 0,48 4.36 



sfc 48 hr sector 
nature of pcpn 
Interaction 
cloud base height 
pcpn amount 



NO, 



3.39 


0.003 


1.91 


0.115 


5.14 


0.0001 


0.22 


0.640 


14.54 


0.0002 



4.24 0.0001 0.49 0.33 
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Table 13.3: The ANCOVA for Southeast region for 1982-83. 

WARM MONTHS 
VARIABLE i P>F F P>F R 



850 wind direction 
sfc 48 hr sector 
Interact ion 
duration of pcpn 
sfc wlmd speed 



synoptic type 
sfc 48 hr range 
Interaction 
duration of pcpn 
afflount of pcpn 



VAEIABLE 



1.02 
5.21 
3.31 
0.13 
3.16 



SO, 



3.70 


0.003 


4.38 


0.005 


1.55 


0.110 


0.78 


0.380 


1,88 


0.175 



2.59 0.001 0.56 4.24 



NO3 

0.428 3.11 

0.0002 

0.0001 

0.718 

0.080 



0.0001 0.65 0.07 



COLD MONTHS 

;i P>F 



P>F 



SO. 



nature of pcpn 
sfc 24 hr sector 
Interaction 
sfc temperature 



850 24 hr sec tor 
nature of pcpn 
interaction 
precipltable water 
amount of pcpn 



4.45 


0.0002 3.79 


0.0001 0.50 


3.46 


0.019 




3.70 


0.0001 




1.66 


0.201 





2.24 



NO, 



2.81 


0.014 


0.57 


0.639 


0.85 


0.575 


4.43 


0.038 


13.8 


0.0003 



tntZ 0.005 0.31 0.39 
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Table 13.4: Results of ANCOVA for northwest region for 1982-83. 

WARM MONTHS 



VARIABLE 



F 



P>F 



F 



P>F 



SO, 



sfc 48 hr range 
alrmass 
interaction 
cloud base height 
sfc temperature 



alrmass 

sfc 48 hr sector 

interaction 

stns with pcpn 

sfc temperature 



4.15 


0,0001 2.54 0.0004 0.54 0.76 


5.12 


0,0001 


0.94 


0,538 


0.00 


0.988 


1.84 


0.179 



NO, 



VARIABLE 



5.07 0.0001 


2.35 0.0003 0.54 0.02 


2.10 0.052 




1.59 0.052 




1.72 0.193 




1,31 0.256 




COLD MONTHS 




f P>F 


F P>F R E^ 



sfc 24 hr range 
synoptic type 
interaction 
sfc temperature 
precipl table water 



3.55 
2.13 
0.77 
0.10 
4.84 



SO, 

0.003 1.77 
. 065 
0.714 
0.749 
0,032 



0,033 0.561 0.69 



alrmass 

850 48 hr sector 

interaction 

pcpn amount 



NO, 



2.61 


0.027 1.60 


0.075 0.44 


1.42 


0.226 




0.99 


0.469 




3.97 


0.052 





0.03 
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14.0 Annual Deposition Prediction Results 

islng the statistical model developed in the Section 13, predictions of 
annual deposition using the weighted average approach for regional averages 
of SO, and NO^ were made for both years. The days used for the selected 
sample In the 1982 analysis were exactly the same as those used in Kell, For 
1983 the following days were used for the selected sample: 





NO. of DAYS 


January 10-14,29-31 


8 


February 1-6 


6 


March 18-21 


4 


April 6-17,28-30 


15 


May 1-4 


4 


June 3-6 


4 


July 28-31 


4 


August 8-11 


4 


September 15-23 


9 


October 11-14 


4 


November 27-30 


4 


December 3-7 


5 



Table 14.1 and 14.2 |ive the results of the deposition prediction for 1982 
and 1983 respectively. As in Kell, two predictions were made for annual 
deposition. Prediction One (Pi) was based on the categorizations of the 
statistical model while Prediction Two (P2) was made assuming no 
categorizations. Table 14.3 gives the absolute fractional errors associated 
with the deposition statistics. Overall, the errors were 10 to 15X in the 
predictions with approximately one-third of the days being sampled. In 
keeping with the results of Kell, PI did not perform noticeably better than 
P2 for the data as a whole or any particular subset of the data. The N0» 
predictions for the northwest region showed PI to be better than P2, but this 
is attributed to statistical fluctuations. The one exception appears to be 
whem the number of days In the selected sample has less than one-third of the 
total number of days. In this case PI was better 11 times and P2 was better 7 
times. When the fraction of days sampled was greater than one-third, then f2 
performed much better than Pi, having the better prediction 11 out of 14 
times . 



15, Summary and Conclusions 



From the data presented In this addendum, the following conclusions can be 
draiin. 

JL. With the exception that some of the very high readingi In 1982 were 
not seen in 1983, the SO. and N0» measurements in 1983 appear to be 
comparable with the 1982 readings, 

2. For the precipitation collected in 1983, conditions were warmer and 
more moist with the rain being more widespread than in 1982. The 
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Table 14.1: Deposition results using weighting averaging for 1982. 



REGION 



MONTHS 



OBSERVED PREDICTION PREDICTION NO. of DAYS 
DEPOSITION ONE TWO SAMPLE SET 



southwest 



central 



southeast 



northwest 



SO, warn 

NO- warm 

SO, cold 

NO cold 

SO. warm 
N0_ warm 
SO, cold 
N0_ cold 

SO, warm 
NO- warm 
SO, cold 
NO- cold 

SO, warm 

««4 

NO- warm 

SO, cold 
N0„ cold 



1699.530 


1922.500 


1754.060 


27 


50 


271.020 


312.510 


285.100 


29 


60 


1087.860 


1174.820 


1015.400 


32 


85 


236.777 


271.394 


264.987 


31 


84 


1391.590 


1304.690 


1104.570 


23 


66 


154.667 


169.881 


131.031 


21 


61 


702.783 


835.454 


688.642 


25 


77 


183,780 


174.262 


196.468 


17 


64 


1004.670 


905.640 


987.834 


17 


41 


139.341 


159.480 


143.585 


18 


46 


939.055 


1084.560 


877.920 


16 


62 


155.186 


152.910 


147.631 


12 


49 


646.690 


543.240 


539.745 


15 


58 


105.258 


98.300 


100.621 


16 


69 


245.207 


148.838 


244.141 


10 


47 


53.576 


49.730 


72.145 


12 


51 



Table 14.2: Deposition results using weighting averaging for 1983. 



REGION 



MONTHS 



OBSERVED PREDICTION PREDICTION NO. of DAYS 
DEPOSITION ONE TWO SAMPLE SET 



southwest 



central 



southeast 



northwest 



SO, warm 
4 

NO warm 
SO, cold 
NO cold 

SO, warm 

N0_ warm 

SO, cold 

4 
N0_ cold 

SO, warm 

««4 

NO- warm 

SO, cold 

NO- cold 

SO, warm 
NO^ warm 
SO, cold 
N0» cold 



1548.840 


1520.490 


1433.490 


27 


48 


165.424 


155.737 


134.230 


23 


45 


730.594 


629.815 


695.282 


22 


57 


133.634 


110.604 


111.290 


22 


57 


1136.850 


889.920 


914.942 


19 


51 


141.235 


96.441 


101.769 


17 


45 


626.560 


604.796 


507.079 


16 


64 


145.667 


163.950 


148.751 


13 


58 


1495.060 


1930.680 


1681.020 


23 


53 


209.851 


239.801 


248.445 


23 


57 


1010.660 


845.205 


1032.380 


18 


69 


217.789 


292.370 


244.205 


22 


66 


627.375 


695.720 


511.796 


15 


54 


105.629 


92.270 


87.985 


16 


67 


206.600 


162.240 


155.729 


6 


36 


25.718 


23.230 


18.889 


3 


27 
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Table 14. 3i Deposition errors using the statistical model based on the 
combined 1982-83 data set 

ABSOLUTE FRACTIONAL ERIOI 



SO. 



NO. 



DATA SET 


PREDICTION 


PREDICTION 


BIST 


PREDICTION 




ONE 


TWO 


ONE 


TWO 


ONE 


TWO 


whole 


0.15 


0.10 


0.13 


0.14 


14 


18 


southwest 


0.13 


0.15 


0.15 


0.13 


*: 


*. 


central 


0.09 


0.06 


0.13 


0.13 


■i 


6 


southeast 


0.18 


0.06 


0.16 


0.10 


1 


i 


northwest 


0.22 


0.15 


0.09 


0.21 


* 


i 


warm months 


0.14 


0.12 


0.14 


0.14 


8 


8 


cold months 


0.17 


0.08 


0.13 


0,14 


6 


10 


1982 


0.16 


0.07 


0.09 


0.11 


i 


10 


1983 


0.15 


0.14 


0.17 


0.17 


s 


3 


<0.35 sampled 


0.16 


0.12 


0.11 


0,14 


11 


7 


>0,35 sampled 


0.14 


0,08 


0.16 


0.15 


1 


It 
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warmer airmasaes (mT, mP) were more In evidence in 1983. Also 
there were more easterly winds during the precipitation in 1983 
than 1982. 



1« High episodes of SO- and NO^ concentrations were associated with 24 
hr 850 mb trajectories originating in the WSW at ranges less than 
400 km and 24 hr surface trajectories originating In the SSW to 
ESE. Also, high SO, readings were associated with warm fronts, mT 
and mP airmasses , convective rain situations, surface and 850 rab 
winds from the WSW to SSE, high surface temperatures, dewpoints, 
preclpitable water, and total totals index, and low wind speeds, 
rainfall amounts and area of precipitation. High NO. readings were 
also associated with low surface temperatures, dewpoints ,wlnd 
speeds, rainfall amounts, and area of precipitation. 

km The statistical model performed more poorly with the 2 years data. 
The amount of variance explained by the model dropped from about 
two- thirds to one-half with the addition of the second year of data 
indicating that the extra data added more noise to the set since 
the relationship among the SO, and the NO^ measurements and the 
meteorological variables was not sufficiently well established with 
the data at hand. 

5. The experiments to predict annual deposition showed that overall 
there was no advantage to using the categorizations of the days 
based on the statistical model and that the errors were of the 
order of 10 to 157.. However, predictions of annual deposition using 
the use of categorizations of the days when the selected sample 
contained less than one-third of the total number of days was 
marginally better than without using the categorizations, but not a 
statistically significant level. 
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Appendix 6 



The foUowlnig is a summary of the data avallablt for the 4 regions (Southwest 
(S¥), Central (CE), Soytheast (SE) , and Northwest (NW)). Sulfate, nitrate, and 
APIOS rainfall data are averages from the APIOS stations which reported 
precipitation on a particular day. For the discrete variables the first number Is 
the number of observations of a given category followed by its percentage of the 
available observations in brackets. For the continuous variables, the mean and 
standard deviation (SD) is given. Missing data (Indicated by MSG) refers to either 
the fact that data were not available, a reliable determination of the variable 
could not be made with the data available (e.g. predominant cloud type), or there 
were no data (e.g. no rainfall reports). 



Table Git A summary of the meteorological data collected. 



TOTAL NO OF OBS 



137 



CE 
153 



SE 
132 



NW 
131 



SYNOPTIC TYPE 



WARM FRONT 36 (26) 29 

COLD FRONT 24 (18) 34 

STATIONARY FRONT 9 (7) 14 

OCCLUDED LOW 43 (31) 45 

COLD LOW 3(2) 7 

AIRMASS 18 (13) 21 

WARM AND COLD FNTS 4(3) 3 



(19) 

(22) 

(9) 

(29) 

(5) 

(14) 

(2) 



26 
22 
14 
41 

8 
15 

6 



(20) 

(17) 

(11) 

(31) 

(6) 

(11) 

(5) 



21 

21 

11 

66 

2 

9 

1 



(16) 

(16) 

(8) 

(50) 

(2) 

(7) 

(1) 



AIRMASS 



CA 

CMA 

MA 

MP 

MT 

CA CMA 

CA MA 

CMA MA 

MA MP 

MP MT 

MS,C 



13 (9) 


19 (12) 


16 (12) 


21 (16) 


5 (4) 


9 (6) 


4 (3) 


6 (5) 


16 (12) 


20 (13) 


15 (11) 


20 (15) 


5 (4) 


5 (3) 


3 (2) 


9 (7) 


2 (1) 





3 (2) 


2 (2) 


2 (1) 


7 (5) 


2 (2) 


5 (4) 


10 (7) 


13 (9) 


11 (8) 


21 (16) 


19 (14) 


19 (12) 


20 (15) 


6 (5) 


46 (34) 


42 (27) 


43 (33) 


21 (16) 


19 (13) 


19 (12) 


15 (U) 


19 (15) 








M 


;t 
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850 MB TRAJECTORY 







S£ 


CE 


SE 


^ 


24 HR SECTOR 










NN¥ 




5 (4) 


3 (2) 


5 (4) 


15 (12) 


UNU 




30 (25) 


39 (27) 


17. (14) 


24 (19) 


WSW 




49 (40) 


52 (36) 


45 (37) 


39 (30) 


SSV 




28 (23) 


32 (22) 


30 (25) 


30 (23) 


SSE 




4 (3) 


8 (6) 


13 (11) 


9 (7) 


ESE 




3 (2) 


5 (3) 


6 (5) 


2 (2) 


ENE 




1 (1) 


3 (2) 


3 (2) 


6 (5) 


NNE 




2 (2) 


3 (2) 


2 (2) 


4 (3) 


MSG 




15 


8 


11 


2 


24 HR RANGE (KM) 










<200 







3 (2) 


4 (3) 


11 (9) 


200-i 


:»oa 


17 (14) 


18 (12) 


10 (8) 


39 (30) 


400-* 


600 


39 (32) 


40 (28) 


34 (28) 


41 (32) 


600-1 


300 


31 (25) 


42 (29) 


37 (31) 


25 (19) 


800- 


1000 


16 (13) 


16 (11) 


19 (16) 


10 (8) 


1000 


-1200 


9 (7) 


18 (12) 


13 (11) 


2 (2) 


1200 


-1400 


10 (8) 


7 (5) 


3 (2) 


1 (1) 


1400 


-1600 





1 (1) 


1 (1) 





MSG 




15 


8 


11 


- 2 


48 HR SECTOR 










nm 




t (11) 


6 (5) 


4 (5) 


7 (8) 


WNW 




20 (24) 


39 (33) 


20 (23) 


24 (26) 


WSW 




42 (50) 


40 (34) 


42 (48) 


23 (24) 


ssw 




6 (7) 


14 (12) 


8 (9) 


24 (26) 


SSE 




4 (5) 


9 (8) 


9 (10) 


5 (5) 


ESE 




1 (1) 


4 (3) 


1 (1) 


3 (3) 


ENE 




1 (1) 


4 (3) 


2 (2) 


4 (4) 


MNE 




1 (1) 


1 (1) 


2 (2) 


4 (4) 


MSG 




53 


36 


44 


37 
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SW CE ,g W. 

48 HR RANGE (KM) 



<200 


2 


(2) 


2 


(2) 


1 


(1) 


5 (5) 


200-400 


6 


(7) 


4 


(3) 


7 


CS) 


9 (10) 


400-600 


7 


(8) 


9 


(8) 


6 


(7) 


12 (13) 


600-800 


15 


(18) 


18 


(15) 


10 


(11) 


29 (31) 


800-1000 


23 


(27) 


23 


(20) 


15 


(17) 


18 (19) 


1000-1200 


9 


(11) 


20 


(17) 


14 


(16) 


14 (15) 


1200-1400 


12 


(14) 


21 


(18) 


13 


(15) 


4 (4) 


1400-1600 


8 


(10) 


11 


(9) 


9 


(10) 


3 (3) 


1600-1800 


2 


(2) 


7 


(6) 


8 


(9) 





1800-2000 







2 


(2) 


4 


(4) 





2000-2200 


t 









1 


(1) 





MSG 


53 




36 




44. 




ai 


SURFACE TRAJECTORY 
















24 HR SECTOR 
















NH¥ 


7 


(6) 


7 


(5) 


8 


(6) 


14 (12) 


wm 


19 


(17) 


20 


(14) 


10 


(8) 


5 (4) 


WSW 


13 


(12) 


24 


(17) 


24 


(20) 


9 (8) 


SSW 


27 


(24) 


29 


(20) 


26 


(21) 


27 (22) 


SSE 


20 


(18) 


20 


(14) 


13 


(10) 


28 (23) 


ESE 


10 


(9) 


29 


(20) 


19 


(15) 


16 (13) 


ENE 


10 


(9) 


4 


(3) 


13 


(11) 


12 (10) 


NNE 


5 


(4) 


11 


(8) 


10 


(8) 


9 (8) 


MSG 


26 




9 




9 




11 


24 HR RANGE (KM) 
















<200 


11 


(10) 


9 


(6) 


8 


(6) 


16 (13) 


200-400 


27 


(24) 


30 


(21) 


26 


(21) 


39 (32) 


400-600 


26 


(23) 


28 


(19) 


24 


(19) 


24 (20) 


600-800 


25 


(23) 


42 


(29) 


32 


(26) 


20 (17) 


800-1000 


12 


(11) 


20 


(14) 


21 


(17) 


16 (13) 


1000-1200 


7 


(6) 


12 


(8) 


12 


(10) 


4 (3) 


1200-1400 


2 


(2) 


3 


(2) 










1400-1600 


1 


(1) 















MSG 


26 




9 




f 




12 
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SURFACE TRAJECTORY 





M 




CI 




SE 




w 


48 HR SECTOR 
















liN¥ 


9 


(10) 


18 


(15) 


8 


(9) 


10 (10) 


WNW 


13 


(15) 


19 


(16) 


10 


(11) 


4 (4) 


WSM 


12 


(14) 


18 


(15) 


17 


(19) 


3 (3) 


SSW 


21 


(24) 


16 


(13) 


15 


(17) 


25 (24) 


SSE 


8 


(9) 


14 


(12) 


4 


(4) 


22 (21) 


ESE 


9 


(10) 


16 


(13) 


10 


(11) 


16 (15) 


ENE 


5 


(6) 


6 


(5) 


11 


(12) 


18 (17) 


NNE 


11 


(12) 


13 


(11) 


14 


(16) 


5 (5) 


MSG 


49 




33 




43 




28 


48 HR RANGE (KM) 
















<200 


5 


(6) 


9 


(8) 


3 


(3) 


8 (8) 


200-400 


11 


(12) 


13 


(11) 


12 


(13) 


16 (16) 


400-600 


21 


(24) 


18 


(15) 


16 


(18) 


22 (21) 


600-800 


23 


(26) 


27 


(22) 


17 


(19) 


21 (20) 


800-1000 


10 


(11) 


23 


(19) 


23 


(26) 


15 (15) 


1000-1200 


7 


(8) 


13 


(11) 


6 


(7) 


18 (17) 


1200-1400 


8 


(9) 


10 


(8) 


9 


(10) 


2 (2) 


1400-1600 







7 


(6) 


3 


(3) 


1 (1) 


1600-1800 


3 


(3) 















MSG 


4f 




33 




43 




11 


TYPI OF PRECIPITATION 


62 


(46) 


62 


(41) 


65 


(51) 




CONTINUOUS 


60 (50) 


CONVECTIVE 


72 


(54) 


90 


(59) 


63 


(49) 


60 (50) 


MSG 


1 




1 




♦ 




tt 


NATURE OF PRECIPITATION 
















RAIN 


97 


(72) 


95 


(62) 


93 


(73) 


82 (68) 


SNOW 


20 


(15) 


43 


(28) 


12 


(9) 


21 (18) 


FREEZING 


5 


(4) 


1 


(1) 


1 


CD 


5 (4) 


MIXED 


12 


(9) 


13 


(9) 


22 


(17) 


12 (10) 


MSG 


3 




1 




4 




11 



lis 



sv 



ei 



SE 



PREDOMINAMT CLOUDS 



TOWEIING CUMULUS 


7 


(5) 


9 


(6) 


2 


(2) 


11 (11) 


ALTOCUMULUS 


4 


(3) 


3 


(2) 


3 


(2) 


2 (2) 


ALTOSTRATUS 


1 


(1) 


2 


(1) 










CUMULONIMBUS 


22 


(17) 


20 


(13) 


21 


(17) 


17 (18) 


CUMULUS 


10 


(8) 


10 


(7) 


2 


(2) 


3 (3) 


NIMBOSTRATUS 


21 


(16) 


10 


(7) 


11 


(9) 


1 (1) 


STRATOCUMULUS 


48 


(36) 


63 


(42) 


75 


(59) 


49 (51) 


STRATUS 


19 


(14) 


32 


(21) 


13 


(10) 


13 (14) 


MSG 


5 




4 




1 




35 


SURFACE WIND SECTOR 
















NN¥ 


9 


(7) 


13 


(8) 


11 


(8) 


12 (9) 


WNW 


16 


(12) 


22 


(14) 


11 


(8) 


15 (12) 


WSW 


24 


(18) 


21 


(14) 


30 


(23) 


15 (12) 


SSW 


22 


(16) 


22 


(14) 


21 


(16) 


12 (9) 


SSE 


12 


(9) 


29 


(19) 


13 


(10) 


18 (14) 


ESE 


33 


(24) 


26 


(17) 


11 


(8) 


15 (12) 


ENE 


16 


(12) 


10 


(7) 


18 


(14) 


31 (24) 


NNE 


3 


(2) 


10 


(7) 


16 


(12) 


9 (7) 


CALM 


1 


(1) 







1 


(1) 


2 (2) 


MSG 


1 














2 


850 MB WIND SECTOR 
















NNW 


8 


(7) 


13 


(10) 


10 


(9) 


17 (15) 


WNW 


23 


(19) 


27 


(21) 


19 


(16) 


30 (27) 


WSW 


32 


(27) 


48 


(37) 


51 


(44) 


26 (23) 


SSW 


23 


(24) 


28 


(21) 


25 


(21) 


14 (12) 


SSE 


10 


(8) 


9 


(7) 


7 


(6) 


9 (8) 


ESE 


11 


(9) 


1 


(1) 


1 


(1) 


7 (6) 


ENE 


2 


(2) 


2 


(2) 


1 


(1) 


4 (4) 


NNE 


4 


(3) 


3 


(2) 


3 


(3) 


6 (5) 


MSG 


19 




22 




15 




18 
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§w 



CLOUD BASE HEIGHT (M) 



MEAN 
SD 

nsG 



SmFACE TEMPERATUIE ("C) 



MEAlf 
SD 

MSG 



SURF ACE DEWPOINT ("C) 



MEAN 
SB 



SURFACE WIND SPEED (m/s) 



MEAI 
SD 



850 MB TEMPERATUIE (-C) 



MIAN 

SD 

MSG 



B50 MB WIND SPEED (kts) 



MEAN 

SD 

MSG 



PREGIPITABLE WATER (mm) 



MIAN 

SD 

MSG 



GE 



SE 



6.2 


5.0 


6.5 


9.2 


9.0 


8,0 


1 


t 


t 



N¥ 



671 


651 


689 


784 


548 


472 


441 


501 


1 


A 


I 


15 



8.4 


7.4 


8.7 


8.2 


t.l 


9.2 


8.6 


9.1 


a. 


i 


S 


2 



ftvJk 



9,3 


8.6 


8.0 


12.7 


4.9 


4.2 


4.4 


4.2 


t 


S 


« 


1. 



4.3 


3.3 


4.7 


5.2 


9.7 


9.5 


8.1 


9.0 


15 


i8 


m 


12 



24 


24 


iJ 


19 


12 


11 


ii 


i 


19 


22 


ii 


18 



22.4 


19.9 


22.5 


20.4 


12.9 


11.3 


11.2 


11,0 


m 


25 


41 


12 
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sw 



CE 



DIJEATION OF PRECIPITATION (hr) 



SE 



m 



TOTAL TOTALS 










MEAN 


17.4 


15.8 


20.7 


17.6 


SD 


18.7 


15.0 


12, § 


14.7 


ESQ 


18 


23 


17 


12 



HEAN 

SB 
MdG 

PRECIPITATION RATE (mm/hr) 



7.8 


7.8 


5.0 


5.4 


42 


70 



Fit 
1,1 

ii 



8.1 

6,3 

77 



MEAN 

SD 

MSG 



1.1 


1.2 


.1.2 


1.0 


1.0 


1,5 


1,2 


1,0 


42 


71 


51 


77 



DAILY AMOUNT (mm) 



MEAN 

SD 

MSG 



7.3 


6.7 


5.8 


5.2 


§ 






¥'■1 



7.1 

6.1 

7 



% OF STNS WITH PRECIP 



MEAN 
SD 



0.76 
0.22 



0.78 
0,21 



0.7B 
0.22 



0.54 
0.27 



SO, (mg/l) 

MEAN 

SD 

MSG 

NO3 (mg/l) 

MEAN 

SD 

MSG 



3.30 
1.8B 



2.88 
3.02 



3.05 
2.88 



1,14 
1.06 

1 



0.61 


0,59 


0.57 


0.24 


0-44 


0.72 


0.64 


0.22 


%. 


t 


i 


1 



APIOS PRECIPITATION (mm) 



MEAN 

SD 

MSG 



8.4 
8,3 

i 



7.5 
9.2 



8.7 

8.8 
3 



ill 



|*i 



Appendix H 



Table Hi i Sulfate and Iltrate data as a function of the categories of tbe discrete 
variables for 1983. 



SO, 



NO, 



MlAI 



SD 



MEAN 



Si 



SYNOPTIC TYPE 



WARM FRONT 


4.06 


1.78 


0.69 


0,49 


COLD FRONT 


2,98 


1.45 


0.50 


0.34 


STATIONARY FRONT 


2,55 


1.76 


0.44 


0.29 


OCCLUDED LOW 


3.05 


1.84 


0.59 


0.43 


COLD LOW 


4.23 


1.70 


0.66 


0.31 


AIRMASS 


2.80 


2.43 


0.65 


0.59 


WARM AND COLD FRONT 


4.50 


1.38 


0,78 


0.28 


AIRMASS 










CA 


1,90 


2,06 


0.61 


0.53 


CMA 


3.38 


2.32 


0.88 


0.81 


MA 


2.93 


1.70 


0.52 


0.42 


MP 


3.91 


1.93 


0.64 


0.47 


MT 


5.47 


3.22 


0.64 


0.23 


CA CMA 


1.10 


0.14 


0.49 


0.00 


CA MA 


2.76 


1.49 


0.62 


0.49 


CMA MA 


3.74 


2.11 


0.79 


0.5.7 


MA MP 


3.39 


1.71 


0.50 


0.33 


MP MT 


4.08 


1.50 


0.70 


0.39 


850 MB TRAJECTORY 










24 HR SECTOR 










NNW 


1.81 


1.71 


0.71 


0.65 


VNU 


2.45 


1,35 


0.52 


0.42 


WSW 


4.10 


1.92 


0.67 


0.45 


##P 


3.32 


1.61 


0.56 


0.39 


SSE 


1.97 


1.43 


0.34 


0.35 


ESE 


2.48 


1.72 


0.70 


0.77 


ENE 


1.74 


«• 


0.51 


_- 


NNl 


3.97 


4.85 


0.89 


1.24 
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so. 



NO, 



48 HR SECTOR 



ilAH 



SD 



MEAN 



SD 



24 HR RANGE (KM) 

<200 

200-400 

400-600 

600-800 

800-1000 

1000-1200 

1200-1400 



2.91 


1.45 


0.61 


0,48 


3.76 


1,99 


0.67 


0.48 


3.06 


2.01 


0.55 


0.41 


3.30 


1.98 


0.59 


0.45 


3.60 


1.99 


0.77 


0.61 


2.63 


1.09 


0.34 


0.17 



WS¥ 

SS£ 

£|S£ 

48 HR RANGE (KM) 



i.8f 


1.37 


0.60 


0.54 


2.60 


1.44 


0.53 


0.45 


3.90 


2.02 


0,64 


0.42 


3.35 


2-14 


0.66 


0.38 


3.47 


1.21 


0.87 


0.57 


0.54 


-. 


0.01 


-- 


7.40 


«.-.. 


1.77 


'Im-tm 



<200 


3 , 86 


0.36 


0.70 


0.54 


200-400 


2.73 


1.49 


0.56 


0.43 


400-600 


1.98 


1.49 


0.78 


0.76 


600-800 


4.00 


1.76 


0.69 


0,40 


800-1000 


3.87 


2.26 


0.72 


0.54 


1000-1200 


4.01 


2.06 


0.64 


0.49 


1200-1400 


2.37 


1.48 


0.42 


0.24 


1400-1600 


2.54 


2.06 


0.52 


0.28 


1600-1800 


2.95 


0.28 


0,49 


0.18 



SUIFACE TRAJECTORY 
24 HR SECTOR 



NNW 
VNV 

wsw 
ssw 

SSE 

ESE 
ENE 
NNE 



2.19 


2.76 


0.53 


0.70 


2.26 


1.46 


0.46 


0.46 


3.27 


1.40 


0.62 


0.44 


4.12 


1.80 


0.75 


0,43 


3.96 


2.02 


0.62 


0,30 


3.05 


1.7B 


0,49 


0.39 


2.43 


1.53 


0.49 


0,46 


2,51 


2.84 


0.69 


0.75 
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24 HR RANGE (KM) 



48 HR RANGE (KM) 



SO, NO, 

4 3 

MEAN SD MEAN SD 



<200 


3.24 


2.40 


0.60 


0.38 


200-400 


3.72 


1.88 


0.61 


0.45 


400-600 


3.91 


1.86 


0.63 


0.37 


600-800 


2.71 


1.71 


0.64 


0.57 


800-1000 


2.82 


2.07 


0.61 


0.54 


1000-1200 


2.41 


1.49 


0.46 


0.35 


1200-1400 


1.18 


0.82 


0,30 


0.08 


SURFACE TRAJECTORY 










48 HR SECTOR 










wm 


3-12 


2.73 


0.52 


0.34 


WNW 


2.41 


1.64 


0.57 


0.62 


WSW 


2.66 


1.28 


0.56 


0.39 


ssw 


3.89 


1.97 


0.66 


0.51 


SSE 


3.55 


1.93 


0.42 


0.27 


ESE 


3.29 


2.56 


0,73 


0.44 


ENE 


4.07 


1.37 


0.83 


0.46 


NNE 


3.32 


2,64 


0.75 


0.67 



<200 


3.93 


1.92 


0.69 


0.39 


200-400 


3.76 


2.63 


0.60 


0.44 


400-600 


3.35 


1,93 


0.63 


0.54 


600-800 


3.55 


2.06 


0.65 


0.44 


800-1000 


3.45 


1.33 


0.65 


0.42 


1000-1200 


2.71 


2.06 


0.63 


0.65 


1200-1400 


2.41 


2.41 


0.70 


0.70 


1400-1600 


_. 


-- 


.. 


-- 


1600-1800 


0.89 


0.27 


0.18 


0.20 
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so, NO, 

4 3 

MEAN SD MEAN SD 



TYPE OF PRECIPITATION 



CONTINUOUS 


2.83 


1.64 


Q.5f 


0.44 


CONVECTIVE 


3.66 


1.95 


0.65 


0.45 



NATURE OF PRECIPITATION 



RAIN 


3.69 


1.73 


0.60 


0.42 


SNOW 


1.91 


1.77 


0.65 


0.49 


FREEZING 


2.87 


1.09 


0,51 


0.39 


MIXED 


2.35 


1.94 


0.61 


0.61 



PREDOMINANT CLOUDS 



TOWERING CUMULUS 


2.38 


1.59 


0.46 


0.29 


ALTOCUMULUS 


2.67 


1.87 


0.33 


0.24 


ALTOSTRATUS 


6.58 


-- 


0.87 


.- 


CUMULONIMBUS 


3.81 


1.20 


0.52 


0.22 


CUMULUS 


3,34 


1.53 


0.66 


0.50 


NIMBOSTRATUS 


2.42 


1.19 


0.45 


0.28 


STRATOCUMULUS 


3.64 


2.21 


0.73 


0.53 


STRATUS 


2.95 


1.99 


0.65 


0.58 



SURFACE WIND SECTOR 

CALM 

NNW 
WNW 
WSW 
SSW 
SSE 
ESE 
ENE 
NNE 

850 Ml WIND SECTOR 

NNW 
WNW 
WSW 
SSW 
SSE 
ESE 
ENE 
NNE 



2,05 


-- 


0.45 


.-- 


2.03 


1.25 


0.51 


0.52 


2.91 


2.32 


0.54 


0.45 


3.14 


1.94 


0.64 


0,50 


3.79 


1.56 


0.63 


0.47 


4.11 


1.68 


0.57 


0.28 


3.62 


1.58 


0.69 


0.46 


2.80 


2.19 


0.53 


0.41 


2.07 


0.46 


0.43 


0.19 



3.01 


2.46 


0.65 


0.58 


2.33 


1.54 


0.52 


0.47 


4.03 


1.86 


0.70 


0.40 


3.81 


1.51 


0.69 


0.43 


2.68 


1.27 


0.47 


0.31 


2.63 


1.64 


0.38 


0.38 


2.02 


0.24 


0.23 


0.02 


3V00 


2.05 


0.80 


0.67 
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Table H.2: Correlation coef f IclentSj algnlf Icatice probabilities, and number of observations between SO. , M0-, and the 
continuous meteorological variables for 1983. 

CORIELATION COEFFICIENTS/PROB > R UNDER HOi RHO^O/NUMBER OF OBSERVATTONS 



SO, 



SO, 



NO, 



ffftf 
IP* 



CLOUD BASE HEIGHT 



SURFACE TIMPEIATURE 



SURFACE DEWPOINT 



SURFACE WIND SPEED 



850 MB TEMPERATUEl 



850 MB MIND SPEED 



0,229 


0.075 


0.01 


0.40 


130 


131 


0.266 


-0.049 


0.0019 


0.57 


134 


135 


0.250 


-0.058 


0.0036 


0.50 


134 


135 


-0.219 


-0.205 


0.01 


0.02 


134 


135 


0.328 


-0.021 


0.0003 


0.82 


120 


122 


-0.011 


-0.061 


0-91 


0.51 


116 


118 



PEECIPITABLE WATER 



TOTAL TOTALS INDEX 



DURATION OF PRICIP 



AVERAGE PRECIP RATI 



AMOUNT OF PIECIP 



AREA OF PRECIP 



0.232 


-0.068 


0.012 


0.46 


116 


118 


0.250 


-0.050 


0-01 


0.59 


117 


119 


0.147 


0.195 


0.16 


0.06 


'!# 


95 


0.206 


0.267 


0.50 


0.01 


94 


95 


-0.198 


-0.376 


0.021 


0.0001 


135 


136 


-0.105 


-0.261 


0.22 


0.0022 


135 


136 



Appendix I 

Table lit The frequency distribution of tbe categorical variables for the top 307. 
of the SO. and MO^ readings fox 1983, 



SO, 



SO, 



SYNOPTIC TYPE 



FCT 



NO 



PCT 



WARM FRONT 13 

COLD FRONT i 

STATIONARY FEONT 1 

OCCLUDED LOW 14 

COLD LOW i 

AIRMASS S 

WARM AND COLD FRONT 2 



32 
15 
2 
34 
5 
7 
S 



13 
6 
2 

11 
Z 
3 
2 



32 
15 

5 
27 

5 
12 

5 



AIRMASS 



CA 

CMA 
MA 
MP 
MT 

CA CMA 
CA MA 
CMA MA 
MA MP 
MP MT 



1 


2 


3 


7 


1 


2 


1 


2 


5 


12 


5 


12 


2 


5 


1 


2 


1 


2 


1 


2 








JO 





3 


7 


4 


10 


7 


17 


S 


20 


14 


34 


9 


22 


f 


17 


9 


22 



Hi 



850 MB TRAJECTORY 

24 HR SECTOR 

NNW 
UNU 
VSW 
SSW 
SSE 
ESE 
NNE 

24 HR RANGE (KM) 

<200 

200-400 

400-600 

600-800 

800-1000 

1000-1200 

1200-1400 

1400-1600 

48 HR SECTOR 

NNW 
WNW 
WSW 
SSW 
SSE 
NNE 

48 HR RANGE (KM) 

<200 

200-400 

400-600 

600-800 

800-1000 

1000-1200 

1200-1400 

1400-1600 

1600-1800 



SO4 NO3 



NO PCT it PCT 



1 


3 


2 


# 


3 


8 


5 


lA 


2 


61 


20 


m 


S 


22 


6 


17 








1 


3 


1 


3 


1 


3 


1 


3 


1 


■3 



i 





§ 





s 


14 


7 


li 


L4 


39 


13 


m 


9 


25 


7 


It 


3 


8 


4 


11 


#.. 


11 


5 


14 


1 


3 


i 














Q 



1 


4 


2 


7 


3 


12 


4 


m 


15 


60 


15 


m 


I 


12 


3 


li 


i 


8 


1 


f 


1 


4 


1 


k 



1 


4 


1 


4 


t 


i 


3 


11 


1 


4 


1 


7 


8 


32 


1 


21 


8 


32 


$ 


32 


3 


12 


3 


11 


■I 


4 


1 


A 


i 


4 


3 


11 





i 
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SURFACE TRAJECTORY 
24 HR SECTOR 

VNV 
WSW 
SSV 
SSE 
ESE 
ENS 
NNE 

24 HR RANGE (KM) 



SO^ NO3 



PGT NO PCT 



t 


3 


X 


3 


% 


3 


3 


9 


I 


9 


3 


9 


12 


35 


i3 


41 


11 


32 


5 


16 


2 


6 


1 


6 


3 


9 


3 


9 


1 


3 


2 


6 



<200 


3 


9 


4 


12 


200-400 


11 


32 


f 


28 


400-600 


12 


35 


7 


22 


600-800 


4 


12 


7 


22 


800-1000 


1 


9 


i 


12 


1000-1200 


1 


3 


1 


3 


1200-1400 


Q 





i 





SURFACE TRAJECTORY 










48 HR SECTOR 










NNW 


3 


10 


3 


11 


WNW 


2 


7 


4 


IS 


WSW 


1 


3 


2 


7 


SSW 


f 


31 


i 


m 


SSE 


3 


10 


1 


4 


ESE 


1 


10 


4 


15 


ENE 


4, 


14 


i 


ii 


NNE 


4 


14 


4 


it 


48 HR RANGE (KM) 










<200 


2 


1 


1 


? 


200-400 


4 


14 


t 


7 


400-600 


t 


24 


5 


If 


600-800 


t 


31 


It 


41 


800-1000 


3 


10 


t 


? 


1000-1200 


1 


7 


I 


7 


1200-1400 


2 


7 


1 


11 


1400-1600 
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so, 



NO, 



NO 



PCX 



PCT 



TYPE OF PRECIPITATION 



CONTINUOUS 
CONVECTIVE 



13 



32 

68 



16 

23 



m 



NATURE OF PRECIPITATION 



RAIN 
SNOW 

FREEZING 
MIXED 



36 


90 


29 


f4 


2 


5 


5 


13 








1. 


1 


2 


5 


3 


i 



PREDOMINANT CLOUDS 



TOWERING CUMULUS 


1 


2 


1 


1 


ALTOCUMULUS 


1 


1 


i 





ALTOSTRATUS 


1 


2 


1 


a 


CUMULONIMBUS 


8 


20 


5 


13 


CUMULUS 


3 


8 


4 


11 


NIMBOSTRATUS 


3 


8 


3 


i 


STRATOCUMULUS 


19 


48 


18 


4f 


STRATUS 


4 


10 


6 


i# 



SURFACE WIND SECTOR 



NNW 
WNW 
WSW 
SSW 
SSE 
ESE 
ENE 
NNE 



1 


2 


2 


5 


■4 


10 


4 


I.i 


7 


18 


1 


li 


f 


18 


1 


it 


3 


12 


4 


li 


1% 


28 


12 


w 


5 


12 


4 


m 











i 



850 MB WIND SECTOR 



NNW 

WNW 

WSW 
coy 

W' <J! ^1 

SSE 

JuiO' i-d 

ENE 
NNE 



2 


6 


3 


« 


3 


9 


3 


f 


13 


38 


14 


40 


12 


35 


10 


m 


1 


3 


2 


6 


1 


1 


X 


J 





® 





i 


2 


i 


i 


# 
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Table 1*2. Averages and standard deviations (SD) of the C0'ii,tin,uo'us 
ite 

'3- 



eteorological variables for the top 307. readings of the SO. and 
HO', 





SO4 




IO3 






MEAN 


SD 


MIAH 


SD 


CLOUD BASE HEIGHT 


874 


650 


745 


540 


SUIFACI TlMPERATirHE 


10.1 


8.1 


8.7 


9.4 


SURFACE DIWPOIIIT 


7.9 


8.1 


6.3 


9.2 


SMFAGE Wiro SPEED 


8,1 


4.3 


8.0 


4.4 


850 MB T1MPERAT011 


7,1 


8.4 


4.8 


10.3 


850 MB WIID SPEED 


25 


13 


22 


11 


PRECIPITABLE WATIR 


24.8 


12.7 


22.5 


14.6 


TOTAL TOTALS 


21.8 


22.1 


18.5 


22.2 


DURATION OF PCPN 


6.6 


3.3 


7,0 


3.7 


AVERAGE lAIMATE o,m/hr 


0,9 


0.9 


0.9 


0.8 


.AMOUHT OF PCPl 


5.8 


5,5 


5.0 


5.1 


I STN'S WITH PCPN 


0.74 


0.18 


0.70 


0.21 
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